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ABSTRACT 


As interest in high power free electron lasers (FELs) has increased, the FEL and ac¬ 
celerator communities have been faced with the need to develop high bunch charge, 
high repetition rate, low emittance electron sources for use as the driving accelerators 
for FELs. A novel superconducting, radio-frequency (SRF) gun/booster has been 
designed by and built for the Naval Postgraduate School (NPS) FEL Beam Physics 
Lab in collaboration with Niowave, Inc., for studying this electron source regime. 
The NPS SRF gun/booster operates at 500 MHz and is based upon a quarter-wave 
structure. It incorporates many features that make it desirable for studying the cath¬ 
odes and transport regimes necessary to explore high bunch charge beams, including 
adjustable field focusing, short transport out of the gun, and the ability to change 
cathode types and materials. After attaining “first beam” in June 2010, the NPS 
gun has been established as the first SRF electron gun in the United States. Initial 
results show excellent agreement with simulation with bunch charges of 110 pC and 
transverse emittance estimates of ~4 mm-mrad. Additionally, a modal analysis tool 
for the NPS FEL simulation software is developed based upon the Her mite-Gaussian 
basis set. Using a minimization of mode coefficients approach, we decompose output 
optical fields for amplifier FEL designs and experiments for FEL optimization and 
comparison of laser output fields. 
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EXECUTIVE SUMMARY 


Use of the free electron laser (FEL) has blossomed as an analysis tool for science, as 
well as developing research opportunities for possible high power applications. Cur¬ 
rent FELs are generally associated with large research institutions and are attached 
to electron accelerators and storage rings. As interest in high power has increased, 
the FEL and accelerator communities have been faced with the need to develop high 
bunch charge, high repetition rate, low omittance electron sources for use as the 
driving accelerators for the FELs. A novel superconducting, radio-frequency (SRF) 
gun/booster has been designed by and built for the Naval Postgraduate School (NPS) 
FEL Beam Physics Lab in collaboration with Niowave, Inc., for studying this electron 
source regime. This research supports the development of the NPS FEL and the U.S. 
Navy’s FEL Innovative Naval Prototype. 

The NPS SRF gun/booster operates at 500 MHz and is based upon a quarter- 
wave structure rather than the more prevalent elliptical cell design. It incorporates 
many features that make it desirable for studying the cathodes and transport regimes 
necessary to explore high bunch charge beams, including adjustable held focusing, 
short transport out of the gun, and the ability to change cathode types and materi¬ 
als. The cathode system allows for the use of conventional metal and semiconductor 
cathodes and the cavity generates fields of sufficient strength for testing field emission 
cathodes as well. With the cathode stalk removed, the cavity can be used as a booster 
cavity in an injector system with thermionic cathodes. Numerous simulations during 
the design phase are compared with experimental results and show excellent agree¬ 
ment. Experimentation and simulation by the author as part of the collaboration 
team is ongoing and results for bunch charge, beam energy, normalized emittance, 
and proposed improvements to date are presented and discussed. 

The development and construction of the prototype cavity is also an excellent 
case study for rapid prototyping in the SRF accelerator field as the gun/booster cavity 
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proceeded from initial design to first cold-test in less 12 months, shattering previous 
SRF gun development cycles by years. It is the first operational SRF electron gun in 
the United States. 

To develop a common comparison scheme for FELs and conventional lasers, a 
modal analysis tool for the NPS FEL simulation software is developed based upon the 
Hcrmite-Gaussian basis set. In amplifier FELs, there are no external factors deter¬ 
mining the appropriate transverse size and phase scaling factors to specify the basis 
set for modal decomposition. Using a minimization of mode coefficients approach, 
the author demonstrates the ability to decompose output optical fields for various 
amplifier FEL designs and experiments as tool for FEL optimization and comparison 
of laser output optical fields. 
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I. 


INTRODUCTION 


Free electron lasers (FEL) represent an active interdisciplinary research area 
incorporating challenges in materials, optical transport, electron sources, accelerator 
design, and many other areas of interest to the scientist and engineer. As the tech¬ 
nology associated with the FEL has matured, its usefulness as a laser that can be 
designed for many applications has increased demand for these systems worldwide in 
both academia and industry. 

While current user facilities are quite large and typically associated with large 
research universities or government user facilities, technology improvements are mak¬ 
ing many of the required subsystems smaller. There are applications in homeland de¬ 
fense, military settings, and other science fundamental and applied research programs 
that have an interest in smaller FEL designs. Particularly in military applications, 
there is also a desire to increase the output laser power of the FEL to increase its 
viability as a self-defense system. For many of these thrusts, necessary research in 
electron beam sources and linear accelerators is necessary. 

In this work, we report on the building and testing of a new superconducting 
radio frequency (SRF) electron gun/booster. The gun/booster is the first quarter 
wave structure designed to accelerate electrons and is unique in that it can be used 
as an energy booster, or, with the addition of a cathode stalk, as an electron gun. 
The gun/booster is designed to operate at 500 MHz, below the frequency of many 
current RF guns in use. Due to the quarter wave structure, it is still extremely small 
compared to its RF wavelength—the entire cryomodulc is less than ~1 m long. 

The full cavity development to “first beam” cycle was accomplished in less 
than 24 months. The cavity’s design and manufacturing progress has been presented 
in numerous poster and oral presentations as well in conference papers by the author 
and other collaborators. This approach is considerably different than the standard 
cavity design approach in which a cavity design is proposed and redesigned multiple 
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times before the cavity is ever fashioned, sometimes many years later. Our approach 
is based on rapid prototyping—understanding that the initial design will have flaws, 
but in the belief that experience will provide invaluable insight, enabling later cavity 
generations to be far superior than if only paper studies were conducted. Cavity 
performance during the initial testing has been below design criteria, but experimental 
results look promising for this type of cavity design. 

Since the end measure of effectiveness for an FEL is the output laser light, 
we develop a theoretical framework for modal analysis of the output light from an 
amplifier FEL. Based upon a Hermite-Gaussian mode decomposition of the optical 
held, an optimization of the free scaling parameters is performed to minimize the 
number of modes necessary to represent the optical held. By comparing the modal 
content of an FEL and a conventional laser, it is possible to perform a near “apples- 
to-apples” comparison of the expected propagation characteristics. 

This FEL analysis tool developed and coded by the author has been incorpo¬ 
rated into the Naval Postgraduate School (NPS) FEL4D simulation code and is used 
to analyze actual and proposed amplifier FEL designs as published by the author in 
collaboration with the the NPS FEL theory group [1]. Additionally, with planned im¬ 
provements to the FEL code allowing the importation of electron distributions from 
common particle accelerator codes, it becomes possible to couple the performance of 
new gun and accelerator designs with the FEL simulation and analyze directly the 
impact of these changes. 

In this dissertation, we will discuss briefly the cathode technologies that are 
compatible with current experiments using legacy equipment from the former Stanford 
FEL and for future use in the NPS SRF gun. The NPS gun/booster development will 
be traced and important design considerations discussed and analyzed. The details 
of the cavity construction process are documented for reference by future researchers. 
Beam simulation studies performed by the author for the NPS gun/booster in direct 
preparation for experimentation are presented, as well as a possible injector design 
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based upon the NPS SRF cavity configuration. Cold cavity and first beam tests 
developed by the author and collaboration team are conducted and presented showing 
that the gun design is successful and shows promise as a possible beam source for 
a linear accelerator. Finally, after a short FEL tutorial, we present the Hermite- 
Gaussian analysis tool and demonstrate its capabilities through FEL analysis. 
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II. 


ELECTRON SOURCES AND DRIVE LASER 


Generating an electron beam is essential in the operation of a free electron 
laser (FEL). The specific electron source can vary from machine to machine, as can 
the process for liberating those electrons from their source to be accelerated. In this 
chapter, we examine electron sources of interest to the research considered in this 
work as well in future work for the research group. We also describe the drive laser 
developed for cathode studies discussed later in this work. 

A. CATHODES OF INTEREST 

There are many ways to liberate electrons from a material to put into a beam. 
The three primary methods used in conjunction with linear accelerators are photoe¬ 
mission, thermionic emission and held emission. Some common desirable character¬ 
istics for a good electron source would be a low thermal omittance at launch, high 
current density, long lifetime, being robust, and providing a reliable method for beam 
shaping. In order of development, we will discuss the attributes of each type cathode 
and their applicability to our research. 

1. Thermionic Cathodes 

Thermionic cathodes have been in use since the first cathode ray tubes. In 
its simplest form, the thermionic cathode is a hot wire or plate in an electric held. 
The thermal energy imparted to the electrons provides them enough energy to escape 
the metal’s surface. Once free of the metal, the electric held accelerates the electrons 
toward an anode. If the anode is a mesh grid or annular in shape, the electrons can 
how through, creating a beam. A configuration with only a cathode and anode is 
referred to as a diode [2], A simple diode is shown in Figure 1. 

As the name implies, the thermionic cathode relies on thermal energy. If one 
desires to increase the current how from a thermionic cathode, all other variables held 
constant, one must increase the cathode temperature. Obviously, this cannot be done 
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Figure 1. A simple diode drawing. Electrons are liberated from a cathode by heating 
and accelerated toward an anode maintained at a higher potential. After [3]. 

indefinitely, as the cathode will eventually melt. Unfortunately, changing the cathode 
temperature is not a fast process, so attempting to generate a high frequency pulsed 
beam requires altering other variables in the system. One can vary the cathode- 
anode voltage, thereby changing the held seen by the electrons freed by their thermal 
energy. When the held is applied, beam is formed; when the held is oh, electrons 
are not pulled from the cathode and no beam is formed. Unfortunately for linear 
accelerators, the required cathode-anode voltage is typically in the hundreds of kV 
and switching high voltages at high speeds is not something to be undertaken lightly. 
There exists another, more tractable, method for generating a pulsed beam from a 
thermionic cathode. 

If we add another electrode, again meshed so the beam can pass through 
it, near the cathode surface, we can apply a much lower voltage and overcome the 
higher fields provided by the main anode. A diagram of this arrangement, called a 
triode, is shown in Figure 2. This type of thermionic cathode is what was installed 
in the Stanford free electron laser injector that is currently being used by the Naval 
Postgraduate School (NPS) for cathode studies. Since the voltage applied to the grid 
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is small (hundreds of volts or less), it can be switched quickly to provide sharp on-off 
edges for the beam as well as allowing for a high repetition rate. 


Anode 



Figure 2. A simple triode. Electrons are liberated from the cathode and accelerated 
toward a higher potential anode. Between the cathode and anode is a gate electrode 
that can be biased to enhance or stop electron flow across the cathode-anode gap. 
After [3]. 

As these cathodes have developed since the late 1800s, many more improve¬ 
ments have been made. As it was discovered that impurities can alter the electron 
affinity of the cathode surface, manufacturers began developing “dispenser” cathodes 
to purposely put beneficial impurities onto the surface, greatly increasing their sourc¬ 
ing capabilities. Electrode geometries have also improved as it was realized that the 
electric field lines can be used to shape the beam as it leaves the cathode, for ex¬ 
ample the Pierce geometry [4]. Additionally, photo-switching now allows thermionic 
cathodes to be switched at extremely high rates. In photo-switching, the cathode is 
used as a thermally assisted photoemitter, only emitting when a laser shines on the 
cathode surface. 

In general, thermionic cathodes are robust. For example, the cathode installed 
in the Stanford gun currently at NPS was the last operational cathode from the 
Stanford FEL. When the FEL was decomissioned at Stanford, the gun was let up to 
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air and left at atmosphere for two years prior to being put back into operation at NPS. 
While the performance may have been degraded compared to when the cathode was 
first installed, it easily conditioned back to high voltages and successfully generated 
continuous beam on the first attempt. The vacuum requirements for thermionic 
cathodes are not as stringent as those for semiconductor cathodes which increases 
the thermionic cathode’s robustness. Studies have shown thermionic cathodes can 
operate continuously for years without failure [5]. Although thermionic cathodes are 
not currently contemplated in the NPS superconducting, radio frequency (SRF) gun, 
the construction of the cavity does not preclude their use. 

2. Metal and Semiconductor Cathodes 

Metal and semiconductor cathodes operate using the photoelectric effect to 
liberate electrons from the cathode surface. Like the thermionic cathodes, these are 
immersed in an electric held to pull the electrons away from the surface as they are 
liberated. The electrodes or cavity shape are designed to pass the electron beam 
out of the gun and into the rest of the linear accelerator system. Unlike thermionic 
cathodes, a drive laser is required to provide the energy necessary to liberate the 
electrons. 

In the photoelectric effect, a photon is absorbed by an electron. If the absorbed 
energy is sufficient to match or exceed the work function of the material, the electron 
can be liberated from the cathode surface. The work function is a property of the 
material and varies from ~1 - 5 eV. The work function also dictates the maximum 
laser wavelength required for a drive laser to be effective. For example, let us consider 
a niobium cathode that has a work function of 4.3 eV [6]; metal cathodes in general, 
have a work function of ~4 eV. 

For a laser to liberate electrons from niobium, it must have a wavelength such 
that the photon energy is equal to or greater than the work function. The energy of a 
photon is given by E = hc/X, thus, ignoring two photon effects, our drive laser must 



have a wavelength shorter than 


he 

X = UeV = 288 nm ’ (IL1) 

which is in the ultraviolet range of the electromagnetic spectrum. Thus, using metal 
cathodes requires the use of short wavelength lasers. To use longer wavelength lasers, 
preferably in the visible range, it is necessary to move away from pure metal cathodes 
and consider semiconductor cathodes. 

Metal cathodes are nearly as robust as thermionic cathodes, requiring similar 
handling and vacuum requirements. Oxide formation and possible cathode prepara¬ 
tion requirements may make them less ideal for use. Their lifetime is long, but, in 
addition to the requirement for a UV drive laser, metal cathodes typically have poor 
quantum efficiency (QE), the number of electrons generated for each incident photon, 
often much less than 1%. 

Usually, a semiconductor cathode is not a bare semiconductor. They are nor¬ 
mally coated with a low work-function material from the alkali metals, such as cesium, 
which has a work function of 1.95 eV [6]. The addition of the alkali metal makes elec¬ 
tron generation much better with QEs up to 10%. However, the lifetime of these types 
of cathodes is short, requiring reprocessing of the cathode to restore QE. Usually, the 
processing is done outside the gun to prevent putting low work function material 
into the gun cavity where reflected laser light from the cathode could cause electron 
emission where emission is not desirable. An additional benefit of the lowered work 
function is that a visible laser can be used as the drive laser, ensuring that electron 
emission only takes place on the cathode, not the electrode surrounding it. Long life, 
high QE cathodes remain an area of active research at multiple facilities. 

Both types of cathodes are considered as candidates for the NPS SRF gun. 
The metal cathodes have already been demonstrated in both constant voltage (DC) 
and SRF systems. As the prevalence of alkali coated semiconductor cathodes in DC 
guns increases, there will be increased desire to try them in an SRF gun, especially if 
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QE lifetimes improve. At issue may be the higher vacuum requirement generally as¬ 
sociated with this type of cathode, and migration of reactive metals from the cathode 
into the cavity. 

3. Field Emission Cathodes 

Field emission sources have been explored in other applications, but are just 
beginning to show promise as electron sources for linear accelerators. To date, no 
held emitter cathodes have been installed in a realistic gun environment, either DC 
or SRF. Of the three types of cathodes discussed, these are the least mature. 

The cathode of greatest interest with respect to the NPS SRF gun is the dia¬ 
mond held emitter array (DFEA) being developed at Vanderbilt University. DFEAs 
are arrays of pyramidal structures that generate high current electron beams from 
their tips, shown in Figure 3. The thermal conductivity, mechanical strength, and 
potentially negative electron affinity of diamond make it ideal for this application [7]. 
Diamond held emitter tips have been seen to generate 100 /iA per tip with 32 MV/m 
held on the tip (turn on is about 15 MV/m) [7]. Data for these devices comes from 
a lab setting in which low voltages of about 1 kV are applied across small gaps (on 
the order of //rri) to generate the required turn-on helds. 

Anecdotal evidence from conferences indicate that these DFEA are extremely 
robust, having been left in drawers unprotected for up to six months and returning 
to full current when put back under held. If this can be quantified, it would bode 
well for many types of electron guns. When coupled with diamond’s resistance to 
contamination and oxide formation, the DFEA probably has the potential to be the 
most robust cathode candidate considered. 

The held emitting cathode is also an interesting consideration for RF guns, 
as the only mechanism necessary to generate an electron beam is an applied held 
of sufficient strength. While the general RF helds will be strong enough for beam 
generation, the peak of a sinusoidal RF wave is not necessarily ideal for transport. 
However, frequency mixing schemes have been proposed for generating fast turn on, 
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5. 0 kV X50.0k' ' 600 nil 


Figure 3. A scanning electron micrograph of a gated diamond field emitter. The 
pyramidal structure in the middle is the emitter with most of the electron emission 
being localized to the extremely sharp tip. 

short duration pulses from held emitter cathodes using RF guns [8]. DFEA cathodes 
are being considered as part of the NPS SRF gun experimental plan. 

B. DESIGN OF A UV DRIVE LASER SYSTEM 

As part of the establishment of the NPS FEL group’s laboratory spaces, it was 
necessary to develop a laser and transport system to conduct initial photoswitching 
experiments with the former Stanford DC gun. Additionally, the laser is to be used for 
separate test-stand experiments in the exploration of cathode performance in guns and 
gun-like structures. It is desirable to operate at multiple wavelengths to discriminate 
between emission materials and mechanisms where practicable. The system should 
also be transportable with minimal setup required before starting operation. 
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1. Minimalist Approach 

To minimize starting costs and development time, the initial drive laser system 
was developed using a minimalist approach. The drive laser chosen was the Contin¬ 
uum Minilite II laser. It is a Nd:YAG (neodymium doped yttrium-aluminum-garnet) 
that operates primarily at 1064 nm with nonlinear crystals to convert the wavelength 
to 532 nm, 355 nm or 266 nm at 50 mJ, 25 mJ, 8 mJ, and 4 mJ, respectively. The 
pulse length is about 4-5 ns with a low, 1 - 15 Hz repetition rate, commensurate 
with our current capabilities to collect data. It is certainly not the drive laser for a 
high-power FEL system, but it is more than sufficient for the current experimental 
needs. 

The initial beam line was designed to be attached to the Stanford DC gun and 
operate only at 266 nm to ensure photoemission from the metallic thermionic cathode 
surface. The optical path, shown in Figure 4, meets the bare minimum diagnostics 
and controls to be useful experimentally. There is one steering mirror for controlling 
the location of the laser spot on the cathode. The virtual cathode is set at the same 
distance from the steering mirror as the cathode, allowing the operator to track the 
laser spot size and motion via a video camera. Focusing of the laser spot is handled 
manually by inserting a focusing optic in a flip-in mount. For quantum efficiency 
measurements, an energy meter is used to measure the average pulse energy. 

2. Optimized Design 

As experience with the minimalist design was accumulated, the details devel¬ 
oped for the additional diagnostics that would help in quantifying the laser beam as 
compared to the electron beam that was emitted. Prior to commencing further, more 
detailed cathode studies and a design for an upgraded laser table to fit in the same 
footprint was developed. 

In the original design, to change the operating wavelength of the laser, it was 
necessary to open the laser head and change or remove crystals. The new table puts 
the doubling crystals on the table and intercepts the green laser light not upconverted 
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Figure 4. A minimalist drive laser design consisting of the bare minimum diagnostics 
and controls for performing QE experiments with the Stanford DC gun. The laser 
path is shown in yellow, terminating at the entrance to the vacuum system. 


to UV passing through the dichroic mirror to use as a second color laser. This allows 
for the green and UV lasers to be coaligned for easier diagnosing of the laser spot 
position inside the test stands and for two-color experiments. 

As the number of different experiments grew, the need to better quantify the 
laser pulse became more apparent. To add to the diagnostics of the laser beam, a 
photodiode for timing and temporal pulse shape has been installed. A beam prohler 
for measuring the transverse laser beam profile has also been installed. The need 
to frequently adjust beam spot sizes drove the decision to include a telescope for 
each beam line, allowing for spot size (diameter) adjustment from 4 mm to 50 /irn. 
Rotatable polarizers were installed, allowing us to adjust the output beam energy of 
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each beam line individually from 0% to 100%. The final beam transport line, shown 
in Figure 5, is much more complicated, but also more agile and capable for future 
experiments. 



Figure 5. The updated NPS drive laser table showing the incorporation of two-color 
operation and the addition of more diagnostics and controls to better characterize 
and control the laser beam. In this photo, not all items have been mounted to the 
table. 
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III. SRF GUN/BOOSTER 


Regardless of the electron source, we desire it to be located in a region of 
high accelerating held in order to bring the beam to high energy quickly. The longer 
the beam is left at lower energies, the greater the impact of space charge on the 
beam, increasing its emittance and making transport more difficult. This accelerating 
gradient can be applied in one of two ways—through a constant voltage or through 
the application of a periodic varying voltage. 

Constant voltage, or direct current (DC), electron guns are well developed 
and used in many electron sources today (e.g., [9], [10] and [11]). In these systems, 
a constant voltage is established between the cathode (usually at negative potential) 
and the anode (usually at ground) and possibly some intermediate electrodes for elec¬ 
tric held shaping. Figure 6 shows a drawing of the former Stanford Superconducting 
Accelerator electron gun as an example. Modern versions of these guns are typically 
designed to a standoff voltage of between 300 kV and 750 kV [9], [10]. They achieve 
fairly high helds at the cathode face, but the overall gradients are typically much lower 
than can be achieved with a varying accelerating held as applied by a radio-frequency 
(RF) held. 

The RF approach has a significant advantage over the DC case, as the time- 
averaged helds are much lower than the peak helds. As a result, the overall gradient 
can be increased, greatly reducing the time the electron beam is at low energy. Unfor¬ 
tunately, these improvements come with some signiheant costs. The beam now does 
not see a constant accelerating gradient, so timing of the injection becomes critical. 
Instead of a DC power supply, RF power supplies must be used, and waveguides must 
be used to port the RF power into a resonant cavity. The cavity itself must be tuned 
and designed to a specihc frequency. The fact that the accelerating voltages in an RF 
cavity oscillate mean that there are only certain times when the helds are conducive 
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Figure 6. This is a drawing of the former Stanford Superconducting Accelerator DC 
gun from files delivered with the gun. This gun operated with a thermionic cathode 
that was gated with an applied grid voltage. It routinely operated at voltages > 200 
kV and had a maximum gradient (with cathode at -240 kV and first electrode at 
ground) of 2.3 MV/m. It typically operated at some lower voltage, allowing the first 
electrode to provide beam focusing. 

to acceleration of charged particles, leading to pulsed beams. Higher gradients also 
necessitate higher wall currents, leading to heating of the cavity and parasitic power 
loss. 

In the RF accelerating scheme, there are two choices that can be explored in 
terms of cavity material type—normal conducting and superconducting. Both types 
of conductors provide the structure for resonant cavities; however, the power dissi¬ 
pated in normal conducting structures is substantially higher than that dissipated 
in superconducting ones [12], due to increased ohmic heating in normal conduct¬ 
ing cavities. While there are additional engineering and refrigerator requirements 
for superconducting cavities, they are more than compensated for by decreased RF 
power demands, especially in continuous wave applications. Table 1 shows the power 
comparison between copper and superconducting cavities for modest accelerating gra- 
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clients. We see immediately that, even for moderate gradients, the superconducting 
cavity provides us with significant power savings. 


Table 1. AC power requirements for superconducting and normal conducting cavities 
at 1 MV/m and 5 MV/m. Q 0 is the quality factor of the cavity, r a /Q 0 is the geometric 
shunt impedance, and E acc is the accelerating gradient. From [12], 


Option 

Super¬ 

conducting 

Normal 

Conducting 

Q o 

2 x 10 9 

2 x 10 4 

r a /Q 0 [Ohm/m], RF frequency = 500 MHz 

330 

900 

P/L [Watt/m] for E acc = 1 MV/m 

1.5 

56,000 

AC Power [kW/m] for E acc = 1 MV/m 

0.54 

112 

AC Power [kW/m] for E acc = 5 MV/m 

13.5 

2,800 


As of this writing, there is only one superconducting RF electron gun in use as 
an FEL electron source in the world. The superconducting RF electron gun used at 
the Forshungszentrum Dresden-Rossendorf in Rossendorf, Germany, shown in Figure 
7, is the electron source for the Free Electron Laser (FEL) at the Electron Linear 
accelerator with high Brilliance and low Emittance (ELBE). The prototype (a single 
half-cell) first developed beam in 2003 [13] and led to the production model gun, a 
3.5-ccll version, which made first beam in late 2007 [14], As this gun (referred to 
from this point as the “Rossendorf gun”) is the only superconducting RF (SRF) gun 
currently in use, its parameters are of particular interest for comparison purposes. 
Design simulation results of a similar geometry 650 MHz Rossendorf gun are shown 
in Table 2, which is a more relevant frequency for comparison. 

The Naval Postgraduate School (NPS) has established an FEL experimental 
program to complement its robust theoretical program, and desires to explore inno¬ 
vative accelerator, bcamlinc, and undulator designs to further FEL technology [17], 
[18]. To this end, it was proposed to develop and build a superconducting electron 
gun for use in our lab. As a brief overview, the gun is designed to operate at 500 MHz 
for frequency compatibility with the main linear accelerator. The gun should be a 
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Table 2. Characteristic design parameters of a 650 MHz Rossendorf SRF electron 
gun, from simulation. Parameters are taken from an unpublished report provided by 
Dr. J.W. Lewellen. The symbols for the parameters are defined in the list of symbols. 


Parameter 

Value 

fo 

650 MHz 

Qb 

1 nC 


24.6 ps, FWHM 

E 0 

26.48 MV/m 


1.19 mm-mrad 

K 

0.97 mm 

n 

9.3 ps 

KE 

13.8 MeV 

5 

0.23% 


* At the end of the simulated beamline, z = 1.42 m (from the cathode surface) 


compact, superconducting RF gun with high gradient and good beam quality. Having 
accomplished the design, development and construction of an initial prototype, NPS 
has the first SRF electron gun in the United States and is undergoing testing. The 
gun design, manufacture, and testing will be discussed extensively in the remainder 
of this chapter. 
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Figure 7. The Rossendorf gun is a superconducting photoelectron gun consisting of 

3.5 elliptical cells, operating at 2 K and 1.3 GHz [15]. The diagram also shows all the 
pertinent components for an SRF gun, including the cryostat. The cathode stalk is 
inserted from the left. For length scale, the 3.5 cell structure is approximately 40 cm 
from cathode surface to the end of the last cell, yielding a cryomodule approximately 

1.5 m in length and 80 cm in diameter. From [16]. 
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A. RF CAVITY TYPES 


There are many different types of cavities that can be used in RF accelerating 
structures. The cavity shape determines the accelerating field shape the electrons will 
experience in the gun. Changing the RF frequency, or introducing a second frequency 
into a cavity can also alter the cavity’s behavior. The ideal cavity mode for an electron 
gun is one in which there is a large axial electric field with small radial electric fields 
at an easily obtained RF frequency. The axial electric held requirement can be met 
by designing a cavity for a transverse magnetic mode (the TM 0 io mode, for example). 
The oscillatory nature of the RF fields requires that any charged particle beam will 
be pulsed in order to properly match with an accelerating held in the cavity. 

The hrst RF cavities were based upon the simple pillbox structure. The pillbox 
is a simple zeroth order cavity shape, allowing for simplicity in manufacture and 
modeling. This structure provides the jumping-off point for the much more advanced 
structures that are in use today. 

1. Right Cylinder Pillbox Cavity 

As mentioned, the simplest RF cavity shape is the pillbox structure. This is a 
right cylindrical structure with a cut-out in the end caps for beam entrance and exit 
as shown in Figure 8. The magnetic held has the desired transverse behavior and the 
electric held lies along the long axis of the cavity with the strongest helds on axis. 



Figure 8. A generalized pillbox cavity showing the transit direction of a particle 
beam. 
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As the pillbox cavity’s fields can be derived analytically, we include some of 
the results as they are illustrative for considering some of the more advanced cavity 
shapes. We follow the general process for solving for waveguide modes as can be 
found in several electromagnetism textbooks, [19] and [20], for example. We assume 
a plane wave of the form 

E (x, t) — E (r, 0) S kz ~ ut) , (III.l) 


where r is the radial distance from the propagation (z) axis, 0 is the angular distance 
from the x — z plane, k is the wavenumber, and u is the angular frequency. We also 
desire to separate the transverse effects from the longitudinal effects, so we define an 
operator, Vj_ = (d/dx,d/dy,0), which captures the transverse effects such that 


d 2 


(IH-2) 


We apply the wave equation to this field, 


1 0 2 


V 2 - 

c 2 d t 2 

2 d 2 1 d 2 

Vi H- 

dz 2 c 2 dt 2 


E 

E 


V 2 + ^ ) ES kz ~^ 


dz 2 


dt 2 


vi + ( ^ -1 


E = 


0 

0 

0 

0 . 


(HI-3) 


Since the wave equation applies to both the E and H fields and we can define 
them in a similar fashion (for a plane wave propagating in the z direction), then E 
in equation (III.3) can be replaced with H to describe the magnetic field. 

Since our expected cavity material can be considered a nearly perfect conduc¬ 
tor, our boundary conditions are the usual boundary conditions for a waveguide, 


E z \ s = 0 and H ± 0 = 0, 
where | s denotes “on the surface.” 


(HI.4) 
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Confining our discussion to the fundamental mode, and assuming we intend to 
use the cavity to accelerate beams, then the solutions we desire must have a non-zero 
E z . This confines our solutions to the transverse magnetic (TM) modes. The E and 
E[ holds in TM modes are related through 


H ± = ±-^~ (z x E x ) . (111.5) 

k 

The transverse electric fields are related to the accelerating held by 

E ± = ±p 2 V ± E z (III.6) 

where T 2 = (c o 2 /c 2 ) — k 2 and E~ satishes equation (III.3). These are traveling wave 
solutions, but we desire standing wave solutions to be applicable to a cavity [12]. 

To make our solution applicable to a cavity instead of a waveguide, we im¬ 
pose the boundary conditions necessary at the end caps, E± = 0 at z = 0, a. This 
constrains the z components, allowing us to find the accelerating held, 


E z (; x , t) — E (r, </>) cos e~ iuJt , (III.7) 

where the boundary conditions imply only certain wavelengths can satisfy the wave 
equation and be present in the cavity (k = pir/a) and p G {I > 0}, where J is the set 
of integers. Thus 


T 2 = 

mnp 




jury 

a J 


(III.8) 


where n G {J > 0}. Putting these together into equation (III.6), it can be shown that 


E ± = sin (^) V ± E (r, 0) e"^ . (III.9) 

ar mnp V O J 

For TM modes, we know H z = 0. However, we still need to find H±. Using 
equation (III. 5) and our E helds, we find 

Hi = cos (^) z x V ± E (r, 0) , (III.10) 

V C *- mnp \ CL y 

where r/ is the impedance of free space {pPo/co). 
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Finally, stipulating that our pillbox cavity has radius R, then the lowest order 
mode (TMqio) for our pillbox cavity has held components 

E r , E(j>i H T f H z 
E z 

H* 

where J n is the n-th order Bessel function. The 2.405 comes from the first zero of the 
Bessel function. 

In general, the radial electric fields can be related to the longitudinal electric 
held by using V ■ E = 0, to see how the radial helds behave as E~ changes along the 
cavity. This is true also for more complicated cavity shapes. We assume here that 
the helds are axisymmetric (E^ = 0), 
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We take the derivative of equation (III.11) with respect to r, 


d_ 

dr 


-E r 
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-1 1 dE r 
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(III.11) 


ignoring terms involving second derivatives because we assume we are near the electric 
axis of the cavity where the curvature should be small, 


1 fdE r 
r y dr 



dE r E r 

_ _ 

dr r 


(III.12) 
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Substituting equation (III. 12) into equation (III. 11), we find 


— E r H— E r + 


d E z 
dz 


E r 


—r d E z 
~Y~dz 


(III.13) 


which can be extremely useful in estimating the radial fields for cavities for which the 
axial fields are known. 


2. Tesla Cavity 

The first historical modifications to the pillbox structure for SRF applications 
were to round the corners where the end caps meet the longitudinal structure to 
improve surface current conduction by removing the 90-degree turn (and weld) and to 
improve the ability to clean the cavities. Removing surface contaminants is of utmost 
importance in superconducting cavity construction, as will be discussed later in the 
subsection on cavity cleaning and assembly. Cavity design modifications continued 
to improve manufacturability and performance. 

Probably the largest driver in early alteration of cavity cell designs is multi- 
pacting, an effect where field-emitted electrons are ejected into the cavity fields and 
then impact the cavity wall very near their origin, liberating secondary electrons that 
then are ejected into similar paths cascading into large populations that pull energy 
from the cavity fields. Once enough electrons are liberated, the cavity’s accelerating 
gradient is limited. Figure 9 shows a trajectory trace from a multipacting simulation. 
The first successful attempt to limit multipacting was to curve the corners of the 
pillbox, as in the HEPL (High Energy Physics Laboratory) accelerator at Stanford 
University [21], These cavities can be seen in Figure 10. 
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Figure 9. Single point multipacting trajectories in a muffin-tin cavity [22], This 
diagram demonstrates how a single electron can, over multiple RF cycles impact very 
near its origin causing localized heating, and for sufficient impact energy, liberate 
a significant population of secondary electrons. These secondaries will follow very 
similar paths, compounding the heating and cavity power limiting effects. 
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Figure 10. A seven cell HEPL structure. From [12]. 


Accelerator cavities continued to evolve, eventually being modified to an ellip¬ 
tical shape, as seen in Figure 11, allowing the cavity to be manufactured from two cell 
halves, with only one equatorial weld. Additionally, the continuously curved cavity 
walls also significantly reduce multipacting as the liberated electrons are forced to 
migrate to the cell equator, as shown in Figure 12, where the electric field nearly 
vanishes. The Tesla cavity geometry, developed for the TESLA linear collider, has 
become the predominant cavity type used in most modern superconducting electron 
accelerators, usually operating at frequencies above 1 GHz and gradients at or above 
10 MV/m. 

3. Quarter-Wave Cavity 

The cavities that have been discussed to this point have all had a hidden 
assumption built into their structures. That assumption is that the cavities are 
designed to accelerate speed-of-light particles (/3 ~ 1), where (3 = v/c. The structures 
to accelerate low /3 particles can have exotic shapes to provide for optimal acceleration. 
Of interest to us is the so-called quarter-wave structure, an example of which is shown 
in Figure 13. The components of the quarter-wave structure are a pillbox type outer 
skin and a coaxial line that intrudes into the interior of the cavity. 



Figure 11. A pair of 1500-MHz 5-cell, niobium cavities. The structures in the middle 
are RF waveguides. From [12]. 

As seen in Figure 13, this particular cavity is designed so that a beam is passed 
not through the end caps of the cavity, but through the side of the cavity and the end 
of the coaxial line. The coaxial structure allows for a TEM mode to be established 
along its length, providing an alternating accelerating gradient (a radial electric field 
in the cavity frame) across the two gaps the beam would see. While the accelerating 
gradient is reversing, the beam is shielded in a drift tube until the gradient again 
reverses and provides a second acceleration as the beam exits the cavity. Driving 
frequencies in quarter-wave structures are typically much lower than in (3 — 1 cavities, 
which decreases the thermal load on the cavity due to the alternating currents induced 
in the cavity walls. 

We propose to use the quarter-wave structure for the acceleration of electrons 
and the development of a compact electron beam source. Our usage of the structure 
will not be in its usual mode; instead we will establish the accelerating gap between 
the end of the coaxial line and an exit port through the opposite endcap. This gives us 
a cavity that should be easy to manufacture, will operate at a frequency low enough 
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Figure 12. The typical mutipacting trajectory in an elliptical cavity is shown on the 
left. On the right is a different form of resonant multipacting that is not eliminated 
by the use of the elliptical shape. After [12], 


to operate at 4.2 K, and can be used as both a booster and, with the addition of 
a cathode stalk, as an electron gun. The remainder of this chapter will address the 
major components’ manufacture, modeling, and testing. 
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Figure 13. The diagram on the left shows the generalized quarter-wave structure. On 
the right is a drawing of a quarter-wave structure used for the acceleration of heavy 
ions. The central coaxial line in this device is filled with liquid helium for cooling. 
After [12]. 


29 




B. FABRICATION 

The NPS quarter-wave, superconducting cavity was initially conceived as a 
booster cavity to be nsed in conjunction with the former Stanford DC electron gun. 
Together, the two devices would make up the injector section of the NPS FEL, pro¬ 
viding electrons nominally at 1.5 MeV (~250 keV from the gun and 1.2 MeV from the 
booster). The booster cavity is designed to handle significant bunch charge (1 nC) 
at fairly low repetition rate (initially 10s of Hz building to kHz as the experiments 
progress) to minimize the liquid helium cooling and RF power required. To also de¬ 
crease the cryogenic cooling and auxiliary support equipment necessary, the cavity is 
designed to be operated at 4.2 K, which is above the superfluid helium temperature 
of 2 K, but sufficiently cold to cause niobium to become superconducting (T c = 9.2 
K). The fundamental frequency was chosen to be 500 MHz in order to be compatible 
with the linear accelerator to which it is to connect. In this case, the accelerator 
operates at 1.3 GHz, meaning that this system’s highest possible repetition rate (lim¬ 
ited by when the various frequency RF fields throughout the beamline are properly 
in phase) is 100 MHz. Since a bunch can only pass through the entire line (due to 
the required phase alignment) once every 13 cycles, this limits the maximum duty 
factor of the accelerator sections to about 8%, reducing RF and cooling power. As 
the NPS FEL is designed to be a low duty factor system (on the order of Hz to kHz 
for the foreseeable future), the 100 MHz bunch frequency does not place a significant 
constraint on system operation. For reference in the remainder of this section, the 
final NPS 500 MHz superconducting gun/booster schematic is shown in Figure 14. 
In this layout, electrons are injected from the left and are accelerated as they pass 
through the nosecone and out the cavity to the right. 
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Figure 14. The NPS 500 MHz superconducting gun/booster as built by Niowave, Inc. 
Shown in the drawing are the ports for flowing liquid nitrogen and helium as well as 
the various enclosures inside the cryostat. Not shown are the locations of various 
internal temperature sensors, the cathode stalk assembly, or coupler assembly. 
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1. Cavity 

Initial cavity design work was performed using the Superfish 1 family of codes. 
The cavity structure was optimized to obtain the desired fundamental frequency 
and held strength by altering the cavity dimensions, including coaxial line diameter, 
overall cavity diameter, cavity length, and coaxial line length. To explore how the 
coaxial parameters affect the fundamental frequency, a right-cylindrical pillbox cavity 
with the same length and radial dimensions of the intended NPS gun, as shown in 
Figure 15, is examined as a right-cylindrical coaxial line was inserted from 0% to 95% 
of the cavity length. Figure 16 illustrates how these parameters affect the cavity’s 
fundamental frequency. As can be seen in Figure 16, the coaxial line radius has 
significant impact on the cavity frequency, depressing it until the coaxial line radius 
reaches half of the cavity radius. Increasing the coaxial line radius beyond 50% of the 
cavity radius only serves to drive the fundamental frequency back up. Insertion of 
the coaxial line further into the cavity depresses the cavity frequency almost linearly 
once inserted beyond 25% of the cavity length. By varying both parameters, a desired 
frequency and accelerating gradient and gap can be designed for a particular cavity. 

While performing this design exercise for the production cavity, we also en¬ 
sured that the fields expected in the cavity did not exceed fundamental limits on the 
maximum surface magnetic held. For niobium, the peak surface magnetic held must 
be less than the superheating held, 0.23 T (183 kA/rn). This H p k limits the peak 
electric held that can be obtained in the cavity to approximately 55 MV/m as well 
[12], Based on the quarter-wave structure, we expect that the peak magnetic held 
will occur in the coaxial line section of the cavity while the maximum surface electric 
held (which determines the amount of held emission we expect) should occur along 
the nose cone. 

After optimization, we hnd via Superfish held maps of the booster configura¬ 
tion that, for the nominal 1.2 MV gap voltage, our maximum magnetic held is on 

1 The Superfish family of simulation codes will be discussed in Chapter IV.A.l.a 
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500MHz General Cavity Search, 30% Indentation F = 707.21996 MHz 



Figure 15. An example pillbox used to explore fundamental frequency shifts due to 
insertion of a coaxial line. Coaxial lines are inserted from the right and insertion 
lengths and coaxial line radius are varied. The pillbox structure is 20.305 cm long 
and has radius 12.000 cm. Pink lines indicate electric held lines and red arrows show 
relative electric held strength. 


the order of 65 mT occuring approximately two-thirds of the way back from the nose 
cone along the coaxial line, while the maximum electric held of 49 MV/m occurs just 
inside the tip of the nose cone. Also of interest is the held that will be experienced at 
the cathode surface in the gun conhguration. Since the cathode stalk can be adjusted 
to the compression and extension limits of the bellows, we can adjust the cathode’s 
surface from a retracted position 14.23 mm behind the nose cone to about 1.51 mm 
beyond the nose cone. Since extending the cathode into the accelerating gap results 
in the emitted electrons experiencing significant defocusing helds upon emission, we 
limit onr consideration of positions to those from insertion hush with the nose cone 
to retraction still sufficient for held emission cathodes (\E\ > 10-15 MV/m) as 
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Figure 16. Fundamental frequency effects in a pillbox cavity seen when varying the 
coaxial line radius and penetration depth (expressed as a percentage of the pillbox 
length) into a simple pillbox structure. The pillbox structure is 20.305 cm long and 
has radius 12.000 cm. 


discussed in the previous chapter. Figures 17 and 18 show the accelerating held, E z 
(on axis), and transverse held, E r (r = 1 cm), for various cathode positions and the 
booster configuration. 

We see from Figures 17 and 18 that there is an obvious trade-off between peak 
accelerating helds on the cathode surface and radial focusing as the beam is launched 
from the cathode surface. From an applications standpoint, assuming the cavity 
achieves its designed gradient, all cathode positions forward of 6.5 mm retraction 
would be usable for held emitter cathodes which require approximately 10 MV/m for 
turn-on [23], [7]. 
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Figure 17. Plots of the axial electric field, E z , along (r, 9, z ) = (0, 0, z ) as the cathode 
position is changed within the cavity. Retraction distances are measured relative to 
the cavity nose cone (z = 0). Data is from Superfish simulations. 


Having a radial focusing field at the cathode surface is beneficial for many 
reasons. Most important is countering the space-charge effects tending to increase 
the beam’s size and emittance until the beam can be accelerated significantly. Figure 
18 provides radial field strength curves for various configurations of the NPS quarter- 
wave cavity. For cathode retractions less than 3.0 mm, there is no focusing as the beam 
leaves the cathode. Maintaining a retraction beyond this point provides focusing for 
the beam as it leaves the cathode until it can be accelerated, countering the effects of 
space charge, but at the expense of accelerating field, E z . This indicates the possibility 
that this gun configuration could be useful for high bunch charges. Also interesting is 
that when used in the booster configuration, there are significant focusing fields (> 2 
MV/m) as the beam passes through the cavity. In many injectors, acceleration after 
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Figure 18. Plots of the radial electric field, E r , along ( r,9,z ) = (0.5 cm, 0, z) as the 
cathode position is changed within the cavity. Retraction distances are measured 
relative to the cavity nose cone (z = 0). Negative E r indicates a radially focusing 
held. Data is from Superfish simulations. 


the gun is delayed until after a focusing element in order to preserve bunch charge 
prior to further acceleration. It may be possible to eliminate that focusing element 
and decrease the size of an injector using this cavity as a booster. 

Actual manufacture of the cavity was performed at Niowave, Inc., facilities 
in Lansing, MI. Its facility includes almost all machines, materials, and facilities for 
manufacture and testing of accelerating structures. With the exception of electron 
beam welding of the cavity, all manufacture and processing of the NPS cavity was 
performed on-site at Niowave. 
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The cavity began as the four primary components shown in Figure 19. The 
nose cone cylinder and end plates were cut from a high residual resistivity ratio 
(RRR), large grain niobium ingot. The hollow cylinder was formed by bending small 
grain sheet niobium. Higher values of RRR material are desirable, as they result 
in decreased heat load generated through the high alternating currents in the cav¬ 
ity surfaces when operating and increased heat transport (thermal conductivity) of 
the generated heat. The ingot material used had a RRR value of 180, while the 
sheet material had a RRR value of 300. These are representative of “typical” values 
of commercially available niobium and no significant attempt was made to procure 
extremely pure niobium for this cavity. 



Figure 19. The component parts of the NPS 500 MHz quarter-wave cavity. The two 
large disks at the bottom are slices from a large grain niobium billet and form the 
end plates of the cavity. The cylinder on the left is rolled sheet material that makes 
up the outer conductor. The solid cylinder on the right is cut from the same billet 
as the end plates and makes up the inner conductor. The two remaining parts in the 
picture are niobium-titanium flanges that mate the cavity to the beamline. 
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The two end plates and nose cone were milled to final thickness using standard 
machining processes. The nose cone is hollowed between the inner beam pipe and 
outer coaxial surface to allow liquid helium to reach as close as possible to the high 
current areas of the cavity. 

Once the cavity was assembled, the first set of verification experiments were 
performed. These experimental results are presented and discussed in the simulation 
and experiment chapter. 2 

2. Coupler 

The primary function of the coupler is to transport RF energy from a source 
into the stored fields of the cavity. In the case of many cavities, this is accomplished 
either through waveguides or RF antennae. As our cavity is extremely compact and 
our design philosophy included minimizing the number of cavity penetrations, it was 
decided not to use a cavity penetrating RF antenna for the first prototype. 

Our design choice was to use an axial waveguide coupler. This has the benefit 
of not requiring an additional cavity penetration as the output beam tube is used 
to convey the evanescent wave from the end of the coaxial waveguide to the cavity. 
Additionally, the hollow coaxial line allows for electron beam extraction. The coupler 
drawing is shown in Figure 20. The two primary drawbacks to this design are that 
the coupler decreases the available aperture for extracting the electron beam and 
the required drift space occupied by the coaxial coupler allows space charge time to 
expand the beam. Simulations described later indicate the resulting aperture decrease 
is not a major issue except at lower cavity gradients. 3 However, the increased drift 
space required by this design and other factors may indicate that a different form 
factor may provide substantial benefits in this cavity’s use as a booster or gun. 

The configuration shown in Figure 20 has two external connections. The 
connection on the bottom is an RF connection for coupling RF power onto the coaxial 

2 See Chapter IV. 

3 See Chapters IV.A.2., IV.A.3. and IV.B.6. for explanation and simulation/experimental results. 
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Figure 20. The NPS cavity RF coupler drawing. The pink portion on the left is the 
500 MHz SRF cavity. The superconducting solenoid is shown near the coupler tip. 
The vertical structures are input and output antennae for exciting the waveguide and 
for reading the return RF signal from the cavity. 


antenna. This is done through capacitive coupling. The connection directly opposite 
of the RF input connection is a receiving antenna, also capacitively coupled. This 
antenna receives two signals: one is a direct path transmitted signal from the input 
antenna and the other is a signal from the cavity that is proportional to the power 
contained in the cavity. As the cavity fields “fill,” this return signal is small. Once 
the cavity has reached its maximum capacity, the return signal reaches a steady state 
level. 

This coupler configuration also allows for adjustment of the coupling constant 
to the SRF cavity. By compressing/extending the bellows around the coupler attach¬ 
ment to the beam tube wall, the coupler tip location relative to the cavity can be 
adjusted. The coupling constant is a measure of the overlap between the evanescent 
coupler fields and the cavity fields [12]. Ideally, the coupler is adjusted such that 
“unity coupling” is achieved. Unity coupling indicates that the cavity is acting as a 
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matched load to the coupler, allowing for efficient transfer of power from the coupler 
to the cavity. To measure the coupling constant, one can observe the behavior of the 
transmitted and return power levels of a square wave gated RF pulse on a frequency 
analyzer. Figure 21 shows illustrative plots for this concept. 


0 < 0.5 0 = 1 



0 > 2 



Figure 21. The left column shows weak coupling. Pf is the RF feed power, P t is 
a measure of the power contained in the cavity, and P r is the power reflected back 
to the feed, but also contains the power emitted from the cavity. When 0 — 1, we 
note that the emitted power and reflected power exactly cancel, while for 0 > 1, the 
emitted power dominates and for 0 < 1 the reflected power dominates P r . The first 
peak in P r is associated with the ratio of reflected to forward power and the second 
peak is associated with the ratio of emitted to forward power. After [12], 


The coaxial coupler antenna is detuned from the fundamental frequency of the 
cavity, 500 MHz, so that it acts as a transmission line rather than a resonator. The 
coupler has a length of 60.5 cm corresponding to resonant frequencies of 123.9 MHz 
(A/4), 371.6 MHz (3A/4), and 619.4 MHz (5A/4). Driving the antenna at 500 MHz 
allows for RF coupling to the cavity without the antenna being resonant. The input 
and sensing antennae are located away from the shorted end of the transmission line 
in order to enable driving the RF amplitude on the coupler as well as allowing for 
reception of the return signal on the output antenna. However, with the required drift 
length of about 75 cm from the cavity (this length includes the bellows section after 
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the coupler termination), there is significant opportunity for electron beam growth 
and omittance degradation at low energies. While the superconducting solenoid pro¬ 
vides focusing, from a beam line layout perspective, it would be advantageous to 
position a booster cavity as close as possible to the gun. In the booster configuration, 
the level of desirability of having another booster cavity located close to this one is 
dependent on the beam energy prior to entering the 500 MHz booster. 

In actual practice, this coupler arrangement turns out to be extremely sensitive 
and difficult to diagnose. In the simpler booster configuration the coupler arrange¬ 
ment is really three coupled systems: the input antenna, the coaxial transmission 
line, and the SRF cavity. An adjustment to one of the antennas does not have easily 
predictable results on the entire system and the arrangement is not easily modeled 
for predictive simulation purposes. In the experiments reported later, the input an¬ 
tenna was inserted in as far as possible without shorting to the coaxial line and all 
antenna adjustments (to adjust (3) were made by changing the compression of the 
coaxial coupler termination bellows. 

In the next iteration of this cavity, it is highly recommended that an intracavity 
loop antenna or transverse coaxial coupler be used, an example of which is shown in 
Figure 22. A rectangular waveguide at the drive frequency desired would be quite 
large relative to the small diameter pipe connecting to the cavity, making it less 
practical. An intracavity loop antenna requires an access port through the cavity 
wall requiring a vacuum tight weld. While introducing a penetration complicates 
the construction of the cavity and cryomodule and introduces a possible increase 
in the number of surface defects and their associated effects on cavity performance, 
the increase in control of the RF coupling and decrease in beamline length may 
compensate these drawbacks. The transverse coaxial coupler allows for greater RF 
power input than the antenna version while maintaining the cavity with only beam 
tube penetrations and achieving a much reduced beamline length over the current 
design. 
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Figure 22. A transverse coaxial coupler utilizing a similar coupling mechanism (al¬ 
though waveguide fed in this case) that could be used to provide power to the NPS 
cavity while greatly reducing the drift length to the next beamline element. The 
antenna is the central vertical structure that projects slightly into the beam tube. 
From [12]. 
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3. Cathode 

To take advantage of the benefits of the cavity design, the cathode assembly 
required a significant amount of design work. As the gun is intended as an experimen¬ 
tal set-up for multiple cathode types and materials, the assembly must be removable 
to exchange cathodes for testing. In most guns, this process is handled by a load-lock 
system allowing for insertion/removal of a cathode without breaking vacuum in the 
gun or cathode assembly, [24] for example. As these load-lock systems are very ex¬ 
pensive and one-off designs, and the cavity cannot be tested in the gun configuration 
to verify simulation until some form of cathode is installed, it was determined that 
using a standard metal cathode was a suitable, low risk first step while the remainder 
of the cathode assembly interface requirements were fleshed out with other research 
groups. To prevent foreign metal contamination in the prototype cavity, we chose to 
use niobium as the cathode material. It has a work function similar to copper (Nb ~ 
4.3 eV, Cu ~ 4.7 eV [6]). In order to liberate electrons via the photoelectric effect, 
the incoming photon must have energy equivalent to the materials work function; 
for niobium, this requires a laser with a wavelength shorter than 288 nm. Niobium 
and copper also have similar quantum efficiencies, on the order of 10 -5 , depending 
on cathode preparation techniques used. As we would also like to test some of these 
cathode preparation techniques[25], [26], [27] and there exists the possibility for ab¬ 
lation of cathode material during cathode processing, the niobium cathode material 
choice presents the least damage risk to the cavity. Using another metal could result 
in plating of that material inside the cavity, thus, degrading the superconducting 
properties of the cavity. The use of an alkali metal on a semiconductor substrate 
requires better vacuum levels than currently available on the system. Alkali metal 
cathodes also present issues (possible contamination of the cavity) that we wish to 
avoid until the cavity has been properly characterized. 

Once it was decided that the cathode would be metal, the next major decision 
was how the cathode should interact with the cavity electrically. One possibility 
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is to mechanically short the cathode stalk to the cavity, preventing the RF power 
from passing down the cathode stalk out of the cavity. Recognizing that one of the 
inherent benefits of the cavity’s geometry is the adjustable radial fields on the cathode 
surface, any simple RF short would inhibit the use of this characteristic, making this 
approach less desirable. Additionally, any short would require physical contact in the 
high held region of the cavity. Any metal particulate introduced here has a very high 
probability of generating a held emission site, degrading gun performance. The other 
possibility is a cantilevered stalk, not shorted directly to the cavity, to support the 
cathode in the desired position. The simplest non-shorted cathode stalk is a right 
cylinder inserted into the nose cone area. 

This simple design comprised the hrst iteration of the cathode stalk. It became 
apparent quickly that there was a singular haw with this design. Just like the in¬ 
put coupler, the cathode stalk becomes an RF transmission line allowing RF energy 
to how down it as a coaxial waveguide. A simple calculation can show how much 
energy this configuration will pull from the cavity and turn into resistive heating. 
This heat is a significant concern because the cathode stalk is not ideally connected 
to the outside environment for conductive cooling, and without active cooling, the 
cathode stalk could present a non-trivial heat load to the cryogenic system, possibly 
impacting the ability of the cavity to remain superconducting in the gun configura¬ 
tion. Additionally, any RF power pulled from the cavity decreases the cavity holds 
used for beam acceleration unless additional RF power is supplied (which would also 
increase the power lost to the cathode stalk). 

To model this system analytically, we simplify the cavity to a pillbox with a 
cathode set at the input side of the cavity. The cathode stalk extends inside a vacuum 
pipe to a short at the base where it connects to the beam pipe, as seen in Figure 23. 

The cavity helds cause a surface charge, sourced from the stalk, to build up 
on the cathode. We can determine this charge using Gauss’ law, forming a Gaussian 
pillbox around the cathode surface. If we place the pillbox boundary near the surface, 
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Figure 23. A resonant cavity has E z as shown. This held impinges on a cathode 
attached to a right cylindrical stalk. The cavity resonates at a fundamental frequency, 
/o- 


the non-z components can be neglected, leaving only Eqz. We calculate the surface 
charge 

(j E ■ n d A = — , 

Js e o 

where S is the surface of the Gaussian pillbox, A is the area of the cathode surface, h 
is the outward normal vector of the Gaussian pillbox, and q is the charge enclosed by 
the surface. Solving this equation, we find the charge that builds up on the cathode 
surface is 

q it) — E 0 e 0 A sin (ut) , (III.14) 


where we have included the oscillatory nature of the electric field with amplitude Eq 
and angular frequency u = 2i r/ 0 . If we take the derivative of (III.14), we get a current 


d q 
dt 


E 0 €qAu cos (ut) = Iq cos (cot) , 


(III.15) 


where J 0 = E 0 e 0 Aui. This current is conductively attached to the cathode assembly, 
which acts as a coaxial transmission line relative to the cavity. The current that pop¬ 
ulates the cathode “sloshes” along exterior surface of the stalk, and we can state that 
the current density at each point along the stalk (assuming the stalk is approximately 
3A/4 relative to the cavity fundamental so that the stalk couples well to the cavity) 
is 

J stalk (A t) = cos (ut) cos f ^ j . (III. 16 ) 
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We integrate out the z dependence in preparation for finding the power, 
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(III.17) 


Since we desire to find the average power lost to the cathode stalk and we have 
the current already, we use the ohmic power equation 


(P) = k ■ 


(III.18) 


Starting with the resistance, we know 

n = , (hi. 19 ) 

where p is the resistivity of the material, l is the length of the conductive path, and 
A ,5 is the this annulus shaped cross-sectional area through which current is flowing. 
The cross-sectional area in this case is a small ring at the outer surface of the stalk 
whose depth is equal to the skin depth of the material, S, at the driving frequency, 
/ 0 . Thus the area is 

A s = nR 2 - 7T (R - 5f , (III.20) 


where R is the radius of the cathode stalk and 


5 [m] « 503.3y j- , 

where /o has units of [Hz] and p has units of [H m]. Substituting (III.20) into (III. 19), 
we obtain 

where l is assumed to be a longer length, A = c// 0 , for power dissipation purposes. 
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Returning to equation (III.18), we now have 
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(III.22) 


Choosing some representative numbers for the gun design: i? 2 = 20 MV/nr, / 0 = 500 
MHz, R — 1 cm, and room temperature copper as the stalk material, we find that 
the power dissipated in the stalk is 1.94 kW. Comparing this with a Superfish model 
of the actual NPS cavity with comparable parameters and a right cylindrical cathode 
stalk, we find 2.04 kW dissipated in the cathode stalk. 

Analyzing equation (III.22), we see that there are several parameters that can 
be used to control the power loss along the stalk. Reducing fields on the cathode 
surface, decreasing the operating frequency (which also increases the skin depth and 
further decreases the power), choosing a material with a lower resistivity, and decreas¬ 
ing the area of the cathode surface area (controlled by the stalk radius). Of these, the 
frequency of the cavity is fixed; however, the material choice for the cathode stalk, 
cathode surface area, and cathode surface fields are all adjustable within the confines 
of the cavity design. 

Additional losses can be attributed to RF power being coupled onto the cath¬ 
ode stalk as a transmission line. This problem is not a new one and has been ob¬ 
served/predicted in both the Rossendorf gun and the proposed Brookhaven SRF gun. 
Both of these guns were envisioned as guns from the outset and were designed to in¬ 
corporate an RF choke, essentially a resonant cavity that preferentially gathers the 
RF power and prevents it from flowing down the cathode stalk. Figure 24 shows the 
choke designs incorporated into these guns. 
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Figure 24. RF choke designs for two SRF gun designs. On the left is the Rossendorf 
gun and the choke cavity just behind the cathode surface. From [15]. On the right is 
the proposed BNL/AES SRF gun with a different choke design placed very near the 
cathode surface. After [11], 

Since the NPS design was initially contemplated as a booster, there is no RF 
choke incorporated in the superconducting section of the cavity. To reduce RF power 
flowing down the cathode stalk, our approach, proposed by our Niowave collaborators, 
has been to shape the cathode stalk so that the cathode stalk acts as a Bragg reflector, 
greatly reducing the RF power coupling down the cathode stalk. Figure 25 shows the 
final version of the cathode assembly. Note that it incorporates several features to 
reduce RF coupling. The cathode holder has a narrow fit into the nose cone to reflect 
as much RF power back into the cavity as possible. The multiple diameter changes 
along the support stalk act as impedance mismatches (or a Bragg reflector) to reflect 
any power that has passed down the stalk back toward the cavity. In addition to 
changing the stalk geometry from the simplified model, we explored the impact of 
changing the conducting material inside the beam tube and the application of active 
cooling. These strategies were modeled in Superfish and the results are summarized 
in Table 3. 

The geometry changes alone decrease the losses sufficiently to allow for oper¬ 
ation as a gun. Should further heat reduction become necessary, additional measures 
can be taken. Anodizing the interior of the vacuum pipe with a layer of copper and no 
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Bellows Cathode Assembly Cathode Tip 

Nosecone 

Figure 25. The assembly used for insertion and positioning of the cathode surface in 
the NPS SRF cavity gun configuration. 

further changes results in a reduction in ohmic heating to half of the original design. 
The thickness of this copper coating would need only be 10 /rm to achieve a thickness 
greater than 3 skin depths at 500 MHz. If active cooling of the cathode assembly is 
desired, making the tip end hollow and passing a liquid nitrogen feed would reduce 
the radiated power to the cavity significantly as well as decreasing the material’s re¬ 
sistivity. Active cooling could also be applied to the exterior of the vacuum pipe as 
well, although the effect would be smaller than cooling the cathode stalk. Applying 
all the improvement measures considered results in a reduction of ohmic losses to 
2.7% of the original design. 

One of the best attributes of this cathode/cavity combination is the ability 
to tune the fields on the cathode surface by adjusting the cathode’s position relative 
to the nose cone. The bellows assembly at the shorted end of the cathode assembly 
allows for the longitudinal positioning of the cathode surface from a position 1.51 
mm beyond the nose cone (with fields E Zt max = 38 MV/m and E r> max = 6.8 MV/m 
defocusing at r = 5 mm measured 1 mm in front of the cathode face) to 14.23 mm 
retracted behind the nose cone (with field E z max = 5.8 MV/rn and E rt max = 483 
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Table 3. Cathode assembly and beam pipe power dissipation due to coupling with the 
NPS SRF cavity as calculated by Superfish. Various combinations of cathode stalk 
design, materials (304 stainless steel and copper), and active cooling are considered. 
The simple model uses a right cylindrical cathode stalk. The complex model uses the 
changing radius cathode stalk. For each model, power dissipated in the vacuum tube 
wall (pipe) and along the cathode stalk (K) are shown. 


Component Arrangement 

Simple Model 

Complex Model 


Pipe [W] 

K [W] 

Pipe [W] 

K [W] 

SS Pipe (293 K), Cu K (293 K) 

3061 

2043 

121.5 

334.1 

SS Pipe (77 K), Cu K (293 K) 

1190 

2043 

47.21 

334.1 

SS Pipe (77 K), Cu K (77 K) 

1190 

794.1 

47.21 

129.9 

Cu Pipe (293 K), Cu K (293 K) 

469.0 

2043 

18.61 

334.1 

Cu Pipe (77 K), Cu K (293 K) 

182.3 

2043 

7.234 

334.1 

Cu Pipe (77 K), Cu K (77 K) 

182.3 

794.1 

7.234 

129.9 


kV/m (defocusing) at r — 5 mm and -58 kV/m (focusing) at r = 2.5 mm measured 
1 mm in front of the cathode face). Figure 26 shows these two cathode position 
extremums and the associated field lines as calculated by Superfish. It can easily 
be seen between Figure 26 and Figures 17 and 18; there is a distinct trade between 
the main accelerating gradient and the focusing fields experienced by the beam as it 
is emitted from the cathode surface. Figure 26 also plainly shows the cathode face 
is a blunt, flat-faced surface. Additional focusing during electron emission could be 
obtained by dishing this surface. 
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Figure 26. By compressing or extending the cathode assembly bellows, the cathode 
position relative to the nose cone can be adjusted. This changes the radial and 
longitudinal fields experienced across the cathode face. The upper diagrams are 
from Niowave design drawings. The bottom diagrams are held plots of the nose 
cone/cathode region of the cavity taken from Superfish. 
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4. Superconducting Solenoid 

A superconducting solenoid is included in the cavity design as a focusing ele¬ 
ment, mounted as close as practicable to the SRF cavity. Figure 20 shows the relative 
position of the solenoid to the cavity. As the electron beam exits the cavity, space 
charge effects will be causing the beam to expand radially. By placing the solenoid 
as close as practicable to the SRF cavity, we maximize the bunch charge that can be 
extracted through the cavity coupler. 

The solenoid is chosen to be superconducting as it is located near the helium 
dewar and inside the liquid nitrogen shield. A heat source at this location would be 
a significant heat load on the cryogenic cooling systems. Using a niobium-titanium 
winding (niobium wires in a titanium matrix) and conductive cooling from the helium 
dewar, when the wire temperature is below the critical temperature for niobium, a DC 
current can be passed through the solenoid with negligible ohmic heating. Conductive 
cooling is provided to the solenoid casing through two copper bus bars, one of which 
is shown in Figure 27. 



Figure 27. Copper bus bar providing conductive cooling from the liquid helium dewar 
to the superconducting solenoid. There are two identical bars with the second one on 
the opposite side of the beam pipe from the one shown. 
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To prevent magnetic fields from becoming trapped in the accelerating cavity, 
the solenoid’s iron flux clamp will be demagnetized after use and the solenoid will 
only be energized after the cavity has reached superconducting temperatures. The 
solenoid is designed for a peak field of 238 mT with a current of 5.0 A. The solenoid, 
after manufacture, was tested to handle a maximum current of twice the design value 
without quench. 

5. Assembly and Cavity Cleaning 

Once the cavity pieces have been manufactured and machining completed, 
they must be assembled in a fashion that minimizes possible defects affecting cavity 
performance. Once the cavity is assembled, the difficulty in correcting any defects that 
arise increases drastically. Defects that can affect cavity performance include weld 
spatter, dust, metal shavings, and hydrogen contamination, among many others. 

a. Cavity Assembly 

The assembly process starts with welding the cavity into its final state. 
All preliminary testing to verify the machining steps 4 must be completed prior to this 
step, as disassembly of the cavity is extremely difficult and costly. Normal welding 
processes can be used, but experience with processes such as tungsten inert gas (TIG) 
welding has been poor [12]. Electron beam welding offers a much more controllable 
environment for ensuring weld quality and minimizing incorporation of gases in the 
weld bead. Figure 28 shows the NPS cavity in an electron beam welder at Sciaky, 
Inc., in Chicago, IL. 

The weld procedure must be well planned to ensure all welds are full 
penetration and have a smooth surface at the cavity inner surface. The process must 
be controlled to prevent detuning of the cavity due to weld shrinkage (excessive heat) 


4 See Chapters IV.B.l. and IV.B.2. for results of these steps. 
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Figure 28. The NPS cavity in the electron beam welder at Sciaky, Inc. Top left is the 
nose cone being welded to the base plate. Bottom left is the base plate being welded 
to the outer conductor of the cavity. The photo on the right is a conflat flange being 
welded to the anode side of the cavity. 


or blow out (focused electron beam). Electron beam welding is conducted under 
vacuum usually less than 10 //Torr [12], which minimizes the gases present to absorb 
into the cavity materials and welds. 

b. Cavity Surface Etching 

The manufacturing process introduces contaminants in the form of ox¬ 
idation layers, foreign metal contamination, and lubricants that must be removed 
for optimal cavity performance. Removal of lubricants is best accomplished through 
use of the usual solvents used in vacuum systems. Accomplishing the removal of the 
oxidation layer and foreign metals require more aggressive methods. The standard 
method adopted for this process is the buffered chemical polish (BCP). BCP is an 
acid mixture consisting of a solution of hydrofluoric acid (HF), nitric acid (HNO3), 
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and phosphoric acid (H 3 PO 4 ) in a ratio of 1:1:2 respectively. Figure 29 shows the 
setup used to perform the BCP on the NPS cavity. 



Figure 29. This station was used for performing the buffered chemical polish of the 
NPS cavity. I11 the inset at lower right is the cavity placed in the rocking cradle that 
mates with the stand shown just to the right of the cooling tub. The hole shown in 
the inset is the connection to the liquid helium fill tube and leads to the exterior of 
the niobium cavity. 


All acids were premixed and maintained under the vent hood during 
the processing. The acid was connected to the cavity via plastic tubing attached 
to the beamline connection points (white connectors in Figure 29 inset) including 
a drain line. The cavity was filled in the vertical position with the nose cone end 
down and drain line open at the top to ensure proper filling of the entire cavity. 
Once filled, the drain and fill lines were clamped and the cavity placed in the rocker 
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frame which was in the cooling tub. The cooling tub was filled with a mixture 
of water and ice to provide cooling of the cavity during the etching process. A 
thermometer was used to measure the temperature of the acid solution inside the 
cavity and maintain it below 18 °C to control the etch rate and minimize the chance 
of hydrogen contamination [12]. During the etch, the cavity was agitated by hand to 
prevent settling of niobium removed from the cavity surfaces. The etch was halted 
after the allotted time sufficient to etch approximately 150 /jrn of surface material by 
dumping the acid followed immediately by flushing with ultrapure water (resistivity 
of 17.5 Mfi-cm as measured at the discharge of the storage tank). After flushing, the 
cavity was filled with ultrapure water and capped for transport to the clean room to 
prevent any airborne particulates from having access to the cavity. 

c. High-Pressure Rinse 

Once inside the class-100 clean room, the cavity was drained and mount¬ 
ed to the high-pressure rinse station shown in Figure 30. As has been seen in previous 
cavities, high-pressure water rinse immediately following cavity etching or electropol¬ 
ishing has shown significant reduction in field emission allowing for attainment of 
higher field gradients [28]. The water rinse removes any last traces of the acids used 
to chemically process the cavity as well as flushing any remaining particulate matter 
from the cavity walls. This must be performed in a “clean” environment to prevent 
any particulates in the air from adhering to the liquid surface and being deposited as 
the liquid evaporates, negating the benefits of the high-pressure water rinse. 

The wand consists of a single, small pipe with a series of small holes 
drilled near the tip, perpendicular to the surface. The tip of the wand is angled 
toward the exterior of the cavity and has a nozzle drilled in it to provide an angled 
spray as well. High-pressure water is provided at the opposite end of the wand 
connected to a high-pressure pump outside the clean room and taking suction from 
the ultrapure water tank. During rinsing, the wand is cycled up and down while 
the cavity is simultaneously rotated. The cavity is mounted with the nose cone at 


56 




Figure 30. The Niowave high-pressure rinse station consists of a high-pressure, ul- 
trapure water feed to a rinse wand. One the left is the stand in the vertical position 
without the cavity showing the wand in the fully extended position. On the right is 
the cavity mounted to the station with the wand in the angled position to allow for 
cleaning of the exterior of the nose cone surface. 

the top (pointed toward the mounting ring). This allows for canting of the wand 
to reach into the rear corners of the cavity and clean the exterior of the nose cone 
and coaxial line. These positions are of utmost importance as they see the highest 
electric and magnetic fields. The cavity was rinsed in each position continuously for 
approximately 1.5 hours and then allowed to air dry in the clean room. Once the 
majority of the water had evaporated, the cavity was placed under vacuum to verify 
that the cavity was still vacuum tight and to prevent any particulate intrusion in the 
now clean and prepared cavity. The cavity at this stage is shown in Figure 31 with the 
superconducting solenoid attached to the beam tube and wired for final installation. 
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Figure 31. The NPS SRF cavity in the class-100 clean room after chemical processing 
and high-pressure rinse. The cavity is under vacuum to remove the remaining water 
and to verify the cavity is still vacuum tight. Note that the beam tubes have also 
been attached to the cavity in preparation for the next stage of assembly. 






d. Cryomodule Assembly 

The cavity and beamline pipes were then brought back to atmospheric 
pressure under dry nitrogen purge and the beam tubes capped with blank flanges so 
that the cavity could be encapsulated in the cryomodule. The cryomodule consists 
of multiple layers of insulation and cryogenic shielding in order to minimize heat 
transport from the outside environment to the liquid helium. As part of the assembly 
process, multiple temperature sensors were installed to monitor temperatures in the 
various parts of the cryomodule. The most important sensors are those located on 
the solenoid casing and at the top of the dewar. The first indicates whether the 
solenoid has reached superconducting temperatures and the second can be used as 
a liquid helium level sensor. Each sensor is wired to a remote monitoring device for 
presentation to the operators. The wiring and temperature monitors are shown in 
Figure 32. 



Temperature Diode Locations 1 thru 9 


1. Thermal Intercept 

2. Cathode End Cavity Flange 

3. Solennold - Top 

4. Solenoid - Side 

5. Beam Tube 

6. Beam Tube 

7. Thermal Intercept 

8. Beam Tube 

9. He Transfer Tube 


Figure 32. Nine temperature monitoring sensors are placed inside the cryomodule to 
monitor temperatures of the cavity and support structures. These sensors read out 
to a remote monitoring station for use by operations personnel. 
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Once prepared, the cavity was moved into place under a support struc¬ 
ture for the “conning tower.” The conning tower supports the cavity inside the cry- 
omodule via the helium fill tube, allowing the cavity to hang from the conning tower. 
The conning tower is welded to the main portion of the cryomodule to complete the 
vacuum tank, providing sufficient support for the cavity. The cavity is hung from a 
truss for alignment relative to the support foundation for assembly of the remaining 
portions of the cryomodule as seen in Figure 33. 



Figure 33. The cavity is in place and aligned to the support foundation. 


After a test fitting of the various portions of the cryostat, the cavity 
was enveloped step by step. The first thermal layer after the helium dewar is three 
thicknesses of insulating blanket. The insulating blanket material is shown in Figure 
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34 and consists of a layer of expanded fabric backer attached to a metalized mylar 
surface. The backer’s expanded structure limits conductive thermal transport while 
the metalized surfaces act to reflect radiated heat from the outer layers of the cryostat. 



Figure 34. Insulating blanket material (shown on left) is wrapped around the helium 
vessel. 


Surrounding the insulating blanket, a liquid nitrogen shield is installed. 
The liquid nitrogen shield, shown in Figure 35, is made from copper with copper 
tubing brazed to the exterior for conductive cooling of the shield to liquid nitrogen 
temperatures. During operation, liquid nitrogen is constantly flowed through the 
shield to maintain it at 77 K. To ensure the liquid nitrogen system is kept liquid 
throughout, a constant drip is maintained from the outlet where it is observable by 
the operators. Reducing the temperature decreases the radiated power incident on 
the insulating blankets by a factor of >230 over room temperature. 

To minimize stray magnetic fields affecting the cavity’s superconducting 
properties, the entire helium vessel is encased in a MuMetal shield: see Figure 36. 
MuMetal is a nickel-iron alloy (~ 80% - 20%) with small amounts of molybdenum 
that has an extremely high magnetic permeability [29]. The MuMetal acts to shield 
the SRF cavity from the earth’s magnetic fields and any stray magnetism in the steel 
vacuum vessel or local stray fields from the SRF cavity. This is extremely important, 
as magnetic held lines can become “pinned” within the niobium as the SRF cavity 
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Figure 35. The copper liquid nitrogen shield in place around the cavity on the left. On 
the right is shown the attachment of positioning strips to keep the cavity suspended 
within the shield and minimizing conductive thermal contact. The black covering 
seen on the left was in place only for initial fit testing and was removed for final 
installation (removed on right). 


goes through transition. These trapped flux lines then contribute to the cavity’s 
residual resistance, increasing the ohmic losses and limiting cavity performance. 

Surrounding the MuMetal and liquid nitrogen shield is a nonmagnetic 
vacuum vessel. This vessel, shown in Figure 37, provides structural support for the 
cavity and beam tubes as well as providing significant thermal isolation from the lab 
environment. The vacuum vessel is welded to the collars at each end of the beam tubes 
connecting to the cavity. In case rework or reprocessing of the cavity is necessary, 
this encapsulation can be removed by removing a welded band and cutting the skip 
welds holding the two halves together. While simple sounding, the effort required 
to remove the cavity from its shields and rebuild the cryomodule is significant and 
should be avoided if at all possible. 

Completing the vacuum vessel is the “conning tower” shown prior to 
assembly in Figure 38. The conning tower is the access point for the cryogen fill and 
exhaust as well as the temperature sensors. Additionally, the conning tower houses 
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Figure 36. MuMetal shield installed around the liquid nitrogen shield. Note that 
superconducting solenoid is encapsulated by MuMetal. The black covering indicates 
this picture was taken during test fitting. 

ports for the guard vacuum burst disk, vacuum pumping, and vacuum monitoring. 
The insulating vacuum for this cryomodule is separate from the beamline vacuum. 
Vacuum levels in the vessel are maintained at 10 -4 to 10~ 5 Torr. 

To support the cavity when warm and during eventual shipment from 
Niowave to NPS, two small stands are installed. These stands allow the helium dewar 
to rest upon them when warm, preventing the cavity from swinging on the helium fill 
tube and exit beam pipe. The feet are small hollow tubes to minimize heat conduction 
through the cryomodule’s thermal shields (see Figure 39). When the cavity is cooled 
to liquid helium temperatures, the contraction of the fill tube is sufficient to lift the 
helium dewar from the feet which, in addition to removing the thermal path to the lab 
environment, also provides some small isolation of the cavity from external sources 
of mechanical vibration. 
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Figure 37. The vacuum vessel is arranged around the cavity and shields. 



Figure 38. The conning tower assembly prior to assembly and mating to the vacuum 
vessel. 
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Figure 39. When the cavity is warm, it rests upon two “shipping feet” that support 
the cavity in addition to the liquid helium fill tube and exit beam tube. When 
cold, the contraction of these tubes lifts the cavity from these supports, isolating it 
thermally and mechanically. 
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e. Coupler Installation 

Once the cryomodule construction was complete, the entire cryomodule 
and foundation returned to the clean room for installation of the RF coupler. The RF 
coupler was cleaned and assembled separately. The coupler was kept under vacuum 
in a surrogate beam tube for cleanliness purposes until the cryomodule was ready 
to receive the coupler. Figure 40 shows the coupler and pumping station ready to 
be mated to the cryomodule. The pumping station includes a leak-off valve for 
introducing gases into the vacuum, a residual gas analyzer, and an “up to air” valve 
for connection to a dry nitrogen purge when vacuum must be released. The vacuum 
pump is a combined roughing and turbo pump capable of maintaining vacuum in 
small volumes of 10~ 9 torr. What is seen in Figure 40, in combination with the 
cryomodule, constitutes the entire system for the initial RF testing. 



Figure 40. The coupler was maintained in a clean environment, under vacuum, while 
the cryomodule was assembled. The coupler is shown on the left. The vacuum 
pumping stand is on the right. 


After re-entry into the clean room, the cryomodule and foundation 
were thoroughly cleaned and prepared for mating with the coupler assembly. Using 
standard clean room procedures to minimize particulate generation, the surrogate 
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beam tube was removed from the coupler and the output beamline of the cryomodule 
was opened. Figure 41 shows the coupler after the surrogate beam tube is removed. 
The coupler is polished, annealed copper cleaned with ultrapure water prior to being 
placed under vacuum. 



Figure 41. The RF coupler assembly with surrogate beam tube removed. The tube 
is supported at the visible flange by a teflon spider. 


The coupler stand is set upon teflon rails of the same thickness as the 
aluminum base plates that will mount the coupler stand to the foundation. Teflon rails 
are used to facilitate manual movement of the coupler into the cyromodulc beam tube 
while minimizing particulate generation. Any particulate contamination of the cavity 
will increase the risk of degraded performance due to held emission sites. Particulate 
contamination of the coupler could result in multipacting, limiting the obtainable 
cavity fields. Figure 42 shows the insertion process. Technicians must ensure the 
coupler remains centered in the beam pipe and parallel to the axis during insertion to 
prevent impacting the beam tube and possibly damaging the coupler surface. Once 
in place, the flange between the sections is tightened, and the coupler assembly and 
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cavity are placed under vacuum. Final alignment of the coupler requires RF power 
and will take place outside of the clean room. Until the cathode assembly is ready 
for installation, the upstream end of the cryomodulc will remain closed with a blank 
flange. 



Figure 42. Insertion of the RF coupler into the cryomodule beam tube. 


After leaving the clean room, the cryomodule received a coat of paint 
and remaining ancillary systems were connected to prepare for initial RF testing. 5 
Figure 43 shows the Niowave facility set up for initial RF testing of the cavity. Figure 
44 shows a close up of the cryomodule and beamline with the ancillary equipment 
connected. On the left is the pumping station as described previously. The metal 
shroud around the shorted end of the coupler is supporting structure for compressing 
and extending the bellows section for adjustment of the coupling strength to the cav¬ 
ity. Under the RF coupler is the input RF power feed from a 100 W amplifer. The 
amplifier is driven by a 500 MHz signal generator controlled at the operator station. 


5 Initial RF testing results are discussed in Chapter IV.B.4 









Directly under the cryomodule is an ionization chamber for measuring radiation gen¬ 
erated during RF testing. The black tubes coming from the top of the cryomodule 
are insulated liquid nitrogen lines connecting to a liquid nitrogen tank seen in Figure 
43. 
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Figure 44. The completed 500 MHz SRF gun/booster ready for RF testing. 
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/. Cathode Assembly Installation 

After initial cyrogenic and RF testing, the cathode assembly was in¬ 
stalled in preparation for further RF processing and first beam test. Prior to installa¬ 
tion, the cathode assembly was cleaned, and the cathode stalk placed in a surrogate 
beam pipe while in the clean room to maintain cathode cleanliness. The cathode 
end of the cryomodule was encapsulated in a portable clean room to minimize in¬ 
troduction of particulates into the cavity during the cathode assembly installation 
process. Figure 45 shows the cathode assembly inside the portable clean room prior 
to installation. 



Figure 45. The cathode assembly prepared for installation into the NPS cavity. The 
assembly is under a portable clean room. The stainless steel surrogate beam tube 
and angle support are removed prior to installation. 


The cathode assembly support structure is mounted upon the same 
foundation as the cryomodule. Aluminum pads will be placed between the transverse 
foot of the assembly support structure after it is in its final position, post alignment. 
To assist in minimizing friction during alignment and sliding the cathode assembly 
into the cavity, the support structure is placed on two teflon rails (similar to the ones 
used during the coupler insertion) with the same thickness as the aluminum pads, as 
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seen in Figure 46. The structure is top-heavy and extreme care must be taken during 
insertion to prevent the cathode assembly from striking the beam tube or interior of 
the cavity. 



Figure 46. Two thin teflon rails are used to minimize friction when sliding the cathode 
assembly structure into place. The step behind the foundation is to allow the structure 
to slide behind the cryomodule with room for the cathode to clear the beam tube prior 
to insertion. One of the aluminum pads for mounting the structure to the foundation 
can be seen near the top of the photo, between the rails. 


To prepare for removal of the surrogate beam tube, the bolts in the 
flange are removed using standard cleanroom procedures and cleaned. To maintain 
a seal between the flanges, the top and bottom bolts are reinstalled after cleaning, 
but only torqued to finger tight. Figure 47 shows the flanges being prepared for 
separation. Even though the flanges have been kept in a clean environment since 
removal from the clean room, the process of removing the bolts from the flanges 
creates particulates that could enter into the cavity once the gasket seals are broken. 
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Prior to opening the cavity and cathode to the environment of the portable clean 
room, all mating surfaces are carefully cleaned. 




Figure 47. The surrogate beam tube on the cathode assembly and the blank flange 
on the cyromodulc are prepared for removal. All surfaces are wiped down prior to 
opening the beam tubes to minimize intrusion of particulates generated by removing 
bolts from the mating flanges. 


The first major issue in the assembly of the NPS gun and beamline 
was encountered as the surrogate beam tube was being removed from the cathode 
assembly. Recalling Figure 25, there is a teflon spider located approximately halfway 
down the cathode assembly to provide support. The spider provides a nonconducting 
supporting surface for the cathode assembly, decreasing the length that is cantilevered. 
In this case, as the surrogate beam tube was being slid along the cathode assembly 
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axis, the teflon spider impinged on the weld bead between the beam tube wall and the 
connecting flange. Significant force was necessary to remove the surrogate beam tube, 
resulting in deformation of the spider. Figure 48 shows a series of photos documenting 
this issue. 



Figure 48. Top left is the initial onset of the binding of the spider. Top right is a 
close-up of the spider as technicians attempt to slide the surrogate beam tube off the 
cathode assembly. Bottom left is a picture of technicians cleaning the spider as it 
will be inserted into the cryomodule along with the cathode assembly. Bottom right 
shows the deformation of the spider after it was removed from the surrogate beam 
tube. 


Once the teflon spider was successfully extricated from the surrogate 
beam tube, it was inspected to ensure it was still safe to insert into the cryomodule. 
As no damage other than deformation of the spider’s arms was noted, it was deemed 
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safe to proceed. The spider was gently cleaned and the cathode assembly was moved 
into position for insertion into the cavity. 

Figure 49 shows the cavity beam tube being opened for insertion of the 
cathode assembly. This figure also shows the cathode button that is installed in this 
assembly. The cathode button is a press-fit piece of small grain niobium whose surface 
was prepared using a buffered chemical polish. No further preparation of the cathode 
surface was performed to this point other than maintaining it in a clean environment 
prior to insertion into the cavity. 



Figure 49. The last seal before opening the cavity beam tube is broken to insert 
the cathode assembly. Note the niobium cathode button is visible at the tip of the 
cathode assembly. 


Technicians carefully maneuvered the cathode assembly along the teflon 
rails, ensuring the cathode stalk entered the cavity beam tube coaxially until the 
cathode tip was visible via a transit from the end of the beamline (described shortly), 
as seen in Figure 50. As the teflon spider began to enter the beam tube, the process 
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was again disrupted. We also see in Figure 50, the spider, due to its deformation, 
was now an interference fit in the beam tube, preventing the assembly from smoothly 
entering the cavity. Upon initial insertion, the spider, rather than entering the beam 
tube, merely slid along the cathode stalk. 



Figure 50. On the left, technicians begin to move the cathode assembly into the cavity. 
On the right, the teflon spider, due to deformation, now no longer has clearance inside 
the beam tube to easily slide without contact. In this photo, we can see that after 
pulling the cathode assembly back out slightly, the spider has slid along the assembly 
from its intended position (where the holes are visible). 

To ensure the spider was carried into the tube as near as possible to 
its intended position, a flexible PVC welding rod was used as a pusher. Using this 
method, as shown in Figure 51, one technician continually attempted to keep the 
teflon spider in its detente while the other slowly inserted the cathode assembly 
further into the cavity. It is estimated that the spider is within 2-3 cm of its intended 
position based on the flange gap when the pusher had to be removed. While less than 
ideal, this location should provide adequate support to the assembly. It is unknown 
but unlikely that the spider will shift position as intended when the cathode position 
is adjusted using the bellows. 

Once the cathode assembly was inserted far enough for the flanges to 
mate, seen in Figure 52, the system was made vacuum tight and the final setup of 
the adjustment mechanism for the cathode position completed. Unfortunately, this 
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Figure 51. A PVC pusher is used to position the spider during insertion. 


version of the adjustment system used custom bellow sections and both flanges at the 
cavity interface are nonrotatable. To obtain satisfactory bolt line up, a minimal torque 
has been applied to the adjustable bellows sections of the cathode assembly. This 
should not present a limitation in the system operation, but it should be noted that a 
rotatable connection to the cavity for follow-on builds would greatly ease connection 
of the cathode assembly to the cavity. Once under vacuum, the cavity and cathode 
assembly achieved vacuum in the microtorr level within an hour, indicating that all 
connections were leak tight and no out-gassing contaminants had been introduced. 

Alignment of the cathode to the nose cone is performed using an optical 
transit, seen in Figure 53, through an optical window in the end of the beamline. 
During installation of the cathode assembly, direction from this station was provided 
to the technicians to ensure the cathode button was properly centered in the nose 
cone. 
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Figure 52. The cathode assembly is mated to the cavity and a vacuum tight connec¬ 
tion is made. The cavity was then placed back under vacuum. 

By providing supplemental lighting to the interior of the beam line, 
features that can be used for determining where the transit is focused become visible. 
This supplemental lighting is provided by two small flashlights, one illuminating the 
beam pipe just downstream of the laser mirror and the other shining through the 
laser port and bouncing off the laser mirror to the cavity. By adjusting the focus of 
the transit, one can orient oneself to the interior of the beamline as shown in Figure 

54. The cathode button was visible (it had a very high contrast ratio) almost from 
the moment it was inserted into the beamline. 

Once the cathode was in place, the drive laser needed to be aligned to 
the cathode. In previous work, the laser technicians had performed alignments up to 
the laser window on the beamline. To assist in verifying alignment, NPS provided a 
UV sensitive camera and UV telephoto lens to image the cathode, as seen in Figure 

55. First, a helium-neon (HeNe) laser operating in the visible spectrum was used 
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Figure 53. Dr. John Lewcllen observes the positioning of the cathode assembly relative 
to the cavity nose cone using an optical transit. 

to align the final steering mirrors to put the HeNe spot onto the cathode surface. 
Once the HeNe was seen to be properly centered on the cathode, a flip-in mirror was 
adjusted to take the HeNe out of the final steering enclosure to a distant point and 
the UV drive laser aligned to shine coincident with the HeNe spot. Once the flip-in 
mirror was removed, both spots should shine on the same spot on the cathode with 
a high degree of confidence. With the HeNe laser interrupted and the focal length of 
the telephoto lens adjusted (it is chromatic), we were able to verify the UV laser was 
striking the cathode in the correct position. 

After completion of the assembly and alignment activities, the NPS 
gun/booster is ready for RF testing and conditioning leading to an initial beam test. 
The various modeling and testing of the cavity will be discussed in the following 
chapter. 
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Figure 54. Changing the focus of the transit allows one to see all the interior com¬ 
ponents of the beamline or to focus on the cathode face. In the upper left is the 
in-vacuum laser mirror. In the upper right is the shorted end of the RF coupler. In 
the picture in the lower right, the grey area is the niobium cathode button, the bright 
ring around the button is the edge of the copper recess in which it fits, and the faint 
bright ring to the left of the button is the cavity nose cone. 
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Figure 55. After aligning the cathode in the cavity, verification that the laser spot is 
striking the expected point on the cathode is still required. These photos are taken 
of both a HeNe spotting laser and the UV drive laser on the cathode. The top photo 
is the HeNe laser. The bottom photos are the UV laser spot with an open aperture 
(left) to show the multiple reflections of the scattered light and a closed aperture 
(right) for fine locating of the primary beam. 
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IV. 


MODELING AND EXPERIMENT 


Modeling and experimentation are essential to any research project, especially 
as the size or complexity of the project increases. In the case of the NPS 500 MHz 
gun/booster, the prototype is a one-of-a-kind device that is also the first of its kind. 
There is no previous experience using the quarter-wave cavity type for the acceleration 
of electrons, nor for any other particle in an axial configuration. Significant time 
and resources were spent in developing accurate representative models and running 
numerous simulations to understand how the systems should react before attempting 
the same setup and tests on the actual hardware. 

This combined approach means that during the experiment, researchers have 
a much better understanding of the results they should be seeing. Experimentation 
also provides grounding for the models, allowing them to increase their representative 
value to the experimental researcher. The NPS research team has followed this ap¬ 
proach throughout the development and testing of the NPS gun/booster. Results of 
both simulation and experiment are presented in the chapter that follows, and where 
applicable, comparisons are made between the two. 

A. MODELING 

Modeling provides an environment where the researcher is able to bring almost 
any variable into the design space. It can allow researchers to examine parameter 
spaces beyond what their experimental environment allows. On the other hand, 
it is possible to develop a model that violates or fails to incorporate fundamental 
rules such that simulations no longer have a basis for comparison to experiment. 
The researcher must be extremely careful to ensure his models are grounded to the 
experiment of interest and that simulations do not exceed the capabilities of the 
programming developed to run them. 
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1. Simulation Codes Used 

Numerous simulation and analysis software packages were and continue to 
be used in the design and development of the NPS superconducting radiofrequency 
(SRF) gun/booster. During the design phase of the cavity, the majority of the simu¬ 
lation work was performed using the Poisson Superfish family of codes developed by 
the Los Alamos Accelerator Code Group. 1 Particle interactions with the cavity have 
been modeled in Parmela and General Particle Tracer (GPT) and benchmarked occa¬ 
sionally against other accelerator codes. We have also occasionally used new software 
packages to evaluate their usefulness by comparing their results to data already ob¬ 
tained in previous simulation work on this cavity. We will primarily focus on Poisson 
Superfish, Parmela, and GPT as these software packages were used for the bulk of 
the work, but will mention other codes used for benchmarking where appropriate. 

Sub-applications within each code for massaging of input decks or external 
data manipulation for analysis will not be discussed as they are not integral to under¬ 
standing the simulations performed. Additional applications, such as the SDDS (Self- 
Describing Data Set) toolkit from Argonne National Labs, 2 MATLAB, and Microsoft 
Excel were used to sort, manipulate, and present data at various times throughout 
this research. Optimization and parameter searches were either performed by simple 
self-generated scripts or through the use of Tel 3 (a cross-platform scripting language) 
optimizer and sequencer scripts provided by Dr. John Lewellen of the Naval Post¬ 
graduate School. 

All simulations were performed on Intel-based Apple Macintosh computers. 
Large simulations were handled using a Mac Pro with two 2.8 GHz Quad-Core In¬ 
tel Xeon processors and 32 GB of RAM. Where necessary to use a Windows-based 

1 Superfish documentation and software can be downloaded from the LAACG at 
http://laacgl.lanl.gov/laacg/services/download sf.plitml. 

2 The latest versions may be downloaded from Argonne National Lab at 
http://www.aps.anl.gov/Accelerator_SystemsJDivision/OperationsWnalysis/software.shtml. 

3 ActiveState at http://www.activestate.com/activetcl provided the free Tel installation used in 
this work. 
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environment, Parallels version 5 was used to host a Windows XP Pro 64 bit environ¬ 
ment with access to 6 processors and 8 GB of RAM. Additional small simulations 
and parameter searches were performed on a MacBook Pro with a 2 GHz Intel Core 
Duo processor and 2 GB or RAM. Parallels version 5 hosting a Windows XP Pro 
32 bit environment with one processor and 1 GB of RAM was used where necessary. 
Simulation comparisons between the two computers were excellent except in the time 
required to complete the computation. 

a. Poisson Superfish 

Poisson Superfish is a family of codes generally used for modeling of 
electromagnetic systems [30]. Poisson is typically used for electro- and magneto¬ 
static models, while Superfish is designed to handle RF electromagnetic fields. Both 
of these solver programs use either two-dimensional Cartesian or axially symmetric 
cylindrical coordinates. Poisson was used to model the magnetic fields of the super¬ 
conducting solenoid and Superfish was used to model all RF cavities examined in this 
work. Other elements, including a normal conducting solenoid, window-frame beam 
steering magnets, and quadrupoles were modeled in Poisson as well. Since Poisson 
and Superfish solve different problems, we discuss their approaches separately. 

Poisson, the static held solver, takes a text input hie describing the 
geometry of a problem and the relevant properties such as boundary conditions and 
materials and solves for the vector potential, A. The hrst step of the algorithm is to 
mesh the area of interest using an irregular triangular mesh. Each triangle touching a 
boundary will have two vertices touching the boundary, altering the boundary shape 
to line segments approximating the intended curve. Poisson then solves Ampere’s 
Law at each mesh point, 

(j H (r) ■ dl = l ( V X A — M) d 1= f J(r)-ndA, (IV. 1) 

JC JC \ ho J Js 

where C is the contour of integration, H is the auxiliary magnetic held, M is the 
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magnetization, J is the current density, S is a surface of area A enclosed by C, and 
h is the outward normal. It approximates that over each triangular area the vector 
potential, A , is linear, the reluctivity is constant, and the current density is constant. 
From H, the program can solve for B , the magnetic held by 

B(r) = VxA(r). (IV.2) 

Additionally, the program uses the Coulomb gauge, 

V • A (r) = 0 . (IV.3) 

The specifics of the algorithm for solving equation (IV. 1) can be found in [31], but 
in general, for each mesh point the program solves a contour integral through the 
surrounding six triangles. The value for the mesh point is updated and the process 
continues through the entire mesh iteratively until a convergence criterion for the 
vector potential is satisfied. Boundary conditions are applied as either Dirichlet (the 
value of the vector potential is given on the boundary) or Neumann (the normal 
derivative value of the vector potential on the boundary is given) boundary conditions. 
These conditions are specified in the input hie. 

Superfish also solves two dimensional cartesian (assuming z dimension 
is infinite) and three dimensional cylindrically symmetric problems. It, however, 
focuses on solving RF problems. Primarily, Superfish solves for resonant frequencies 
of cavities and waveguides and then determines the electromagnetic fields resulting 
from the problem geometry and frequencies of interest. 

Superfish starts from the generalized Helmholtz equation 

V 2 0 + k 2 (j) = S , (IV.4) 

where (f>(x) are the eigenfunctions and S' is a source term. The modes (TE or TM) 
of a cavity are determined from the boundary conditions specified by the user. The 
geometry of the problem is meshed exactly the same as for Poisson. The constitutive 
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relations for the assumed homogeneous, isotropic, non-conducting medium in the 
cavity are 


D = eE (IV.5) 

B = fiH , (IV.6) 


where e is the permittivity and /i is the permeability as input by the user in the 
problem definition. 

For solutions presented in this work, a field drive point is provided by 
the user, which allows Superfish to use the auxiliary magnetic held as the eigenfunc¬ 
tions, 

V x (V x H) - k 2 H = 0 . (IV.7) 

Applying Stokes Theorem, we get the integral equations 


f V x (V x H) ■ h dA = / (V x H) ■ dl = k 2 f H ■ h dA , (IV.8) 

J S C s 

where S is the area bounded by the contour C and h is the outward normal of the 
surfaces. The same process of integrating around the mesh triangles surrounding 
the point of interest is implemented in the Superfish solver and iterating until the 
convergence criteria is met. 

To find resonant frequencies, Superfish assumes a fictitious magnetic 
current density K, which is used to drive the fields in the cavity. At resonance, the 
current required to drive the fields goes to zero. The function 


D(k 2 ) = kc 


j H K dV 
fe(H-H) dV 


(IV.9) 


is used to search for zero crossings. A zero crossing is necessary but not sufficient for 
determining a resonance as 


d D (k 2 ) _ 
d (k 2 ) 

must also be true for a particular crossing to be a resonance. 


(IV.10) 
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When the user provides a frequency of interest, Superfish will solve 
for the fields at a resonance near this frequency. Subroutines can also be used to 
find many pertinent values of interest including energy per unit volume in the cavity, 
power losses on the cavity walls, average accelerating held, shunt impedance, cavity 
quality factor, maximum held values along boundaries, and many other parameters 
based on the user’s input assumptions. Of critical importance is the ability to output 
held prohles and maps that can be imported into other codes, such as Parmela and 
GPT for use in modeling particle trajectories. 

b. Parmela 

Parmela is a particle accelerator design and simulation code developed 
at Los Alamos National Lab beginning in the 1980s. As with Poisson Superhsh, the 
code is maintained and distributed by the Los Alamos Accelerator Code Group. 4 
Parmela uses an input hie that describes the accelerator geometry, helds, particles, 
and related subroutines to be invoked. The internal library of accelerator equipment 
including cavities and correctors is extensive, but arbitrary helds may also be im¬ 
ported. Particle distributions of multiple species and charge signs can be invoked 
simultaneously and distributions may be generated within Parmela or imported from 
an external source. The code tracks particles through the accelerator based upon a 
phase angle clock based on a designated RF frequency rather than an explicit time 
reference. Particle transport is handled by various subroutines based on the dynamics 
invoked in the input hie. The number of particles in a simulation is limited only by 
the available memory of the computer performing the simulation. Each specihed par¬ 
ticle is really a macroparticle representing many charged particles, but maintaining 
the appropriate charge to mass ratio for particle dynamics. 

RF held conventions in Parmela are based upon a “sine” convention in 
which 0 = 0 indicates a zero crossing. Since Parmela assumes a particle charge of 

4 Parmela may be downloaded, after approval, from the LAACG at 
http: //laacgl.lanl.gov/laacg/services/services.phtml 



+1 unless designated otherwise, the building RF fields from this point represent an 
accelerating held in the positive z direction. There is no issue assuming a +1 charge 
for electrons, so long as the convention is maintained throughout (i.e. reversing holds 
for bending magnets). 

Particle tracking is handled through an internal database that compares 
a particle’s position with the beamline geometry. When particles are lost due to 
interaction with the beam line geometry or beam line elements, the particle’s database 
entry is annotated as “lost” and is no longer tracked. 

Particle motion is imparted by the code based on local and global held 
elements and interactions with other particles. At each phase step (or specihed num¬ 
ber of steps), each particle receives an impulse resulting from space charge, helds from 
the local element, and distributed helds. Space charge helds are probably the most 
difficult to understand, and Parmela has many implementations to consider. The 
two primary forms are 2D and 3D. The 2D form is used primarily for axisymmetric 
beams with particles distributed upon a 2D mesh and Poisson’s equation solved for 
space charge forces. In the 3D form, a mesh is formed around the bunch and used 
for calculating particle-particle interactions. In all cases, Parmela transforms to the 
moving frame of the bunch to perform the particle-particle space charge interactions. 

Much more detailed information can be found in the Parmela users 
manual [32], There is signihcantly more detail required to learn to use the program 
efficiently, it is only our intent to provide a general understanding of the program 
sufficient to understand the results presented. Data output is to a binary data hie 
for use by Parmela’s dedicated interpreter program. For this work, a translator was 
used to transfer this binary data to SDDS format for use with previously developed 
analysis and manipulation tools. Parmela is a Windows command line program only. 

c. General Particle Tracer 

General Particle Tracer (GPT) is similar to Parmela in that it tracks 
particles through input helds and geometries, but has a more generalized approach as 



its use is not intended to be specifically limited to particle accelerator studies. GPT 
also uses the macroparticle approach identical to Parmela. Many of the features 
and interactions with the executable code are similar in that a geometry and holds 
are specified from the internal libraries or imported from external sources. Multiple 
species and charges may be handled simultaneously. Particle distributions may be 
generated internally or imported. Particle transport is solved in the time domain by 
a 5th order embedded Runge-Kutta integrator with adaptive step size control [33]. 

Data output is to a binary database that can be directly interpreted by 
the Windows GUI environment or translated to text, SDDS format, or other specific 
formats via included translators. Data acquisition is initiated by user-stipulated times 
or positions. In time-specified output, all particle data is collected at each specified 
time. For positional output, non-destructive “screens” are placed at the required 
position and particles recorded as they pass through the desired location. GPT will 
interpolate to determine the coordinates and time when the particle passed through 
the desired screen location. 

Numerous space charge routines are available for use, including 2D 
and 3D models similar to Parmela’s. The two routines predominantly invoked in 
this work are the spacecharge3D() and spacecharge3Dmesh() routines. The mesh 
version calculates space charge forces in the rest frame of the bunch approximating the 
charge density at the nodes of a non-equidistant mesh encompassing the bunch. The 
charge density is then fed into a Poisson solver to obtain the electric potential, which 
is differentiated to obtain the rest frame electric fields. These fields are combined with 
the system geometry fields and applied as impulses to the particles. This particular 
routine was developed and optimized with high charge and high brightness bunches 
in mind [34], The spacecharge3D() routine is a particle-particle solver using the 
fully relativistic held equations. It makes no approximations, but has a significant 
cost in CPU time over the meshed approach [33]. This routine was used only when 
modeling low energy beams or when very few particles were being modeled. 
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Of particular interest is that GPT is fully customizable. If a particular 
beam line element is used repeatedly, it is possible to write a custom element that 
can be compiled into the GPT program for use. GPT can be executed in a Windows 
graphical user interface, Windows command line, and Linux environments. 

2. Beam Simulation Methodology 

To ensure comparability between beam simulations using different software 
packages and different component models, a standard methodology was employed. 
In general, it is very difficult to perform multidimensional parameterization by hand. 
To simplify this, beam lines were approached in sections, limiting the number of free 
parameters available for analysis and adjustment and terminating the propagation 
at specific points associated with the sections to limit processing time. Parameters 
controlled by sections previous to the section of interest are held fixed and an ini¬ 
tial parameter search through the current chunk is performed. After analysis of the 
resulting data, a carefully chosen starting point and limited parameter space is pro¬ 
vided to the optimizer code. The optimizer attempts to adjust the free parameters 
it is given to minimize a figure of merit (FOM) provided by the user. The resulting 
optimized variables are incorporated and the next section is taken into consideration. 

As an example, let us consider optimization of beam radius in the context of the 
500 MHz gun. At the start, a cathode position, bunch charge, beam shape, etc., have 
already been determined. The first section is the accelerating cavity itself (terminating 
prior to the next major beam line element—in this case, the superconducting solenoid) 
and the variable of interest is the cavity launch phase. For a FOM, we desire to 
maximize the average energy of the beam, minimize the spot size at the termination, 
and the transport the maximum charge, with transport being more important than 
energy or spot size. To account for this, we set up the FOM, 

FOM= A , (IV.ll) 

S (tf) 

where a r is the beam radius, E is the average kinetic energy of the beam, N in is the 
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number of particles in the beam at instantiation, and N out is the number of particles 
at the last data recording location. Note that the transmission of particles is raised 
to the fifth power to account for the severe penalty applied for lost particles. The 
number of particles chosen is dependent on the task being performed—in general, a 
parameter search is performed with fewer particles than an optimization run since 
many more runs are typically performed for a parameter search. 

The output data from this first section is analyzed and a range of viable 
phases is provided to the optimizer with an estimated best value and an initial step 
size. The optimizer then looks for a minima in the FOM by adjusting the allowed free 
parameters. Once an optimal launch phase is found, that value is held fixed for the 
next chunk. In this example, that would be propagation from the superconducting 
solenoid to the normal conducting solenoid. 

For the second section, the applied magnetic held of the solenoid is allowed to 
vary to optimize transport and spot size (meaning the system will attempt to focus 
at the final data recording point). The same parameter search process is applied to 
the new variable and an optimal magnetic held determined. The hnal section is to 
include the normal conducting solenoid and terminate at some down stream point 
(taken in the following simulations to be the assumed location of the next beam line 
element). 

Once parameter values have been determined through the entire beamline, 
these are used as starting values for an entire system optimization. In general, the 
values found in the above process provide good estimates of the optimized values, 
however some mod if ication or interaction by the user is necessary. For instance, 
attempting to determine the optimal focusing of a solenoid prior to a booster cavity 
requires looking at both the magnetic held strength of the solenoid and the cavity 
phase simultaneously, as there can be multiple local minima in the resulting parameter 
surfaces. This requires external bounding to ensure the optimization routine remains 
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near an optimal solution. Where deviations from this general process occurred, they 
are described in the simulation description. 


3. Space Charge 

The primary force making transport difficult in high bunch charge, low energy 
beams is the self-force of the electrons in the beam interacting with each other. This 
is normally termed “space charge.” As one attempts to push the charges closer 
together, creating smaller/shorter bunches, the forces acting to expand the beam 
increase significantly. For high power FEL applications, it is desirable to push the 
peak current of the beam as high as possible. 5 This requires a beam with high-bunch 
charge, and short temporal extent. Since this concept is so important to simulations 
performed in this research, it is worthwhile to review some of the equations and some 
simplified models to develop some instinct into how space charge affects beams of 
interest. 

We can use continuous beam theory as a zeroth order approximation to the 
moderate bunch lengths possible with the current NPS drive laser, as the bunch length 
is long compared to the gaps between the electrodes in the Stanford gun. Additionally, 
insight developed here will be useful when our focus shifts to RF guns. We start with a 
simplified derivation of the envelope equation [4] and compare its results with a simple 
simulation that includes longitudinal effects. To start the derivation, we assume a 
uniform charge density, laminar (no large divergence), steady state current with the 
geometry as shown in Figure 56. 

This current has two self fields in the lab frame—a radial electric field described 
by Gauss’ law and a magnetic held described by Ampere’s law in SI units, 


E r 

Be 


I r 

2-jie 0 v R 2 
ho I r 
2vr R 2 ' 


(IV.12) 
(IV.13) 


5 FEL theory explaining this will be detailed in the next chapter. 
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Figure 56. A current distribution of radius, R, moving with velocity, v, having current, 
/. 


Applying the relativistic Lorentz force, we find that the radius of the beam changes 
over time as 


d 2 r 
d t 2 
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Changing from time to propagation distance, z = c/3t, we rewrite equation (IV. 14) as 

5 - 2 ^W^ • 

The leading coefficient can be defined as the generalized perveance, 

K _ 

27re 0 777c 3 /? 3 7 3 ’ 

which is a measure of the space charge forces present in a relativistic charged particle 
beam. 

If we now include a driving term (such as from a solenoid field) to equation 
(IV. 14) with coefficient k 2 } proportional to the particle beam radius and an omittance 
term, we are very close to the envelope equation. To get to the transverse envelope 
equation, we only consider how a particle at the outer edge of the beam will behave, 
and we have 

if" + -A - A = 0 . (IV.16) 
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The normalized omittance, e n , can be viewed as a measure of the internal disorder 
of the beam. A higher emittance beam will diverge faster than a lower omittance 
beam. In particular, the normalized emittance is a statistical quantity of trace space 
that describes the area encompassed by the group of particles making up the beam. 
Mathematically, the normalized emittance is defined 

e n = Afi'yyJ(x 2 ) (O') 2 } - (xx') 2 , (IV.17) 

where x is the transverse coordinate displacement from the axis and x' is the trans¬ 
verse velocity of a particle. 

The envelope equation gives the researcher insight into how the maximum 
extent of a well-defined, long beam behaves under space charge forces. But with the 
simulation tools available, we can also look at some more representative beams and 
how they behave as compared to the envelope approach. To start, we define bunch 
geometries of interest that will be used in further simulations in this work. Figures 57, 
58, and 59 show pertinent data for three basic beam shapes, a uniform distribution— 
uniform distribution longitudinally and radially, a Gaussian distribution—Gaussian 
longitudinally and radially, and an asymmetric shape—a fast rise and slower fall 
longitudinally and Gaussian radially. To model the envelope equation, we use a 
simple MATLAB script, setting the emittance to zero so that radial growth is due 
solely to space charge forces. 

In Figure 60, we see three beam comparisons. The solid lines are best repre¬ 
sentative of the assumptions made when we derived the envelope equation—that the 
beam is long compared to its radius. As the beam is shortened significantly in the 
longitudinal direction as compared to the radial direction, the envelope equation sig¬ 
nificantly overestimates the radial beam expansion. The same value of charge is used 
for each calculation, however the GPT simulations take into account the increased 
longitudinal forces in the “pancake” beam case, quickly decreasing the charge density 
with a resulting decrease in radial beam expansion. Additionally, we see between 
the left and right plots the effects of the y -3 dependence in the space charge term as 
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time=1e-006 




Figure 57. A distribution that is longitudinally and radially uniform. The histograms 
show particle distributions in z (upper) and x (lower). The distribution is shown after 
generation using GPT. 

there is almost a two order of magnitude reduction in beam expansion of the envelope 
equation estimates with a two order of magnitude increase in energy (KE goes from 
10 keV —■> 1 MeV, goes from 0.942 —■> 0.0387). 

The longitudinal forces are important as photoinjectors have become fairly 
common as electron sources to linear accelerators. Even with the 500 MHz operating 
frequency of the NPS cavity, a 40 ps pulse represents 7.2 degrees of phase for the 
cavity fields. As the cavity frequency increases, the pulse length (in phase) must 
decrease to avoid energy spread in the beam. Ultrashort pulses with high bunch 
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time=1e-006 



time=4e-011 



Figure 58. A distribution that is longitudinally and radially Gaussian. The his¬ 
tograms show particle distributions in z (upper) and x (lower). The distribution is 
shown after generation using GPT. 


charges have been shown to lead to elliptical electron bunches having extremely de¬ 
sirable propagation characteristics [35], thus, being able to generate these type pulses 
is desirable. 

Knowing the longitudinal expansion of the pulse is important, and in the short 
pulse length regime, possibly desirable for propagation characteristics. We compare 
the three pulse types with that of the long pulse generated envelope equation. Using 
a common beam energy of 100 MeV, a bunch is started from the origin and allowed 
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time=1e-006 



nme=4e-d09 



Figure 59. A distribution with a fast rise and slow fall longitudinally and a Gaussian 
shape radially. The histograms show particle distributions in 2 (upper) and x (lower). 
The distribution is shown after generation using GPT. 


to propagate in free-space under only the influence of internal space charge forces. 
Table 4 shows the resulting maximal radial, cr r , and longitudinal, a z , dimensions of 
the beams as simulated in GPT. 

As expected, the long, thin beam expanded significantly less than the other 
higher charge density beams. What is surprising is how little difference there is 
between the Gaussian and asymmetric shapes. They both have central, or near 
central, high density areas that quickly expand under space charge forces. Accelerated 
radial growth is slowed by the rapid reduction of charge density as the short pulses 
“blow out.” When compared to a long pulse with an equivalent charge density (using 
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Figure 60. Space charge effects on different aspect ratio beams at low and high 
energies. Blue plots are derived from the envelope equation and red plots are results 
from GPT simulations. The plot on the left is at 1 MeV energy, and the plot on the 
right is at 10 keV energy. A uniform bunch charge distribution is used. 

the short, uniform pulse as a basis) of ~1.1 pC/mm 3 , we see a radial expansion beyond 
that of the short pulse. This is probably due to the fact that the charge in the middle 
of the beam, where the radius increases the most, cannot decrease the bunch charge 
density as quickly due to the longitudinal extent of the beam preventing the pulse’s 
longitudinal expansion. 

It is apparent from our studies of space charge that external forces are re¬ 
quired to keep the electron pulses we generate from expanding significantly. For FEL 
purposes, the peak current is of critical importance to the FEL interaction. As the 
beam expands radially and longitudinally, not only does transport become more dif¬ 
ficult (due to nonlinearities in the focusing and bending element holds), but the peak 
current decreases as well (due to having the same charge in a longer, wider electron 
pulse). External forces for transport, such as solenoids and quadrupoles, to refocus 
the beam are required, especially when the beam is at low energy and space charge 
forces are significantly more important. Additionally, we can mitigate some of the 
longitudinal expansion of the beam by accelerating to higher energy as soon as pos¬ 
sible, thus increasing 7 . If one can induce a negative energy chirp (the head of the 
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Table 4. Many photoinjectors operate in a “pancake” regime where the electron pnlse 
is short compared to its radius. Comparisons are presented between a 1 nC long pulse 
with length = 100 R and 1 nC, 40 ps pulses of various pulse shapes at a kinetic energy 
of 100 MeV. All “final” measurements are taken after propagating 100 m. 


Beam Type 

a r Start 

<j r Final 

a z Start 

a z Final 

Long Thin 

0.005 m 

0.436 m 

0.50 m 

0.62 m 

Short Uniform 

0.005 m 

2.23 m 

0.0113 m 

1.14 m 

Short Gaussian 

0.005 m 

2.99 m 

0.0085 m 

1.48 m 

Short Asymmetric 

0.005 m 

3.16 m 

0.0072 m 

1.85 m 

Long, high charge density 

0.005 m 

3.60 m 

0.50 m 

1.49 m 


beam at lower energy than the tail), a dispersive section, such as a chicane, can be 
used to decrease the pulse length, increasing the peak current. 
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4. DC Gun Analysis 

The first experimental system that was modeled during the process of devel¬ 
oping an injector section for the future NPS FEL was the Stanford DC gun. Initially, 
the plan was to use the Stanford gun as the electron source with the SRF 500 MHz 
cavity as a booster to bring the energy to near 1.5 MeV prior to injection into the 
main linear accelerator. As part of the experimentation process, models were devel¬ 
oped for the gun and associated diagnostic beam line by the author in preparation 
for experimental verification. 

a. Stanford DC Gun 

The former Stanford DC gun was initially installed as an electron source 
for a small diagnostic beam line. The hope of the experiment was to develop experi¬ 
ence running high voltage systems as well as testing the viability of the gun for use 
as the NPS FEL gun. The electron source is a gridded thermionic cathode embedded 
within a static voltage accelerating structure. Recalling Figure 6, the cathode sits 
at the electrical axis of the gun surrounded by a Pierce electrode to provide initial 
focusing as the beam is launched [4], The three electrodes are electrically isolated. 
The Pierce electrode and cathode are maintained at large negative potential and the 
final electrode is at ground. The ratio of voltages between the Pierce electrode and fo¬ 
cusing electrode is step-wise variable by choosing its potential from steps on a voltage 
divider seen in Figure 61. 

In operation at Stanford, the DC gun was used in thermionic mode. By 
energizing a resistive heater behind the cathode surface, the cathode temperature was 
increased (T > 1000 K) such that electrons were “boiled” off the surface continuously. 
This is an extremely reliable method for generating high currents, but the current 
is continuous and thus not ideal for RF accelerator applications. To pulse the beam 
a secondary held is generated between the grid and the cathode surface opposing 
the accelerating fields in the gun structure. Since the grid is much closer to the 
cathode surface, the voltage on the grid can be less than 100 V and hold off the gun 
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Figure 61. The DC gun voltage divider showing the resistive ladder on the right 
connecting ground to the cathode potential. The white wire coming off to the left is 
an adjustable connection to the focusing electrode. 


fields. By switching the grid voltage, the electron beam can be pulsed. Unfortunately, 
the switching cannot be performed arbitrarily fast, so generating temporally short, 
high charge bunches requires additional tools. In the case of the Stanford FEL, a 
sub-harmonic buncher cavity was used—the beam passes through the cavity during 
a field zero crossing such that the head of the beam is decelerated and the tail is 
accelerated compressing the beam longitudinally. 

For NPS purposes, we desired to use the gun in a combined photo- 
thermal mode: heating the cathode to increase the electron kinetic energy, but not to 
the point of liberating them from the cathode surface, and providing the liberating 
energy via a drive laser [17], [18]. The drive laser also acts as the gating mechanism, 
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allowing for much higher repetition rates and shorter pulse durations. Initial exper¬ 
iments indicate this approach is sound, and “first beam” in the NPS Beam Physics 
Lab facility was achieved 17 April 2009 as captured in the inset of Figure 62. The 
drive laser for these experiments was the first generation version of the minimal drive 
laser design presented earlier. From this system, the first attempts at generating 
beam showed significant discrepancies between the expected charge liberated and the 
charge collected at the dump. 



Figure 62. NPS staff and researchers celebrate “first beam” of the NPS Beam Physics 
Lab, 17 April 2009. In the inset is a plot of the current captured from this first 
measurement. 


The NPS DC gun beamline, shown in Figure 63, consists of a DC gun, 
a four-way cross for laser entry and mirror insertion, a solenoid for electron beam 
focusing, a large six-way cross for vacuum pumping, and an instrument section with 
Faraday cups and a residual gas analyzer. The DC gun can operate with potentials 
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up to 250 kV, provided by a high voltage power supply. The Faraday cup section 
allows for crude beam centering as the two metal blocks are electrically isolated and 
conhgnred similar to Figure 64, with both plates being connected to an oscilloscope 
for integration of the two signals. 



Pumping 


Faraday Cups 

RGA \ Solenoid 


Corona 

Shield 


DC Gun 


Figure 63. The DC gun diagnostic line for photo-thermal electron beam generation. 


Using the best drawings available (the drawings provided with the Stan¬ 
ford gun show a third stage that was apparently never built), the internal dimensions 
and geometries were input in Poisson and solved for static electric fields. The desired 
operating practice is to hold the Pierce electrode and cathode at -240 kV, the focusing 
electrode at -140 kV, and the final electrode at ground potential. The geometry and 
resulting eqnipotential lines are shown in Figure 65. In actual operation, we were un¬ 
able to condition the gun to its design voltage. However, because the voltage divider 
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Figure 64. The Faraday cup design in the diagnostic line also allows for determination 
of beam centering during transport. The first Faraday cup consists of a copper plate 
with a central hole approximately 1 cm in diameter. The second Faraday cup is a 
copper block approximately 3 cm in diameter. Prior to taking measurements, the 
voltage signal on the second Faraday cup is maximized and the voltage signal from 
the first Faraday cup is minimized, ensuring the beam is properly centered as it passes 
into the diagnostic station. 


only allows for stepwise settings between the electrodes, the ratio of cathode/focusing 
voltages remain constant regardless of the operating voltage at the time. By scaling 
the field map in the simulation software, we can model any potential value for a given 
electrode voltage ratio. 

The initial diagnostic beam line was planned for developing procedures 
and initial measurements of quantum efficiency and bunch charge. Using the au¬ 
thor’s experience in naval engineering, these were developed to be easily understood 
and transcribed to fleet language. Additionally, the experience developed in these 
experiments is directly applicable to experiments planned for the SRF gun. To deter¬ 
mine the optimal voltage ratio between the cathode and focusing electrode, various 
DC gun field maps were generated and used as accelerating gradient sections in GPT 
simulations. The beamline layout, as shown in Figure 63, was incorporated into the 
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z [cm] 

Figure 65. The Poisson model for the Stanford DC Gun. The left-most structure 
is the cathode and Pierce cone. The middle electrode is shaped to provide beam 
focusing as the beam is accelerated toward the grounded electrode (final structure on 
the right). The pink lines mark equipotential lines. 


GPT model and various bunch charges propagated to determine the optimal gun 
setup for performing the desired measurements. The results of these simulation runs 
are shown in Figures 66 and 67. 

The simulation runs show us important values as the beam evolves 
moving through the DC gun structure. In Figure 66, we see that for the five electrode 
configurations (where the first number indicates the percent of the total voltage drop 
between the cathode and focusing electrode and the second number is the percent 
voltage between the focusing electrode and the ground electrode) the beam gains 
energy as expected. In these runs the maximum voltage potential was established as 
-160 kV, commensurate with the conditioned status of the system at the time. These 
simulations also reflect relevant parameters expected in actual experiment—a 5 nC 
bunch charge liberated in a 4 ns pulse (due to pulse length of the Minilite II drive 
laser) from a 4 mm radius spot size at the cathode and a uniform beam distribution. 

If we now focus our attention on Figure 67, we see that it is possible to 
transport the entire bunch charge as described through the gun structure. Of primary 
interest is the effect of space charge on the bunch. In the 25-75 case, the bunch 
drastically expands through the focusing electrode (low focusing, and still a lower 
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Electrode Bias Study - Beam Energy 
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Figure 66. As the voltage ratio between the focusing electrode and cathode are 
changed, the beam falls through a different potential field. The plots show the effect 
of changing this ratio for a constant maximum voltage of -160 kV, illuminated laser 
spot, and bunch charge. 


energy beam) and exits the gun at nearly the full aperture of the transport system. 
From this graph, we can see that it is desirable to have as much of the potential drop 
across the cathode-focusing electrode gap as possible. In realistic terms, the voltage 
divider is only capable of delivering approximately 57% of the potential between the 
cathode and focusing electrode, putting the expected beam envelope close to that of 
the 60-40 case. We find this setup acceptable as it transports a beam from the gun 
diverging at approximately the same rate as all cases considered with the exception 
of the 25-75 case. The location of the normal conducting solenoid, immediately after 
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Figure 67. As the voltage ratio between the focusing electrode and cathode are 
changed, the beam envelope changes. The plots show the effect of changing this ratio 
for a constant maximum voltage of -160kV, illuminated laser spot, and bunch charge. 


an isolating gate valve, is as close to gun exit as practicable, and we can expect to 
be able to focus the majority of the beam to the diagnostics further down the beam 

line. 

Our first experiments focused on attempting to measure the cathode’s 
quantum efficiency. Quantum efficiency is a ratio of the number of electrons gener¬ 
ated per incident photon. To measure this quantity, we integrate the voltage signal 
from the Faraday cups that measure the current passing through a known resistance 
providing the bunch charge. Simultaneously, we measure the pulse energy of the laser 
pulse incident on the cathode surface and assume that the transmission properties of 
the transport line remain constant and that the full pulse energy reflected from the 
laser mirror is incident on the cathode surface. As an example using actual data (266 
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nm UV wavelength, pulse energy of 170 n J, 50% transmission, Vfci — 1-68 V, V F C 2 
= 7.3 V, cr FC i = 7.18 ns, cr FC2 = 12 ns, and 77 = 10 kh2, and assuming the voltage 
pulses are triangular), where “FC” indicates the Faraday cup the measurement was 


taken from, we find the quantum efficiency of the system to be 
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As the first measurements were made, the quantum efficiencies and 
bunch charges were much lower than expected for a photothermal cathode [36], [37], 
[38]. Figures 68 and 69 show plots of quantum efficiency versus heater voltage and 
quantum efficiency decay over time for the Stanford DC gun. 

To compile the data shown in Figure 68, the cathode was allowed to 
stabilize at a specified heater voltage for a minimum of 10 minutes prior to data col¬ 
lection. Beam transport was optimized (maximizing V F C 2 + Vfci) f° r a constant gun 
voltage of -145 kV, which was slightly below the maximum conditioned voltage at 
the time. Bunch charge and pulse energy were recorded at 30 second intervals for 3 
minutes and averaged to compile a quantum efficiency for a particular heater voltage. 
Laser pulse energy was constant at approximately 170 /jJ throughout the experiment. 
The illuminated area on the cathode is estimated to be approximately 1 cm in diam¬ 
eter which would illuminate the entire active area of the cathode. In reviewing the 
data and comparing to estimates and measurements of quantum efficiency for com¬ 
parable cathodes, our values were found to be at least two orders of magnitude below 
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Figure 68. Quantum efficiency of the DC gun as heater voltage is increased. 


expectations. Prior to performing further experiments, we focused on attempting to 
increase laser energy transport and further conditioning of the cathode and gun. 

After application of the improvement tasks, the DC gun photothermal 
cathode finally began to approach quantum efficiencies measured in other systems. 
Bunch charges were measured commensurate with the improved quantum efficiency 
and higher pulse energy (now delivering approximately 250 /iJ to the cathode surface). 
Additionally, the 10 kfl resistors previously used to measure voltage on the Faraday 
cups were replaced with 50 D resistors for better signal matching to the oscilloscope. 
Figure 69 shows a quantum efficiency decay plot for the cathode after performing a 
brief surge to full heater voltage (5.0 V) for 30 minutes and then allowing the cathode 
to cool for 30 minutes prior to taking data. We see the expected exponential decay, 
but the cathode did not reach a minimum value as we expected. Our expectation, 
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Figure 69. Quantum efficiency decay of the DC gun cathode held at constant voltage 
below the replenishment voltage. 


given the age and history of the cathode, was that the decay plots should start to reach 
a minimum fairly quickly after the heater was turned off. It is possible that these type 
dispenser cathodes are much more robust than we anticipated. Unfortunately, due 
to external factors, we have not been able to revisit this to extend the measurement 
period in hopes of observing a minimum. 

As experiments were performed using the DC gun, numerous issues 
were identified that may add uncertainty to the results obtained thus far. With 
the drive laser wavelength in the UV, we have been unable to determine where the 
electrons are being generated. When using the UV beam (266 nm), all known metals 
in the region of the cathode (cathode surface, grid, Pierce cone, all other electrodes) 
will photoemit. Removing the second doubler crystal allows for operating the laser in 
the green region (532 nm) which is below the work function of copper. Unfortunately, 
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we were unable to validate laser alignment at 532 nm due to laser crystal issues and 
could not rule out photoemission from the copper structures in the gun. 

Transporting the electron beam a significant distance allows for beam 
expansion and possible losses in transport, especially at low energies. Most quan¬ 
tum efficiency studies are performed with a small cathode-anode gap and extremely 
short transport distances to ensure capturing all generated charge in the collection 
mechanism. With the DC gun, this is problematic due to the high voltages involved. 
However, in this device, it is possible to measure the replenishment current flowing 
into the cathode region. To improve our ability to measure where the electrons in the 
beam were being generated, the Eimac cathode’s separate (and isolated) electrical 
connections, shown in Figure 70, have been instrumented with fast current trans¬ 
formers to capture the replenishment charge flowing into the various parts (Pierce 
electrode, cathode, and grid). The concept has been verified at low voltage without 
the drive laser, but we have not been able to test at design gradients due to an ongoing 
radiation stand down at our NPS facilities. 

Since the Eimac cathode used in the DC gun system is the one that 
was shipped installed in the gun from Stanford and was exposed to air for more than 
two years before being placed under vacuum, it is possible that the cathode has been 
significantly contaminated. No attempt to reinitialize the cathode was made for the 
results presented here. A modified reinitialization was attempted later, but no signif¬ 
icant improvement was seen. To better understand the cathode and its performance, 
and possible effects on an SRF cavity, separate, non-beam generating experiments 
were conducted [41], To incorporate lessons learned and minimize nonessential vari¬ 
ables, a lower voltage (100 kV) test cell with a simplified cathode-anode geometry is 
being manufactured that will incorporate the improved measuring scheme and the up¬ 
graded drive laser system for testing various photothcrmal beam generating schemes. 
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Figure 70. The general layout of an Eimac electron source is shown. The picture 
on the left is an actual device, removed from its vacuum system. Connection to the 
the lab system is by conflat flange with copper seal. The related diagram on the 
right shows the internal wiring and isolation of the cathode sourcing, grid bias, and 
heater systems. Not shown is the Pierce electrode, which is part of the gun internal 
structure—it would be located very near the grid but still electrically isolated. After 
[39], [40], 
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5. Booster Configuration 

The first injector system modeled was the booster configuration. In this con¬ 
figuration, it was intended that the source beam would be generated using the former 
Stanford DC gun at an energy of approximately 150-200 keV. The booster would 
provide an additional 1.2 MeV for merging into the primary accelerator. Of primary 
difference with the previous DC gun simulations, the beam pulse length is greatly 
decreased. In the previous case, a 4 ns pulse was used, but this corresponds to two 
complete RF cycles of the 500 MHz booster and is not ideal. To correspond to com¬ 
mercially available UV drive laser specifications, a “long” drive laser pulse of 40 ps 
is used with a reduced bunch charge (1 nC). There are drive lasers available with 
shorter pulses and higher repitition rates, however this pulse length corresponds to 
the drive laser being used in the current 500 MHz gun experiments. 

Since the booster is of primary concern for this study, we assume the input 
beam is uniform, monoenergetic, and has an initial emittance of zero to remove any 
effects that might be present from the electron source. The beam is initiated at the 
closest point of connection from the booster to the gun section. Space charge forces 
are not applied until after the beam reaches this earliest connection point (62.7 cm 
before the booster cavity nose cone). The simulated beam line is terminated at a 
point approximately 3.0 m from the booster nose cone to compare results from both 
Parmela and GPT, and to correspond to the beam line as originally proposed. These 
simulations also serve to compare the two code results on similar beam lines. 
a. Acceptance 

Acceptance as defined in the booster configuration under consideration 
will be defined as the cavity phases for which an injected beam can successfully be 
propagated through the cavity. Obviously, there are injection techniques that can 
be used to adjust the cavity acceptance, however these methods can only change the 
acceptance by a little compared to the overall acceptance of the cavity. To simulate the 
cavity acceptance, multiple runs, using the same beam parameters and only varying 
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the cavity phase, we can obtain general plots for the optimal phasing of the booster 
cavity for charge transport, beam size, emittance, etc. Figure 71 shows plots of the 
cavity acceptance for a 200 kV input beam as described previously. 
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Figure 71. Acceptance of a the NPS 500 MHz booster cavity for a 200 kV input beam. 
The green shaded area indicates “acceptable” beam performance characteristics while 
red indicates beam that has transport problems through the cavity. 


In Figure 71, the “acceptable” cavity phases have been highlighted in 
green. To come to this conclusion about “acceptability,” there are two main data 
points that are considered. First is whether the beam is transported to the expected 
z position, in this case the center of the superconducting solenoid just after the booster 
cavity at a z position of about 0.79 m. Additionally, there should be a significant 
population of particles that are transported to the desired z location. The simulation 
software will calculate parameters for a single particle if it can successfully transport 
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only one particle through the desired geometry; however, these type results are not 
useful. From Figure 71, we see that the cavity, at best, transports about 22% of the 
initial population (of 30,000 particles) to the desired z position. Within the phase 
range of —40° < <f> < 75°, essentially the entire bunch that makes it through the drift 
to the cavity and through the iris of the nose cone and can be propagated to the first 
focusing element of the beamline. Other than the energy plot plotted in Figure 71 
that reflects the phase of the cavity, we do not yet have any information about the 
beam itself. 

b. Beam Transport 

To obtain data about the beam itself, it is necessary to transport the 
beam through the desired beamline and observe its performance characteristics. In 
our case, we are primarily interested in the beam energy, radius, and emittance. Once 
the beam has left the accelerating structures, its average energy does not change un¬ 
less beam scraping occurs. The internal beam energy will change as space charge 
accelerates the head and decelerates the tail, resulting in energy spread within the 
beam. With respect to the beam’s radius, there are two primary measurements of 
interest. The first is the average beam radius which represents the primary beam that 
is transported. The second is the maximum radius which takes into account the par¬ 
ticles that have been ejected from the core beam and make up the beam halo. While 
the beam size can be adjusted after the primary accelerator fairly easily, emittance 
after the gun/injector, in general, cannot. Unfortunately, without discarding portions 
of the beam, once emittance growth has occurred, it usually cannot be reduced again. 
It, therefore, behooves us to choose a launch phase where the beam does not have 
excessive emittance. Figure 72 demonstrates the effect of cavity phase on both beam 
radius and emittance at the booster exit. Recalling that the peak of the energy plot 
is at 0°, the ideal phase for radius and emittance will probably be less than 20°. 
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Beam Parameters of Booster 
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Figure 72. Emittance and beam radius dependence on cavity phase of the NPS 
booster cavity. The sharp increase at the far right is due to poor transport at that 
phase value. 


Figure 73 shows the 200 kV beam evolution as it is propagated through 
our notional beam line from initiation, through the booster cavity, the superconduct¬ 
ing solenoid, the normal conducting solenoid to a point approximating where the 
merger to an accelerator structure or other diagnostic would be placed (approxi¬ 
mately 4.8 m from the beam launch). Transport parameters (solenoid strengths) 
were determined using the methodology discussed earlier, attempting to preserve as 
much of the beam as possible. After optimization through the beam line, the param¬ 
eters of interest were found to be: </> = 29.41°, Bsc-pk = 280 mT, B NC - P k = 28.7 mT, 
where “SC” indicates the superconducting solenoid and “NC” indicates the normal 
conducting solenoid. The transport envelope using these optimized values is shown 
in Figure 73. We see that the beam radius increases significantly from initiation to 
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the booster nose cone iris, where the beam loses more than 1/3 of the charge. From 
this point, we see the remaining beam is fairly successfully transported to a beam 
focus at the target. The divergence of the maximum radius from the average ra¬ 
dius is an indication that the combination of the booster cavity and superconducting 
solenoid generate significant halo that impact the beam pipe just before the normal 
conducting solenoid. Propagation from the normal conducting solenoid to the target 
is uneventful other than generation of halo from the back of the beam where it is 
overfocused by the normal conducting solenoid at about the z = 4.0 m point. 


Booster Transport -- 200 kV Input 
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Figure 73. Beam transport envelope for the NPS booster cavity with a 200 kV input 
beam. 



In an attempt to increase beam transport and to counter the effects of 
space charge in the initial drift prior to the booster cavity, a second normal conducting 
solenoid with the same field map as the beam line version is placed between 20 and 
30 cm from the gun connection (which is also 20-30 cm before the booster cavity 
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nose cone). We optimize for the solenoid positioning and peak field strength using 
the transport parameters (0, Bgc-pk, and -Bjvc-pfc) from before for comparison. The 
new transport envelopes are shown in Figure 74 and show significant improvement 
in booster cavity performance. Optimization places the center of the input solenoid 
at 22.7 cm from the beam start with a peak held strength of 66.9 rnT. This position 
and peak held result in a maximum beam radius of about 1.3 cm as it passes through 
the booster cavity nose cone. From the radius plots and loss of particles shown in 
in the particle plot, it now appears that the superconducting solenoid is overfocusing 
the beam and that further optimization would result in nearly full transport of the 
beam. Regardless, the addition of this one element increased beam transport by a 
factor of more than four. 


Booster Transport -- 200 kV Input, Input Solenoid 
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Figure 74. Beam transport envelope for the NPS booster cavity with a 200 kV input 
beam and the addition of a normal conducting solenoid prior to booster injection. 
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The impact of adding this solenoid is important to the transverse ernit- 
tance. Figure 75 shows the two optimized beams overlayed for comparison; the only 
difference is the addition of the pre-booster solenoid. The blue plot represents the 
addition of the starting solenoid, which the red plot represents only a drift until the 
booster cavity. From these plots, we see the solenoid initially increases emittance, 
but the emittance growth seen in the red plot going through the booster cavity is 
mitigated by the effect of the starting solenoid. This emittance improvement (al¬ 
most a factor of 2 over the no-focus injected beam) is maintained almost through the 
rest of the entire transport line. The non-compensated beam having an emittance of 
66.3 n m compared to the solenoid compensated beam emittance of 40.8 //m. Neither 
emittance is particularly noteworthy, but considering the less than ideal initiation, 
the fact that it is not larger is of note. 

From this study, it appears that the NPS cavity is suitable as a booster, 
but either the input drift should be minimized or some emittance compensation device 
(such as a solenoid) should be included in the design prior to the booster. However, 
an injection energy of only 1.4 MeV may present a challenging merge problem for the 
main accelerator. This issue is not taken up in this study, but left for future work. 
We now depart the booster concept briefly to examine the use of the NPS 500 MHz 
cavity in the gun configuration. 
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Booster -- Emittance Comparison 



Figure 75. The addition of a solenoid prior to the booster results in higher emittance 
at the gun, but lower emittance at the end of the transport line. The final emittance 
value for the drift only simulation is a significant underestimate as the gun nose cone 
has filtered out the high emittance particles and the transported final bunch charge 
is significantly lower than with the addition of the input solenoid. 
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6. Gun Configuration 

Insertion of a cathode stalk into the booster allows the same cavity to be used 
as a gun. The type of electron source is dependent upon the type of cathode installed. 
The NPS SRF gun has been designed with multiple cathode types in mind, including 
field emitter arrays and photoemitter cathodes (metal and semiconducting). Turn-on 
fields for gated diamond field emitter arrays have been described as anywhere from 
1 MV/m [23] to 15 MV/m [7]. As these are fairly new devices that have typically 
only been operated in test stands with small voltages and small cathode-anode gaps, 
their performance in realistic gun environments is unknown. If we assume the turn¬ 
on voltage in a real application will be bounded by these values, and assume 15 
MV/m is necessary for full current development from a field emitter tip, then we 
see from Figure 17 that fields are acceptable for field emitter array cathode use for 
all retraction positions (relative to the nose cone) forward of 6.5 mm. The 6.5 mm 
retraction position represents the maximal retraction for ensuring turn-on across a 
1 cm diameter array cathode as electric field values drop significantly as they are 
measured radially from the electric axis. For photoelectric cathodes, a drive laser is 
required along with the laser port into the beam line and a mirror to direct the laser 
beam to the cathode. 

For both simulation packages, the beam line was modeled (elements placed) 
based on the system installed at Niowave. Figure 76 shows the layout and provides 
the as-built measurements. The only difference in the model input files is that the 
distance from the nose cone to the cathode surface has to be accounted for as all 
measurements in the beam line use the nose cone as the reference point. 

a. Acceptance 

Using a standard 1 nC, uniform electron beam, we first determine the 
expected acceptance of the gun for various cathode positions. Figures 77 and 78 show 
these plots and pertinent values as a function of the launch phase of the electron beam. 
These figures are fairly dense with information, so we take some time to analyze them. 
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Figure 76. The NPS 500 MHz gun diagnostic line has most of the instrumentation 
necessary to fully characterize the expected beams from the gun. Numerous steering 
coils, diagnostic pots (each of which includes a Faraday cup, an yttrium aluminum 
garnet crystal for beam imaging, and space for future pepper pots for omittance mea¬ 
surements) and quadrupoles are installed for ensuring proper transport. A dipole 
magnet and spectrometer leg are included for beam energy and energy spread mea¬ 
surements as well. 


Looking at the top of Figure 77, we see that all of the cathode positions have about 
the same phase acceptance window (85° < </> < 225°) for the cavity voltage applied 
(1.24 MV). We also see from this plot that the focusing at the cathode provided by 
retracting the cathode behind the nose cone is critical—until the cathode is retracted 
about 5 mm, full transport of the bunch charge is impossible. Once the cathode is 
retracted approximately 6.5 mm, full transport of a 1 nC bunch becomes possible 
at nearly all phase angles. Moving to the energy plots in the bottom of Figure 77, 
we see that as the cathode is retracted, the peak of the energy plots are shifting 
slightly earlier. This can be accounted for by the additional time of flight required 
by the bunch to reach the high held regions in the gap. Peak average beam energy 
is affected by the cathode position, but only varies by about 10 keV over the 9 mm 
of cathode travel. The spikes in two of the plots near the turn off phase are artifacts 
of how GPT counts particles. Since they are far beyond the peak energy where 
we plan to operate the cavity, they can be ignored. The radius plots in the top of 
Figure 78 provide correlation with the transport plots above them. The ragged lines 
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are indicative of particles being lost during transport. As the average radius plots 
remain fairly smooth for retractions positions beyond 5 mm, it can be assumed that 
the losses in these beams are small compared to the entire bunch charge. The spike 
in radius seen in the 9 mm retraction position in this subplot can be attributed to 
being beyond the acceptance cutoff of the cavity and ignored. In the bottom plot 
of Figure 78, we see that the launch phase chosen has an impact on the transverse 
omittance expected out of the gun as well. We also see that it is desirable to operate 
the gun in the region where the fields are building as this results in generally lower 
omittance values. From this plot, we further see that the cathode position again plays 
a role in omittance determination. For the 0 and 2.5 mm retraction plots, the ragged 
lines indicate again that particles are being lost—that cavity transport is acting as 
an omittance filter removing particles contributing to higher emittance as they strike 
the beam tube or fail to exit the cavity. If these plots followed the general trend of 
the other retraction values prior to transport losses, it is likely they would result in 
higher emittance values. 
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Figure 77. Acceptance plots of the NPS gun are plotted for five cathode positions. 
The positions are labeled with respect to their retraction behind the nose cone in mm 
(i.e., “0.0” is the cathode face flush with the nose cone and “0.9” is 9 mm retracted 
behind the nose cone). The top plot is the general transmission of charge through the 
cavity with the z position of the final screen passed included for clarification. The 
bottom plot is the output final energy versus launch phase plots. 
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Radius - Phase Curves 
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Figure 78. Radius and emittance dependence of the NPS gun on phase plotted for 
five cathode positions. The positions are labeled with respect to their retraction 
behind the nose cone in mm (i.e., “0.0” is the cathode face flush with the nose cone 
and “0.9” is 9 mm retracted behind the nose cone). The top plot is maximum and 
average radius at the final output screen as a function of phase. The bottom plot is 
transverse emittance as a function of launch phase. 


126 





























To examine the impact of the pulse shape on the gun performance, 
we choose a single cathode position for analysis. The 6.5 mm retraction position 
was chosen as we know from Superfish studies that the entire cathode face remains 
above the expected turn-on voltage for held emitter cathodes, broadening the overall 
applicability of the beam studies performed. Using these more advanced beam shapes, 
we plot in Figures 79 and 80 the acceptance plots for a single cathode position to 
examine the effects of more realistic beams on the gun acceptance. 

From Figure 79, we can see that the acceptance window for each of the 
beam types is essentially the same. Recalling that for the window to be valid, both the 
z 0 ut position should be greater than 18 cm and the transmission should be close to 1.0. 
We find the cavity phase window for each beam type to be about 85° < (j) < 225°. 
There is little difference in the energy plots. The difference in the rising plots is 
merely a sign of the S(p used in generating the plot (5°) indicating the actual leading 
edge for transport through the cavity is probably near 83° and the difference between 
the beam types is about 1°. Looking to the radius plots in Figure 80, we see that 
all three beams show similar behavior. The average beam size for the Gaussian and 
asymmetric beams appears less than the uniform beam due to the Gaussian shape 
radially—there is more charge in the core of the Gaussian and asymmetric beams at 
the outset. The overall beam sizes, using the maximum radius all follow very similar 
plots. Likewise, there is little difference in the beam omittances, shown in Figure 
80, other than an apparent slight advantage for the uniform beam for earlier phase 
delays, but it only a difference on the order of 1 fim. From these plots, it appears that 
the NPS cavity is fairly tolerant of the beam shape and no major impacts differing 
from the simplest uniform beam occur when more complicated beams are considered. 
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Beam Type Effect on Acceptance 



0.00E+00 5.00E+01 1.00E+02 1.50E+02 2.00E+02 2.50E+02 3.00E+02 3.50E+02 4.00E+02 

Cavity Phase [deg] 


Energy - Phase Curves 



Cavity Phase [deg] 


Figure 79. Various parameters of the NPS gun are plotted for three different beam 
types: Uniform, Gaussian, and asymmetric. The top plot is the general transmission 
of charge through the cavity with the z position of the final screen passed included for 
clarification. The bottom plot is the output final energy versus launch phase plots. 
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Radius - Phase Curves 




Figure 80. Radius and emittance dependence of the NPS gun on phase plotted for 
three different beam types: Uniform, Gaussian, and asymmetric. The top plot is 
maximum and average radius at the final output screen as a function of phase. The 
bottom plot is transverse emittance as a function of launch phase. 
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b. 


Beam Parameters 


Now that the gun has been studied, further characterization requires 
propagation of representative beams to determine the beam characteristics that can 
be expected from the gun. Two primary optimization regimes are considered— 
optimizing for the smallest spot size at a given point (assumed to be a dipole or 
merger entrance 3.66 m from the nose cone in subsequent simulations) or for the 
minimum omittance at the same location. From the plots shown in Figures 77 and 
79, we can expect the optimal launch phase to be near 150° to maximize energy for 
a cathode retracted 6.5 mm for both the omittance and spot size cases. Proceeding 
with the optimization scheme described previously, we see the final optimized beam 
envelopes in Figure 81 and the omittance evolution in Figure 82. 

Looking at the plots in Figure 81, there is surprising similarity between 
the envelopes regardless of the optimization goal. The spot size optimization plots 
appear to delay halo formation in all the beam shapes until near the simulation termi¬ 
nation except for the asymmetric case. In the emittance optimization, halo generation 
takes place for the uniform beam just before termination and the asymmetric beam 
about 60 cm before termination. From the beam visualization software, it can be seen 
that the halo formation is being generated by the lower energy tail of the beam being 
over focused clue by the normal conducting solenoid. Of the three beam types, only 
the asymmetric beam appeared smaller (using the maximum radius as the measure 
of interest) when using the emittance optimization routine, but the terminal primary 
beam size was still smaller in the spot size optimization. 

Beam emittance is a critical parameter for beam transport. In general, 
the smaller the number can be made, the easier transport becomes. Again, there 
was little difference between the two optimization routines when comparing beam to 
beam. The non-uniform beams were clearly better beams for transport arriving at the 
simulation termination with factors of 3 better (for the Gaussian beam) and 4 better 
(for the asymmetric case). From a prediction stand point this bodes well for using 
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the uniform beam as the primary simulation shape as it will tend to overestimate the 
beam cmittance. Since the simulations do not take into account surface imperfections 
of the cathode, differing quantum efficiencies along the cathode surface, and were all 
started as cold beams (no initial energy or emittance), the uniform beam is probably 
the best the experimentalist can hope for in terms of predicted performance. 

By the time the beams reached simulation termination, they had expe¬ 
rienced significant longitudinal growth in all cases. For injection into another accel¬ 
erating structure, we would like to have the beam fairly compact in order to ensure 
the beam does not develop any unwanted energy spread by sampling the accelerating 
holds too far out of phase from desired. To achieve this, we can shorten the transport 
line, however the line considered is probably close to a minimal design already. A 
buncher cavity could be installed to compress the bunch by passing the beam through 
an RF cavity near an RF zero crossing so that the head and tail of the beam receive 
kicks that would tend to compress the beam longitudinally. This was how the beam 
was compressed in the original Stanford FEL injector, however this scheme requires 
additional drift space for the compression to take place prior to injection into the main 
linear accelerator, expanding the total system length. The more compact solution is 
to add a booster to the gun and accelerate the bunch to a higher energy where space 
charge effects affecting the size of the bunch can be minimized. 


131 



Gun Transport -- Spot Size Optimization 



z[m] 


Gun Transport -- Emittance Optimization 



z[m] 


Figure 81. Optimized beam transport envelopes for uniform, Gaussian, and asym¬ 
metric beams initiated from the NPS gun. The top plot is based on optimizing for 
minimum spot size. The bottom plot is based on optimization for minimum emit¬ 
tance. 
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Gun Transport -- Spot Size Optimization 
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Figure 82. Optimized beam emittance evolution plots for uniform, Gaussian, and 
asymmetric beams initiated from the NPS gun. The top plot is based on optimizing 
for minimum spot size. The bottom plot is based on optimization for minimum 
emittance. 
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7. Injector Proposal 

During the exercise of predicting the performance of the NPS 500 MHz gun, it 
became apparent that with this gun design, it was of utmost importance to increase 
the beam energy as quickly as possible in order to minimize the space charge effects 
experienced by the electron beam. There are many booster cavity designs already 
developed. However, considering that all fixtures and designs for manufacturing the 
quarter wave structure have already been developed, the author proposes and analyzes 
a similar end-to-end structure be used as a booster. Figure 83 provides a Superfish 
model of such a cavity. 



Figure 83. The first iteration of the double bell quarter wave booster cavity. This is a 
reflection of the NPS 500 MHz quarter wave gun, with simplified nose cone structures, 
taken along the end plane of the cavity. The fundamental frequency is 500.30 MHz 
and the pink lines follow electric held lines. 


a. Cavity Development 

To develop this booster structure based upon the NPS gun/booster ge¬ 
ometry, there are two primary design goals: accelerate the beam generated from the 
NPS 500 MHz gun to an energy > 3 MeV for merge properties with the main accelera¬ 
tor and do so in as compact a package as possible at the same fundamental frequency. 
To accomplish these goals, some assumptions were made to aid the analysis. The first 
assumption was that the coaxial coupler currently installed in the NPS gun would be 
replaced with a coaxial coupler perpendicular to the beam line as recommended in 
the previous chapter. This provides RF power to the NPS gun. Second, the entire 
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assembly of gun and booster would be immersed in a combined cryomodule with the 
superconducting solenoids at the exit of the gun and the booster immersed in liquid 
helium to ensure they are below T c . Third, to minimize additional drifts, RF power 
is assumed to be applied to the booster cavity via intracavity antennae. No modifi¬ 
cations to the actual NPS gun geometry will be modeled; only basic modifications to 
the booster cavity were modeled. 

As a first attempt at generating a booster cavity, the NPS gun structure 
(minus the cathode assembly) was simply reflected about the anode side of the cavity, 
as shown in Figure 83. The booster cavity was then modeled in Superfish, using the 
same maximum electric and magnetic field parameters as were used in designing the 
NPS gun structure. This reflection approach worked surprisingly well, resulting in 
a booster cavity with a fundamental frequency of 500.30 MHz. After adding the 
correct internal nose cone structure, the Superfish model showed negligible change in 
the fundamental frequency. The resulting accelerating fields and radial fields on axis 
and 1 cm off axis are plotted in Figure 84. Peak on-axis accelerating fields, E z , are 
found near the two nose cones and peak at 21.0 MV/m. The large dip in on-axis 
accelerating field in the center of the cavity was concerning as it is a drop of 48% 
from the peak fields. The radial fields indicate that the cavity will provide substantial 
focusing (7.5 MV/m) as the beam passes through the entrance nose cone with similar 
defocusing fields at the exit nose cone. The defocusing at the end of the cavity should 
be mitigated by the beam increasing its energy by about 2.61 MeV in passing through 
the booster. 

To mitigate the field dip in the “reflected cavity,” the cavity was short¬ 
ened by 5 cm. After shortening the cavity, the fundamental frequency decreased to 
485.53 MHz. To tune the cavity back to the desired fundamental, the cavity radius 
was decreased by 6.9 mm from the starting 12 cm radius, resulting in a fundamen¬ 
tal frequency of 500.14 MHz. Plotting the same on-axis and 1 cm off-axis electric 
fields, shown in Figure 85, we see success in decreasing the dip in E z to only 75% 
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of the peak fields. Peak E z has also increased to 22.2 MV/m in the on-axis case. 
For off-axis fields, peak values increased for E z and E r as well, to 24.0 MV/m and 
7.9 MV/m respectively. Unfortunately, shortening the cavity results in a decreased 
integrated gap voltage of about 2.27 MV (a decrease of 13% over the longer cavity) 
and increased increased radial electric fields. It is unclear at this point whether it is 
worth mitigating the dip at the expense of beam energy. 
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Electromagnetic field data from file B00STER2.SF 
Problem title line 1: 500MHz NPS Electron/Gun Cavity 



Electromagnetic field data from file B00STER2.SF 
Problem title line 1: 500MHz NPS Electron/Gun Cavity 



Figure 84. Electric fields along the electric axis of the long, double QW booster cavity 
(top) and a line displaced 1 cm radially (bottom). The red line is E z and the green 
line is E r . 
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Electromagnetic field data from file B00STER2.SF 
Problem title line 1: 500MHz NPS Electron/Gun Cavity 



Z (cm) 

Electromagnetic field data from file B00STER2.SF 
Problem title line 1: 500MHz NPS Electron/Gun Cavity 



Z (cm) 


Figure 85. Electric fields along the electric axis of the short, double QW booster 
cavity (top) and a line displaced 1 cm radially (bottom). The red line is E z and the 
green line is E r . 
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b. 


Beam Modeling 


With satisfactory cavity field maps, both versions were imported into 
GPT to determine their performance on a beam generated by the NPS gun. Both 
systems were started with a uniform 1 nC electron bunch 40 ps in length and 5 mm in 
radius. The NPS gun cathode is placed at a 6.5 mm retracted position and the beam 
launched at the optimal phase for spot size or omittance depending on the booster 
optimization desired. The gun cavity and booster cavities are kept at maximum held 
level. The two booster cavities will be compared based on phase acceptance, final 
beam kinetic energy, bunch charge preservation, spot size at the notional merger 
entrance, and emittance preservation. 

We begin with the acceptance plots of the two booster cavities. After 
optimizing the gun for transport for spot size, the superconducting solenoid is opti¬ 
mized to force the beam to a waist at the exit plane of the booster cavity. Figures 
86-89 show the example acceptance surfaces of the the two cavities. In this case, 2D 
parameter surfaces are required for visualization as the phase delay and focus location 
from the gun’s superconducting solenoid are tightly coupled. 

For the long booster cavity, the emittance optimized acceptance plots 
are shown in Figures 86 and 87, where we have truncated the search volume to a 
small launch phase and gun solenoid strength (parameterized from 0 to 1) space where 
transport is expected. From the energy plot, we see that for optimal energy deposition 
into the beam, we should operate the booster at a phase delay of about 105°. The 
range of phase delays chosen (50° < (j) < 150°) all show excellent transmission as 
seen in the right hand plots. For the parameter space plotted, the gun solenoid has 
the most significant effect, with transmission decreasing rapidly above 0.38 (180 mT) 
and below 0.18 (86 mT). The lower left emittance surface seems to indicate that 
the emittance improves as the phase delay increases and the solenoid value decreases 
relative to the range considered. 
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When we look at the acceptance surfaces (spot size optimization) for 
the short booster cavity in Figures 88 and 89, we see that the effects of booster delay 
phase and solenoid focusing from the gun solenoid are fairly decoupled. Note that 
this plot is also only over a region of nearly full charge transport for purposes of 
optimization. If one considers the full range of phase delays, the surface becomes 
much more complicated, but as we are interested in transporting our beam at the 
greatest possible energy, we ignore these other areas. From these plots, we see that 
the booster should operate with a delay of about 100° from the cavity launch phase 
and the gun solenoid should have a peak held of about 150 mT. In general, both 
booster designs have equivalent phase delay acceptance ranges of 50° < (ft < 150° 
with respect to the gun launch phase and gun solenoid peak holds of between 100-200 
mT being satisfactory for either optimization scheme. 
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Figure 86. Acceptance plots for the long, double QW booster cavity, optimized for 
smallest transverse emittance at the booster exit, are shown. The upper plot shows 
average kinetic energy of the beam as the cavity phase delay is changed. The bottom 
plot shows a strong emittance dependence on solenoid strength for higher solenoid 
strengths. 
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Figure 87. Projections of the plots shown in Figure 86. 
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Figure 88. Acceptance plots for the short, double QW booster cavity, optimized for 
smallest spot size at the booster exit, are shown. The upper plot shows the average 
kinetic energy of the beam as the cavity phase delay is changed. The bottom plot 
shows a strong spot size dependence on solenoid strength. 
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Figure 89. Projections of the plots shown in Figure 88. 
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If we then continue the transport optimizations, we can obtain a final 
system for comparison to the pure gun and transport cases considered earlier. Figure 
90 shows the beam envelopes for the spot size optimized booster designs. With the 
compact nature of this booster configuration and assuming an additional 0.5 m from 
the booster solenoid for the cryomodnle containment, the injector cryomodule would 
terminate at about the z = 1.3 m. We allow for an additional 1.7 m propagation dis¬ 
tance similar to the propagation distance in the gun design to the normal conducting 
solenoid. We see that both the long and short booster cavities have an average radius 
of about 4 mm at the termination point, although both suffer from halo generation 
due to lower energy particles ejected at the beam focus. However, due to the dispar¬ 
ity between the average and maximum radius, the halo charge makes up little of the 
total bunch charge for both designs. Both designs achieved full transport of 1 nC for 
all optimizations. We also see in Figure 90 that the long booster results in a higher 
beam energy, by about 210 keV, which is about half of what was expected from the 
integral fields. Both cavities, when optimized, apparently run slight off crest. Both 
the long and short booster designs fall about 10% short of their maximum energies, 
but still above the desired design goaf of a 3 MeV beam at the exit, when optimized 
for the smallest spot size. 

Figure 91 shows the transverse emittance evolution for the injector de¬ 
signs. For both cases, the gun and booster cavity parameters were optimized to 
obtain the smallest transverse emittance at the termination point. Neither cavity 
significantly improves over the emittance values found from just the gun with both 
designs having final emittance values of 20-30 pm. The short booster cavity does ap¬ 
pear to have an advantage over the longer one as it does not add as much emittance 
in the process of acceleration. Of note is that neither booster greatly increased the 
transverse emittance, indicating they are still good candidates for boosters. Addi¬ 
tionally, we see that optimizing the injector system for transverse emittance causes 
both booster cavities to run further off crest, resulting in a decrease in final beam 
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Figure 90. Beam transport envelopes for the long and short double QW booster 
cavity designs. The NPS gun is optimized for smallest spot size at the gun exit with 
the cathode 6.5 mm retracted for both designs. 


energy (though still above the 3 MeV design goal) of about 100-200 keV with the 
longer booster cavity suffering more than the shorter cavity. If we compare these 
optimized omittance values with the transverse emittance values from the spot size 
optimizations ( ei ong = 13.5 //rri, e s h or t = 17.6 /ini) we see that there was no significant 
improvement over the spot size optimization. In fact, the long booster found a worse 
configuration for emittance. 

The final comparison points between the original gun design and the 
proposed injector are in longitudinal emittance and energy spread. Figure 92 shows 
the output transverse emittance trace space for the injector system. From these 
plots, we can observe that there is significant energy spread and bunch length change. 
The spot size optimized beams show 20%-30% increase in bunch length, while the 
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Figure 91. Beam emittance evolution for the long and short double, QW booster 
cavity designs. The NPS gun is optimized for best emittance at the gun exit with the 
cathode 6.5 mm retracted for both designs. 


emittance optimized beams show a decrease in bunch length of about 10%. Table 
5 shows the numerical values related to Figure 92 for the four injector and two gun 
optimization cases. 


Table 5. Longitudinal omittances for the NPS gun and proposed injector designs. 


Design 

Gun 

Gun 

Long booster 
Long booster 
Short booster 
Short booster 


Optimized For 
spot size 
emittance 
spot size 
emittance 
spot size 
emittance 


e 2 [eV s 

i.8ii x nr 8 

1.054 x 10" 8 
1.764 x 10“ 7 
7.672 x 10" 8 
1.688 x 10" 7 
7.053 x 10" 8 


Energy Spread 
1 
1 

7 

8 
9 
6 
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Recalling that the gun cavity parameters were optimized separately 
from the booster to isolate the booster contributions, it is definitely possible that a 
combined optimal configuration between the gun and booster could be found that 
would decrease the omittance values further. No clear “winner” was seen in this 
study, however with a higher output energy and comparable optimized spot size 
and omittances (transverse and longitudinal), the longer cavity may provide slightly 
greater benefit in an injector design. 
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Figure 92. Longitudinal beam emittance plots for the long and short double QW 
booster cavity designs. The upper plots are for the long booster, and the lower plots 
are for the short booster. 


149 






B. SRF GUN/BOOSTER EXPERIMENTS 

Simulation provides many insights into how systems work and can provide 
the experimentalist a platform upon which proposed experiments may be tried with¬ 
out fear of damaging expensive equipment. Simulations also provide the researcher 
the opportunity to explore larger data sets than possible with actual experiments. 
However, there is no substitute for developing hard data by experimentation with a 
physical system. If the experimental results compare favorably with simulation data, 
it provides confidence in both—that the simulation parameters are representative of 
the physical system and that the experiment is behaving as predicted by simulation. 

As part of this philosophy, we have already presented simulation results for 
many possible experiments with the NPS gun/booster. In this portion of the chapter, 
we will discuss the actual experiments conducted during the development and con¬ 
struction of the gun/booster as well as the experimental results from the first beam 
testing of the first SRF gun in the United States [42], 

1. Resonant Modes 

The gun/booster system is fundamentally a resonant cavity for storing RF 
energy and coupling that energy into a beam of electrons. We previously discussed 
how the cavity shape was adjusted relative to the maximum fields in the cavity, but 
the situation is actually a bit more complicated. Every cavity will be able to support 
a sequence of discrete frequencies of RF, usually associated with higher-order modes. 
For a cavity intended to operate a the fundamental frequency, such as the NPS 
gun/booster, it is desirable to have any other higher-order modes be well separated 
from the fundamental frequency. 

The resonant frequency is one of the first characteristics that can easily be 
tested, it was also one of the first to be simulated. Two separate software packages 
were used to simulate the cavity and determine the resonant modes: Superfish and 
COMSOL Multiphysics. 
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a. Superfish Model 

To evaluate the cavity and the separation of resonant frequencies, we 
started with Superfish. Recall from the previous section that Superfish uses the 
meshed area to calculate a fictitious magnetic current density, K , and solves for the 
frequencies at which the driving magnetic current required to maintain the electro¬ 
magnetic fields goes to zero. This is captured in the equation 


D(k 2 ) = kc 


f H K dV 
f e (H ■ H) dR 


introduced previously. As the frequency, and thus the wavenumber k , is changed 
over the range provided, Superfish looks for locations of zero crossings of D{k 2 ) and 
where the first derivative with respect to k 2 is -1. These locations correlate to cavity 
eigenfrequencies and the resonant frequencies are then provided to the user in an 
output hie. 


b. COMSOL Model 

COMSOL Multiphysics is a finite-element code designed specifically to 
couple multiple physics areas together into single simulations. It consists of a core 
set of solvers and physics libraries with the ability to add additional libraries as the 
user requires. We used COMSOL primarily as a validation tool for the model as 
it has the ability to import CAD models directly. After obtaining the CAD model 
from Niowave, it was imported and established in the COMSOL environment. For 
the modal analysis and held solving performed, we used the electromagnetic waves 
package (called TM Wave) in the RF module and the dedicated eigenfrequency solver 
within that package to find the hrst 10 lowest eigenfrequencies. The COMSOL held 
map shown in Figure 93 for the fundamental mode, confirmed the locations of the 
peak electric and magnetic helds predicted by Superhsh. 

c. Comparison with Experiment 

To measure the resonant modes of the cavity, the machined pieces of 
the cavity were assembled prior to welding. To ensure good conductive contact at 
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Figure 93. COMSOL field map of the simplified NPS SRF booster showing relative 
electric held strength and direction (arrows) and relative magnetic held strength (red, 
strongest —> blue, weakest). 


the joints between pieces, weights were placed on top of the cavity. Figure 94 shows 
the cavity as it was tested for resonant modes. An RF antenna was placed near the 
beam exit as the transmitting antenna and a receiver was located near the nose cone 
(this connection is not visible in the figure). A network analyzer was used to scan 
across frequencies looking for strong signals on the receiver lead. 

The network analyzer scans across a user determined window of fre¬ 
quencies and measures the signal loss between the transmitted and received signal. 
For frequencies where the cavity is not resonant, losses are high. When a resonance is 
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Figure 94. NPS gun/booster cavity configuration for eigenfrequency testing. On the 
left, is the cavity with the output end removed to see the nose cone. On the right, 
we see the cavity as tested. The blocks on top are there to provide weight to make 
good electrical contact with the sides of the cavity. 

found, the output signal from the network analyzer spikes, as can be seen in Figure 95. 
If one zooms in on one of the spikes in the wideband sweep, the classic resonance plot 
appears. As we see in Figure 95, prior to welding, the cavity’s fundamental frequency 
is 499.55 MHz. This is about 450 kHz lower than design, but the cavity will shrink 
slightly in the welding process and will contract when taken to cryogenic tempera¬ 
tures. Both processes will increase the resonant frequency. The data collected from 
this experiment and comparisons to predicted resonant modes is reported in Table 6. 

Reviewing the results compiled in Table 6, we see that there is excellent 
agreement between the two simulation codes and the actual measured cavity frequen¬ 
cies. The overestimation seen in the fundamental frequency can be traced to the fact 
that the simulations use the ideal dimensions provided to the manufacturers. The 
measurements are taken on an unwelded, warm cavity, and the fundamental frequency 
can be expected to increase after welding and when the cavity is cooled. Through the 
7th eigenfrequency, the agreement between cavity and simulations continued to be 
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Figure 95. Screen captures of the network analyzer data from the NPS cavity in the 
booster configuration prior to welding. On the left is the fundamental resonance at 
499.55 MHz. On the right is a wide band sweep looking for higher order modes. 


exceptional. Frequencies above 3 GHz were unable to be measured with the frequency 
analyzer on hand at the time of the experiment. 


Table 6. NPS gun/booster cavity eigenfrequencies and comparison with simulation 


Eigenfrequency 

Number 

COMSOL 

[MHz] 

Superfish 

[MHz] 

Measured 

[MHz] 

1 

501.18 

500.85 

499.55 

2 

1,091 

1,091 

1,089 

3 

1,597 

1,597 

1,596 

4 

2,107 

2,110 

2,103 

5 

2,271 

2,270 

2,272 

6 

2,405 

2,402 

2,402 

7 

2,629 

2,628 

2,623 

8 

3,034 

2,988 


9 

3,073 



10 

3,501 
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Bead Pull 


A bead pull experiment is a method for mapping the electric fields of a cavity. 
The basic process is to pass a small conducting bead through the interior of the 
operating cavity. As the conductor passes through the fields of the cavity, it expels 
the fields, slightly changing the cavity phase and resonant frequency. From this phase 
shift, one can calculate the electric held as a function of position along a line through 
the cavity. 

As we are most interested in the on-axis fields, a thin fishing line was fed 
through the center of the cavity using the setup shown in Figure 96. The line was 
connected to a small electric motor attached to a spool for adjusting the bead’s 
position in the cavity. The bead is a small piece of stainless steel hypodermic needle 
threaded onto the fishing line and held in place with knots. A weight is attached to 
the free end of the line to maintain constant tension on the line via a pulley system 
that also ensures the line passes through the cavity as close as possible to the electric 
axis. Using a ruler, adjustments are made in the bead position within the cavity and 
the received RF amplitude and phase shift were measured. 

Using the motor, the bead was pulled straight through the cavity from exit 
to nose cone and the phase shift recorded as a function of time. The start and end 
positions were noted over a total transit distance of 20 cm and recorded for later use. 
We determine the held profile from the frequency shift as a function of phase shift 
recorded at a specific bead position [43], [44], 


s_i 

fo 


$ 4 > 

2Q? 


(IV.19) 


where 5f is the frequency shift, f 0 is the fundamental frequency of the cavity, Q o is 
the cavity quality factor, and 5(j) is the phase shift observed. We can relate this to 
the holds at the bead position [44] through the equation 


\E\ 


tf 


U f 0 e 0 vri?3 ’ 


(IV.20) 
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Figure 96. LCDR Sean Niles stands next to the measurement system for the cavity 
bead pull. A transmitting antenna is inserted from the top of the cavity and a receiver 
below. The pulleys are connected via thin fishing line and attached on one end to 
an electric motor and a weight on the other to keep the line taught and to provide a 
reference point for linear position measurements. 


where \E\ is the magnitude of the electric held (only E z on axis for this cavity and 
frequency), U is the stored energy in the cavity (assumed constant as the perturbation 
due to the bead is assumed small), and R is the radius of the bead. Combining 
equations (IV.19) and (IV.20), we find 

\E\ 2 _ S(j) 1 

~W~ - 2Q e 0 vr R 3 

= 2Q7re 0 A 3 ^ 

\E\ oc y/ty . (IV.21) 
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Taking the data displayed in Figure 97, the phase offset is subtracted to obtain a true 
zero. We also use the magnitude of the phase shift to avoid imaginary values when 
applying equation (IV.21). 


Bead Pull Phase Shift 



Figure 97. Phase shift as a function of position of the bead within the NPS booster 
cavity. The pull required 12.248 s to traverse 20 cm of the cavity and beam tubes. 


To make a meaningful comparison with the field profile from Superfish, a plot 
of E z is obtained along the cavity axis (in the booster configuration) and normalized 
to the peak field value. The peaks of the Superfish plot and the normalize peak of 
the phase data plot are aligned and the relative positions of the phase plot adjusted 
(time reversed since we measured from exit to nose cone) in order to compare the two 
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data sets. The resulting plots are plotted in Figure 98. From this plot, the bead pull 
profile shows good agreement with the Superfish profile. 


Bead Pull Field Profile 



- bead pull 

- Superfish 


z [cm] 


Figure 98. The electric field profile generated from a bead pull of the NPS booster 
cavity. The normalized bead pull plot is compared to a normalized plot generated in 
Superfish. 
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3. SC Solenoid 

A superconducting solenoid is installed just after the cavity to provide beam 
focusing for electron beam transport and for emittance compensation. Due to the 
solenoid being located near the liquid helium dewar, it is chosen to be supercon¬ 
ducting to allow for high amperage currents (and thus, high magnetic fields) without 
generating a measurable heat load to the cavity. The solenoid is an aluminum spool 
wound with niobium-titanium wires embedded in a copper matrix and backed with 
an iron backer to provide a return for the magnetic held line circuit. The solenoid 
has 2,236 turns and was designed for a maximum 10 A current resulting in a peak 
magnetic held of 476 mT (more than twice the anticipated necessary held for beam 
focusing). 

a. Superfish Model 

The hrst models were developed using Superhsh to predict the held 
prohle and maximum held strengths in the iron back and along the beam axis. The 
held maps shown in Figure 99 indicate that within the normal operating parameters 
of the solenoid (I so [ < 5 A), the iron backer should not saturate. Above about 5.5 
A, saturation begins to be a minor concern. The held prohle shape is as expected 
and desired. The held extent is sufficient to affect the SRF cavity when it transitions 
to superconducting, so we must remain vigilant to ensure no residual magnetic held 
remains on the solenoid to prevent cavity held degradation. Should the iron core 
become magnetized, this can be corrected by degaussing the solenoid prior to the next 
cavity transition. The hux trapping process will be discussed later in this section. 
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Magnetic field date from file SOLENOID.POT 
Problem title line 1: Superconducting solenoid tor Niowave SET gun - v2 



Figure 99. The superconducting solenoid field map generated by Superfish. The lower 
plot is a 2D map showing the magnetic field lines of the solenoid. The inset at the top 
is the magnetic field strength along a line through the magnetic axis of the solenoid. 
The peak magnetic field along this line occurs at the center of the solenoid and B p k 
= 4,760 Gauss (476 mT). 
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b. 


COMSOL Model 


Performing the same simulation in COMSOL provided us with confi¬ 
dence about the peak fields and saturation onset currents. Figure 100 shows the held 
map and on-axis magnetic held prohle obtained from COMSOL. The results from 
COMSOL agree extremely well with the results obtained from Superhsh. 
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Figure 100. The superconducting solenoid held map generated by COMSOL. The 
lower plot is a 2D map showing that the maximum magnetic hux density (contours) 
is 2.66 T, occurring at the inner corner of the iron backer and the relative strength of 
the magnetic held (arrows). The inset at the top is the magnetic held strength along 
a line through the magnetic axis of the solenoid. The peak magnetic held along this 
line occurs at the center of the solenoid and B p *. = 476 rnT. 
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c. 


Comparison with Experiment 


To verify the design, the solenoid held was measuring using a Hall 
probe in a tabletop test. This type probe has a Hall effect sensor at the tip of a long 
wand that measures magnetic held strength. The experimental measurement set up 
is shown in Figure 101. The hall probe is mounted in a nonmagnetic plastic block 
that is machined to ht into a U-channel aluminum bar. The U-channel and meter 
stick are passed through the solenoid and the block adjusted such that the probe will 
pass as near to the solenoid axis as possible. As the probe is slowly moved through 
the solenoid, magnetometer readings are taken. Measurements were made entering 
and exiting the solenoid from both sides to ensure any bias was averaged out. 

Since the solenoid measurements were made at room temperature, we 
necessarily performed them at low current. There is no reason to believe that the 
held profile will change when the temperature is decreased and the current is in¬ 
creased where the solenoid is at superconducting temperatures. Figure 102 shows the 
experimental data compared with the appropriate plot generated from Superfish for 
three different applied currents (20 mA, 40 mA, and 60 mA). All plots show excellent 
agreement between simulation and experiment. 

The solenoid was then placed in a cryogen bath to test its supercon¬ 
ducting properties. As the solenoid was cooled, its internal resistance was measured 
using a standard 4-wire measurement. At room temperature, the solenoid had a re¬ 
sistance of 165 Q. Niobium-titanium has a transition temperature of about 15 K and 
we expect the resistance just prior to transition to be about 1.6 D. As expected, the 
resistance dropped to 0.8 H as the sensing wire resistance steadily increased to 1 kh2 
indicating the solenoid was superconducting. The solenoid was then slowly ramped 
to its design current of 10 A and operated there for at least a minute without quench. 
During the ramp up of current, the power supply voltage remained linear with the 
current increase indicating no significant resistive heating in the power leads. To test 
the solenoid at quench, we applied a 10 A current to the solenoid and then lifted it 
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Figure 101. The experiment setup for measuring the magnetic field strength through 
the center of the superconducting solenoid. 


from the liquid helium bath. As the temperature increased to transition, we saw no 
increase in the sensing line voltage until a large spike was encountered, indicating the 
solenoid was no longer superconducting. The solenoid was immersed back into the 
liquid helium bath, but due to resistive heating did not return to the superconducting 
state until current was removed to allow for the solenoid to cool. The solenoid was 
removed briefly, Figure 103 for inspection before returning to the bath to measure its 
inductance. 
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Figure 102. Experimentally determined axial magnetic field maps of the supercon¬ 
ducting solenoid compared to Superfish models for various applied currents. Magnetic 
fields in this plot are measured in gauss. 


To measure the inductance of the superconducting solenoid, we recall 


V a 



(IV.22) 


where V s is the voltage across the solenoid and C is the solenoid inductance. To 
obtain the current as a function of time, the power supply voltage is ramped in order 
to maintain a constant V s . Differentiating and rearranging Ohm’s law provides us 
with a method to calculate the time rate of change of the current as the applied 
voltage changes, 


d 1 _ 1 dV 
dt 7 Z dt 


(IV.23) 


Using the power supply voltage as a function of time recorded in Figure 104, we can 
calculate the solenoid inductance. The line resistance is 7 Z = 100 mO and AV/At = 
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Figure 103. Inspecting the superconducting solenoid after removal from the liquid 
helium bath. 


600 mV / 6.6 s, giving a rate of change in the current of 0.91 A/s. Putting this value 
into our inductance equation, we find the solenoid has an inductance of about 0.55 
H. 
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Figure 104. Power supply voltage as a function of time used to calculate the induc¬ 
tance of the superconducting solenoid. 
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d. Comparison with Simulation 

For a solenoid to focus an electron beam, it must also cause the beam 
to rotate. This rotation is generated in the fringe fields of the beam, causing the 
beam to rotate as it enters the fields and then stop as it exits the fringe fields. 
The imparted transverse motion causes the beam to focus at some distant point 
before expanding again. Using GPT, we can model this motion and compare with an 
analytical calculation as a benchmark of both the simulations and our experimental 
measurements. To visualize the rotation of the beam, a uniform beam is broken 
up into four beamlets, each of which makes up 1/4 of the beam as seen in Figure 
105. As the beam propagates, we observe that the fringe fields begin to rotate the 
beam (top right) and in the final picture, we see the output beam. Particle mixing 
has occurred during propagation due to internal beam Coulomb forces, but the lines 
between beamlets are still easily identifiable. From the simulation, we can estimate 
that the 1.2 MeV beam has rotated approximately 45° as it passed through the 
solenoid with B p = 117 mT. 

We can perform a hand calculation to compare against what we observe 
in the simulation. The rotation angle experienced by a beam as it passes through 
a solenoid is related to the magnetic field strength, held extent, and beam rigidity 
through the relation [45] 


29 


BL 

W)’ 


(IV.24) 


where 6 is the rotation angle, B is the magnetic held strength, L is the magnet’s 
effective length, and (Bp) is the beam rigidity. The rigidity of an electron beam can 
be found by applying [46] 


Bp = 3.3356 p [T m] , 


(IV.25) 


where p is the momentum and is given in GeV/c. To obtain the momentum, we must 
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first find 7 from the total energy equation, 


777 T 0 C 2 = K E + moc 2 
KE 

7 = 1 +-j , (IV.26) 

rriQC z 

where uiq is the rest mass of the electron. We know the kinetic energy of the beam 
is 1.2 MeV and can obtain the integrated held from the Superfish held maps we have 
generated (0.469 T m) at the stated peak magnetic held value. From these values, 
we calculate 7 = 3.35 and can calculate p = 7 (3mc = 1.63 MeV/c. Using these values 
and combining equations (IV.26) and (IV.25) and substituting into equation (IV.24), 
we have an expected rotation angle through this solenoid of 


6 


BL 

2 (3.3356 p) 


43 c 


which is in excellent agreement with the simulation. 


(IV.27) 


168 



1.2 MV Cavity, 10k Particles, Cross-sections 
position=0.006 


1.2 MV Cavity, 10k Particles, Cross-sections 
position=0.18365 




o.oio 


0.005 


-0.005 


1.2 MV Cavity, 10k Particles, Cross-sections 
position=3.66625 


0.03 

0.02 

o.oi 

>* 0.00 

- 0.01 

- 0.02 

-0.03 


, 

’ la ffflii L 

Ys 



Figure 105. Cross sections of a 1.2 MeV beam from the NPS gun passing through the 
superconducting solenoid. Top left is the beam just after launch from the cathode. 
Top right is the beam at the center of the superconducting solenoid. Bottom is the 
beam after propagating to an output screen approximately 3.5 m from the solenoid. 
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e. Thermal Transport 

Recall from the previous chapter that the superconducting solenoid is 
located outside the liquid helium dewar and thermally connected to the dewar via two 
solid, copper bus bars. It was hoped that this thermal connection would be sufficient 
to allow the solenoid to cool below T c . However, during initial cold tests of the cavity, 
it was discovered that this cooling connection was insufficient to cool the solenoid 
fully. In normal operating mode, with no RF applied to the cavity, and no current 
applied to the solenoid, the coolest temperature, as measured on the outside of the 
solenoid case, was 11.3 K. 

Because the solenoid is not superconducting, any current introduced 
into the solenoid will travel through the copper portion of the wire and cause resistive 
heating. This generated heat is tied to the liquid helium dewar and will cause an 
additional heat load to the helium bath, decreasing the operating time available on a 
single fill. Without opening the cryomodule, there is no way to increase the thermal 
connection of the solenoid to the helium dewar. 

It is possible that there are additional external heat loads that are 
preventing the solenoid from cooling sufficiently and becoming superconducting. Our 
initial thought was that the coupler, which passes through the solenoid inside the 
beam tube, could be causing radiative heating since it is not cooled in any way. To 
test this hypothesis, we wrapped the vacuum flange where the coupler short connects 
to the beam tube to the off-gassing line of the cryomodule liquid nitrogen shield and 
insulated it. After an hour of applied cooling, no significant change in the solenoid 
temperature was seen. It is possible that the conductive heating from the beam tube 
could be the additional heat source we seek to mitigate, but this has not been tested 
to date. 

To avoid introducing too much heat into the cryomodule with the 
solenoid, the current operating procedure is to run the solenoid in pulsed mode. 
The present gun status is such that the beam is of low enough energy that the fo- 
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cusing fields generated with the reduced solenoid current are sufficient for present 
experiments. Solving this heat transfer problem is essential in the next generation of 
the NPS SRF gun. 

f. Flux Trapping 

Most physics students have at least a passing familiarity with the Meiss¬ 
ner effect in passing—usually in the context of the simple demonstration of a magnet 
placed on a superconducting material that is taken through the transition tempera¬ 
ture. As the superconductor’s temperature drops below T c , the magnetic held from 
the magnet is expelled and the magnet levitates above the superconductor surface. 
The Meissner effect in a superconducting cavity behaves in a similar way. With the 
simple physics demonstration, both the magnet and superconductor are usually small 
and the effects of trapped flux are typically not observed. In a superconducting RF 
cavity, the raw amount of superconducting material makes the likelihood of impurities 
being present in the bulk material almost definite. These impurities can cause signifi¬ 
cant issues if the superconductor transitions in the presence of a magnetic held. Flux 
trapping is a condition where impurities or lattice defects in the superconductor act 
as trapping sites, resulting in only partial hux expulsion [12]. The pinned magnetic 
hux lines result in increased residual resistance increasing the cavity losses [47]. These 
increased losses result in a decreased cavity Q 0 and associated accelerating helds. 

Since the NPS cavity design has a solenoid capable of developing a 
nontrivial magnetic held near a superconducting cavity, it is essential that we under¬ 
stand how this affects performance. As the solenoid is used, particularly at higher 
held levels, it is possible to induce a residual magnetic held in the solenoid when the 
current is removed. These held lines can pass through the walls of the cavity during 
transition, and, if not fully expelled, result in degraded cavity performance. There is 
no issue if the cavity is maintained superconducting as the magnetic helds from the 
solenoid will not be able to penetrate the cavity surfaces. 
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To test this effect and observe the impact on this particular cavity, we 
purposely attempted to induce flux trapping during the first cold tests of the cav¬ 
ity. In the first step, we degaussed the solenoid prior to cooling the cavity below T c 
to ensure no residual fields from previous experiments remained. We then took the 
cavity superconducting and measured a cavity Q o value as a baseline. The solenoid 
was operated briefly at moderate current (/ = 4.2 A) to magnetize it. After the cav¬ 
ity quenched, we refilled the helium dewar and brought the cavity back to the same 
operating state to measure another Q 0 value. Once back to the previous operating 
temperature, the Q 0 measurement indicates a moderate decrease in cavity perfor¬ 
mance. We then allowed the cavity to quench again and degaussed the solenoid to 
remove any residual magnetic fields. After cooling the cavity again, we were success¬ 
ful in returning to the original Q 0 value, indicating that the cavity can be affected 
by residual solenoid fields, but that the effect can be corrected by degaussing the 
solenoid when the cavity is normal conducting. 

4. Cold Cavity Processing 

When an accelerating cavity is first brought into service, it must be condi¬ 
tioned prior to full use. Conditioning typically involves putting the cavity into the 
expected operating environment (under vacuum, cryogenic temperatures, etc.) and 
placing it under an RF load. As RF is put into the cavity, any held emission sites will 
begin to emit and hopefully condition themselves away. Multipacting barriers may 
also be encountered as RF power levels change. Typically, once a cavity has been 
conditioned to a certain RF power level, unless something in the system changes (such 
as venting the vacuum), the cavity can be operated safely and easily below the max¬ 
imum conditioned level. During conditioning, both held emission and multipacting 
are generating accelerated electrons that run into the walls of the cavity and beam 
tube and generate X-rays as they are stopped by the metal. By monitoring radiation 
levels, we can estimate the level of conditioning that takes place. 
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The NPS gun/booster cavity is designed to be a batch fill cavity: the cavity 
will not, as many other SRF cavities are, be cooled once to liquid helium temperatures 
and maintained there for long periods of time. As the cavity temperature is taken 
from room temperature down to 4 K, it must transition through a temperature region 
where the possibility of “Q-disease” becomes increased. Q-disease is the formation of 
hydride precipitates on the surface of the cavity [12]. These hydrides have much higher 
electrical resistance and contribute to localized heating of the cavity, possibly creating 
localized regions where the temperature of the cavity exceeds T c . The majority of 
the dissolved hydrogen forming the hydride precipitates comes from the acids used to 
etch the cavity surfaces when the bath temperature gets too high. 

The typical cool down process used for the NPS cavity is to flow liquid nitrogen 
through the nitrogen shield as well as into the helium dewar. Once temperatures in 
the cavity approach 150 K, the onset of greatest hydrogen diffusion rate, we allow the 
liquid nitrogen to boil off and switch to putting liquid helium into the dewar. The 
goal is to pass through the temperature region of 66 K < T < 150 K as quickly as 
possible [12]. The result is a decreased risk of Q-disease at the expense of a significant 
boil off of helium. After five cycles between room temperature and 4 K, no decrease 
in cavity Q 0 has been observed to date. 

Some of the experimental values that are typically desired from cold cavity 
testing are estimates of the cavity’s Qo value and the resulting cavity gradient. Plots 
of Q o values are also useful in diagnosing any issues that may be encountered during 
cold cavity testing, such as field emission, Q-disease, and multipacting. 
a. Helium Processing 

As part of the processing of the NPS cavity, the first ever helium pro¬ 
cessing in a commercial facility was performed. In helium processing, a small amount 
of helium gas is introduced into the cavity vacuum to serve as an ionization source 
[48], [49]. As electrons are emitted from field emission sites, the electrons ionize the 
helium gas causing the helium ions to back bombard the emission sites. The helium 
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ions have greater momentum and the resulting back bombardment of the emission 
sites increases the local temperature more than the ohmic heating and processes them 
down faster. Experience at Cornell [48] indicates they have observed factor of 2 in¬ 
creases in accelerating gradient using helium processing on their cavities. Estimates 
of accelerating gradient improvement in the NPS gun were just shy of this mark, 
improving the estimated accelerated gradient from 400 kV to 750 kV. The Qq plots 
calculated by Niowave compares this series of processing runs to those from the initial 
cold testing, as shown in Figure 106. 



integrated voltage (MV) 


Figure 106. Two Q 0 vs. E plots for the NPS booster cavity are plotted. The black 
dots correspond to measurements made during the initial cold testing in September 
2009. The red dots were taken after considerable RF and helium processing in October 
2009. The September 2009 measurements were limited by X-ray generation. The 
October 2009 measurements were field emission limited. 
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5. Full Beam Test 

After the cavity was vented and the cathode stalk installed as discussed in 
the previous chapter, all systems, including the drive laser, were ready for first beam 
testing. First beam testing was conducted at Niowave facilities the week of 07-11 
June 2010 and established the NPS 500 MHz SRF gun as the first operational SRF 
gun in the United States. The gun transported an estimated 30 pC from a niobium 
cathode. Continued testing to fully characterize the gun has been ongoing since that 
date, but that data will not be reported here. Interested readers are encouraged to 
consult the literature for more recent published results. 

The insertion of the cathode constitutes a major change to the cavity and 
could possibly affect how the system operated and behaved. Prior to conducting 
beam tests, we characterized the system and conducted processing runs to process 
any new held emission sites introduced as part of the venting process or present on 
the cathode assembly. While there were a few held emitters and multipacting barriers 
that required processing, the largest impact appeared to be the loss of cavity gradient 
due to coupling to the cathode stalk. With the cathode fully retracted (about 14 mm 
behind the nose cone), Q o measured approximately 5 x 10 8 . Adjustments to the 
cathode position placed it at an estimated retraction of about 9 mm behind the nose 
cone with a Qo of 6 x 10'. However, there is no way to currently describe the cathode 
position with confidence. 

We know that changing the cathode position alters the helds inside the cavity 
and has an effect on the fundamental frequency of the cavity as well. Using the 
cathode itself as the probe, we can simulate the frequency shift as a function of cathode 
position over the range of desired positions, shown in Figure 107. Assuming we can 
verify this plot between beam performance characteristics and measured fundamental 
frequency against the simulations, we will have a useful diagnostic for determining 
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the position of the cathode. Full experimental verification has not yet been performed, 
however the first two data points indicate the slope of the curve appears to match 
the cavity frequency shift expected. 


Frequency Impact of Cathode Position 



Figure 107. Results of Superfish simulations of the NPS SRF gun showing resonant 
frequency vs. cathode retraction. As the cathode position relative to the nose cone 
is adjusted, the fundamental frequency, / 0 , of the gun changes. It is possible that 
this effect could be used to measure cathode position once verified experimentally. 
The first two measured resonant frequencies as the cathode position is changed are 
plotted in red. 


After RF processing, the gun was not able to return to field levels obtained 
in the booster configuration. It is believed that this is due primarily to the coupling 
of power out of the cavity by the cathode stalk. As the stalk was inserted nearer 
to the nose cone surface, the Q 0 of the cavity was seen to decrease, indicating a 
loss of stored power in the cavity. The cathode was returned to the 9 mm retracted 
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position and an estimated accelerating gradient of 200 kV reverified. At this point, 
it was decided that, although the final gap voltage was not as high as desired, it was 
enough for performing initial beam tests. First beam with the NPS 500 MHz SRF 
gun was achieved on the first attempt at 1110 EDT on 09 June 2010. Once beam was 
established, the team proceeded to commission the diagnostic beam line, verifying 
operation of all correctors and transporting electron beam to every instrumentation 
station in the beam line over the course of only 12 hours. 

a. Instrumentation Stations 

The instrumentation stations for the NPS gun beam line are the pri¬ 
mary method for collecting data about the electron beam. There are five instrumenta¬ 
tion stations along the main beam line and an additional two along the spectrometer 
leg. With the exception of the first station, the stations have identical construction. 
Figure 108 shows the pertinent portions of the stations. Only the first instrumentation 
station has the pinholes for performing emittance estimates. 
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Cerium YAG 



Figure 108. Each instrumentation station consists of an electrically isolated rod at¬ 
tached to a stepper motor driven linear actuator for insertion. Each station has 
a Cerium YAG crystal mounted in front of a polished aluminum mirror for beam 
viewing. The aluminum plate below the screen has pinholes on the first station for 
omittance measurements; for the others, this is a blank plate. The blank plate sections 
can be used as Faraday cups for measuring bunch charge. 


Each station can be independently positioned to one of two instruments 
or retracted so as to present a clear aperture. The upper-most consists of a Cerium 
doped YAG screen which scintillates when struck by the electron beam. Cameras 
view the beam via top-mounted, down-looking video cameras via a polished aluminum 
mirror behind the Cerium YAG screen. Pin holes are installed in the first instrument 
station for performing omittance measurements. The bottom position of each station 
is intended for use a Faraday cup. Each station is electrically isolated, so any charge 
impacting the instrument carriage can be captured for bunch charge measurements. 
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Figure 109 shows some representative screens from the first beam test. For relative 
size measurements, the white line in each picture is the side view of the screen and 
each screen is 25.4 mm in length. 



Figure 109. Beam cross sections captured from stations 2 (bottom left), 3 (top left), 
and 4 (right) during the first beam tests of the NPS 500 MHz SRF gun. Brightness 
is proportional to beam intensity. The bright line at the top of each picture is the 
side view of the Cerium YAG crystal. 
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b. Experimental Verification of Electron Beam Energy 

To verify estimates of the cavity gradient, we prepared simulations to 
calculate the beam energy using the “window-frame” beam correctors, window-frame 
correctors are solenoid coils wrapped around a square iron ring, as seen in Figure 
110. The solenoids are operated in pairs; each pair provides a deflecting field in one 
plane. While there are multiple window-frame correctors installed along the beam 
line, the intent is to use the first one to deflect the beam in the vertical plane (to 
avoid impinging on the laser mirror). Based on the same principle of the solenoid 
calculation, a beam of specific energy has a certain rigidity or resistance to deflection. 
For the same magnetic field a beam with a higher energy will deflect less than a 
similar beam with lower energy. 



Figure 110. On the right is the window-frame corrector as installed on the diagnostic 
beam line. On the left is the field map through the axis of the corrector measured in 
the same manner as the solenoid field. A Gaussian fit to the field profile is plotted as 
well, for comparison. 


Based on this concept, a series of simulated beams of multiple ener¬ 
gies, using the NPS gun as the source, were propagated down a reduced beam line 
(terminating at the second diagnostic station). To model the window-frame corrector 
in the absence of actual held maps, a Gaussian approximation to the held was used. 
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This approximation was generated numerically in MATLAB, exported to a text file 
and imported into GPT. Actual field maps were inaccessible for these simulations as 
GPT and Superfish both require solenoidal fields to have their magnetic axis aligned 
to the beam propagation direction. For simplification and expedience, an approx¬ 
imate method was used. Hand calculations using an impulse approximation were 
performed to verify the simulation results. Figure 111 shows the results of the sim¬ 
ulation runs, providing a simple process for determining beam energy based on the 
measured centroid deflection on the YAG screen and the applied corrector current. 


WF Corrector Kick (SCI - IP2) 
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Figure 111. A look-up chart for determining beam energy based on corrector current 
and centroid deflection. The ragged values near 0.035 m deflection is due to beam 
clipping on the beam tube wall during propagation. 


The first estimate of beam energy was obtained after commissioning 
the beam line to the first instrumentation station. With a corrector current of ~0.66 
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A, we observed a deflection of ~7.5 mm. Correcting this deflection for the additional 
propagation distance to the second instrumentation station and consulting the chart 
in Figure 111, we estimated the beam energy at 400 keV. The impulse calculation 
confirmed this result. However, RF calculations by Professor Brian Rusnak estimated 
the maximum beam energy as no higher than ~200 keV. Once the remainder of the 
beam line was commissioned, the dipole magnet was used to perform a second energy 
measurement, finding the beam energy to be approximately 200 keV. 

Attempts to duplicate the measurement on the following day after per¬ 
forming further RF conditioning were complicated by difficulty transporting the beam 
to the termination of the spectrometer leg. Using a sweep of the beam through the 
first instrumentation station after the dipole, along the spectrometer leg, indicated a 
beam energy of approximately 250 keV. 

Possible sources of error in the window-frame energy measurement 
could be misalignment of the beam prior to entering the corrector field. If the beam is 
already pointed off the beam line axis, any deflection measurement would be skewed 
by this error. Additionally, it was determined that the light-tight box around the 
UV drive laser insertion table could be magnetized, which could affect the accuracy 
of measurements performed with the window-frame correctors in this area of the 
beam line. However, with correlation between the dipole beam energy measurements 
and the RF calculations, we feel confident that the NPS SRF gun has developed a 
minimum 200 keV beam of electrons. 

c. Experimental Verification of Bunch Charge 

Calculation of bunch charge in the SRF beam line is performed in a 
similar fashion to that described in the DC gun beam line. In this case, there is 
only a single metal plate upon which the beam impinges. To verify that the beam is 
centered on the Faraday cup section, the electron beam is first steered and focused 
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onto the Cerium YAG screen. Then the Faraday cup is moved into the proper mea¬ 
surement position. Traces on the oscilloscope were almost rectangular and time broad¬ 
ened (about 10 ns), as seen in Figure 112. 



Figure 112. Faraday cup signal plots for three different cases. The blue is the signal 
from 1st beam. In red is the maximum signal strength from the same system set up 
after adjusting the laser position on the cathode. The green is a general plot from 
later in the experimental runs. 


From a photodiode calibration, our Boeing collaborators were able to 
tell us the laser pulse energy for these three plots was about 180 /rJ per pulse. Us¬ 
ing this value with the same process as shown in equation (IV. 18), and modeling 
the pulse as a rectangular pulse, we find the bunch charges and quantum efficiency 
estimates shown in Table 7. Additionally, after adjustment of the laser and further 
transport optimization, the largest voltage pulse obtained during this series of runs 
is shown as well. The QE measurements are about an order of magnitude less than 
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anticipated when compared to values in the literature. However, the majority of the 
literature values are reported after having had some kind of processing of the surface 
beyond the BCP that was performed on ours. Future work includes performing laser 
cleaning/preparation of the cathode surface and these measurements provide us with 
a baseline value with which we can measure the improvement seen. It should also 
be noted that these values are preliminary in that the transport and launch phase 
are only roughly optimized. Any adjustment resulting in increased charge transport 
will show a corresponding increase in QE. It is possible that further adjustments will 
improve these measurement values prior to any cathode surface alterations. 


Table 7. NPS gun bunch charge and QE measurements from the Erst experimental 
run. 


Run Name 

V 

mV] 

Epulse [/^J] 

q [pC] 

QE 

1st Beam 

100 

180 

20 

5.2 x 10~ 7 

Highest Signal 

370 

180 

74 

1.9 x 10~ 6 

General 

110 

180 

22 

5.7 x 10“ 7 

“Best” 

550 

300 

110 

1.7 x 10~ 6 


d. Bunch Charge vs. Launch Phase 

Using the Faraday cup in the second instrumentation station as our 
detector, we can plot an acceptance plot of the gun for the 9 mm retracted cathode 
position. We adjust the launch phase by changing the phase delay between the RF 
signal and the laser pulse by adjusting a manual RF phase shifter. Using this process, 
we performed this measurement twice, once for an estimated 250 keV beam and once 
for a 300 keV beam with the plots shown in Figure 113. We see that the cavity has 
a relative phase acceptance of about 150° for the 250 keV beam and 165° for the 300 
keV beam as measured between roll-off points. Comparing these acceptance window 
widths with simulation results reported earlier (of about 140°) we hnd them to be in 
good agreement. 
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Figure 113. NPS gun relative phase acceptance for beam energies of 250 keV and 300 
keV. Bunch charge was measured at the second instrumentation station. 
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V. 


FEL BACKGROUND 


In the next chapter, a scheme for describing the optical fields of a free electron 
laser (FEL) will be described. This chapter will provide a moderate background for 
the reader unfamiliar with FELs, but will not perform detailed derivations. For the 
interested reader, there are many accessible texts and papers [50], [51], [52] available 
for further study. This background chapter will start with the basic history of the 
FEL leading to a discussion of the major components of most FEL systems and the 
basic physics of the FEL. 

A. HISTORY AND COMPONENTS 

The FEL, as its name implies, is a laser. However, it is unlike other lasers 
in that the lasing medium does not have specific energy levels in which a population 
inversion must be established. In a typical laser, a pumping mechanism, whether a 
chemical reaction, flash lamps, or other technique is used to elevate electrons to an 
excited state. Once in this excited state, photons are collected and amplified through 
stimulated emission as the elevated energy level is depopulated. 

1. History 

The technology of the FEL has its origins in the microwave tubes of the 1930s, 
predating the solid state laser. With microwave tubes, scientists discovered they could 
generate coherent radiation from the oscillations of electrons. This coherent radiation, 
typically in the RF wavelengths was collected and passed into waveguides for many 
purposes. Depending on the size of the resonator cavities used, the microwave tube 
could be used to generate wavelengths of that fundamental size. This limited these 
tubes to wavelengths associated with resonant cavities that could be easily manufac¬ 
tured, typically on the order of centimeters. Although it was recognized that coherent 
radiation could be obtained from a beam of free electrons, the connection between 
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the microwave tube and the laser, first demonstrated in 1960, would take another 10 
years before the concepts were combined. 

The first laser, demonstrated at Hughes Research Labs [53] in 1960, was a 
solid state laser. The lasing medium was a solid laser crystal, pumped with a flash 
lamp. An optical resonator was put around the crystal to provide feedback during 
amplification. Gas lasers, chemical lasers, solid state lasers, and almost all other lasers 
operate in a similar fashion. Each of these lasers is dependent on the lasing medium 
properties, usually requiring shut down or reduced duty cycle to avoid over-heating 
the lasing medium or replenishment of the constituent chemicals. Additionally, every 
laser’s operating wavelength is fixed to the lasing medium. If a different wavelength 
is desired, new constituents must be found. 

As initially proposed [54] and demonstrated [55], an electron beam passing 
through a periodic magnetic held can experience stimulated emission, resulting in a 
laser that has no fixed lasing medium. The fact that the laser wavelength is established 
through an adjustable system, the electron accelerator, the FEL can conceivably be 
designed to operate at any desired wavelength. In fact, today FELs have operated at 
wavelengths from cm (microwaves) to 1.5 A (hard X-ray) and nearly all wavelengths 
in between [56]. There have been many advances in technology since the first FEL 
lased at Stanford University in 1976. 

Over the years of development, the FEL output power level has not received 
as much attention as other laser types. There probably is not a single primary reason 
other than the excitement of being able to design a laser for any particular part of 
the spectrum kept the investment money for FELs focused primarily on driving to 
shorter wavelengths. Since most conventional lasers share the same optical principles, 
once a laser was made to work at a particular wavelength, the main investment focus 
was in improving efficiency and output power. The FEL is finally receiving some 
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of the same attention, due primarily to military interest, driving FEL output power 
from watts, to kilowatts [57], to tens of kilowatts [58], with designs for a 100 kW laser 
underway [59]. 

2. FEL Components 

All FELs have some common components that are important to the overall 
system performance. Some basic description of the important, major portions of a 
notional FEL will be discussed with the next section describing more of the physics 
and mathematical details. Most FELs can be described by the system layout similar 
to that shown in Figure 114 from Jefferson Laboratory, Newport News, VA. 



UV higgler 


Figure 114. A system layout diagram of the JLAB FEL. While specific to one machine, 
this layout can be extrapolated to many other FEL systems. From [60]. 


a. Injector 

All FELs use free electrons as their lasing medium. Typically the source 
of electrons is located in what is called the “gun.” As explained earlier, the two 
primary types of guns are radiofrequency (RF) and constant voltage (DC). Both of 
which can have different types of cathodes and can be normal, superconducting, or 
even a combination. There are trades-offs that must be carefully considered for both, 
but both types of guns are in use. Thermionic and photoelectric cathodes are the 
two primary electron sources in the guns and, again, there are trades necessary in 
determining the best type for the application intended. 
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The injector’s primary purpose is to bring the energy of the bunch 
launched from the gun high enough to avoid slippage with respect to the RF fields as 
the bunch enters the primary linear accelerator (LINAC) sections. The injector may 
also have a dispersive element to adjust the electron bunch phase space for proper 
matching with the main LINAC. 

b. Linear Accelerator 

In the LINAC is where the electrons are given their primary boost to 
high energy. The LINAC can develop average kinetic energies from a few MeV to 
a many GeV depending on the wavelength desired in the end application. The only 
difference between machines is the amount of RF power put into the cavities and 
the operating frequency. The operating frequency will determine, to a certain extent, 
the maximum repetition rate at which the laser can operate—higher frequencies have 
more RF “buckets” for electron bunches. Currently, the standard rule of thumb is 
that the higher the repetition rate, the lower the bunch charge to prevent developing 
instabilities in the RF fields (wake fields) that would tend to break up the beam as 
more bunches pass through the LINAC. Multiple accelerating sections, or multiple 
passes through the same section, can be used to achieve the desired final beam energy. 

c. Undulator 

The undulator is the section of the FEL where the lasing interaction 
takes place. The undulator is a periodic magnetic field, which through the Lorentz 
force, causes the electrons to “wiggle.” This section of the FEL is also known as the 
“interaction region” and is directly analogous to the lasing medium of a standard 
laser. A diagram of a linear undulator is shown in Figure 115. 

In all FELs, the undulator causes the electrons to radiate in the for¬ 
ward direction as they accelerate side-to-side. In an oscillator FEL, the spontaneous 
emission is collected at a mirror and reflected back to a second mirror which sends 
the light back into the undulator. The next electron bunch arrives at the same time 


190 




Figure 115. A linear undulator with alternating magnetic fields cause the electrons 
to oscillate in the opposite transverse plane. From [61]. 

as the light pulse and the electrons now oscillate in the presence of an electric field 
stimulating them to radiate at the same wavelength. Over many passes, this laser 
light builds up to significant levels. 

SASE (Self-Amplified Spontaneous Emission) and amplifier FELs are 
single pass designs. In an amplifier FEL a seed laser of the desired output wavelength 
is chosen and sent through the undulator. The electrons see the electric field of the 
seed laser and amplify it through stimulated emission. In SASE, the electrons feed 
forward their spontaneous emission. As the spontaneous emission from the back of 
the electron bunch overtakes the front of the bunch, the front of the bunch amplifies 
this signal in the same manner described for the other types of FEL. 

d. Beam Dump 

The beam dump is essentially the completion of the electrical circuit 
that was started at the cathode. As the electron beam passed through the undula¬ 
tor and gave up energy to the laser, a significant energy spread is induced on the 
beam. Unfortunately, this means that the beam is unable to used again for the FEL 
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interaction. In some systems, the beam is dumped into a metal block at the full 
remaining beam energy. This can result in extremely large amounts of radiation. An¬ 
other option, as shown in Figure 114 and demonstrated in numerous FEL facilities, 
is to recirculate the electron beam through the LINAC 180° out of phase from the 
accelerating RF. Essentially, this recirculation scheme couples most of the left over 
energy of the electron beam back into the RF fields in the LINAC. Not only does 
this reduce the overall RF power required for the FEL, increasing overall efficiency, 
but it decreases the amount of kinetic energy dumped into the beam dump, which 
decreases the radiation generated. 


B. BASIC FEL PHYSICS 

In this section, we develop and examine the fundamental equations that govern 
FEL physics. We will start with the relativistic Lorentz force equations and apply 
them to an electron moving through the periodic magnetic held of the undulator. 
We will then discuss how an FEL is tunable by examining the resonance condition. 
Finally we will discuss the importance of peak current to the operation of an FEL. 

1. Lorentz Equations 

We start with the relativistic Lorentz equations, recognizing that for an FEL, 
we intend to use highly relativistic electrons, 

F=-e[I! + c(?xB)]=^ = |( 7 ii^) , (V.l) 


where 7 is defined in the usual manner, 


1 

' 


(V.2) 


2. Electron Motion in a Periodic Magnetic Field 

Let us first consider an electron in a periodic magnetic held without an electric 
held. We also assume a helical undulator for simplicity in calculation; the same 
methodology can be used to solve for a linear undulator. We define our undulator 
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field as 


B = B (cos (k 0 z ), sin (k 0 z ), 0) , 


(V.3) 


where k 0 = 27t/A 0 and A 0 is the wavelength of the undulator. Substituting equation 
(V.3) into equation (V.l), we find the equation for the acceleration experienced by 
the electron, 


(3 = —((3 Z sin (M , ~(3 Z cos (k 0 z ), - [f5 x sin (k 0 z) - (3 y cos (/c 0 ^)]) • (V.4) 

We look at the x component of equation (V.4) and integrate to hnd 


Px dt = 
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(V.5) 


If we assume perfect injection and helical motion, then v Xo = 0. Otherwise, the 
electron would wander away from the axis at constant velocity. The f3 y equation 
follows the same process with the same argument for the integration constant, giving 


Py 


-K 

7 


sin (k 0 z) , 


(V.6) 


where 


K = 


eB A 
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2nmc 


(V.7) 


Putting equations (V.5) and (V.6) into the z component of equation (V.4), 
we hnd that for perfect injection, 


P= 0 , (V.8) 

meaning that /3 Z is a constant, /3 0 . We now know z(t), and assuming that t — 0 
corresponds to the start of the undulator and we call that position z 0 = 0, then 


z(t) = C/dot. 

Let us now solve for the electron trajectories. We will solve the x component 
and state the y component when complete as they are solved in the same manner. 
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Starting with the x— c f3 x equation, 
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and similarly, 
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For perfect helical motion x 2 + y 2 = r 2 = a constant. This gives us 
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(V.9) 


(V.10) 


(V.ll) 


lk 0 p 0 

but this can only be a constant if the cross term in square brackets is equal to zero. 
The only way for this to be true is if Xo = 2/0 = 0. We now know the motion of the 
electrons to be 


, A) ct 

3. Resonance Condition 

One of the key characteristics of the FEL is its tunability. To demonstrate 
where that capability comes from, let us examine the resonance condition. The reso¬ 
nance condition is defined as being satisfied when one wavelength of light, A, passes 
over an electron as the electron traverses one wavelength of the undulator, Ao- This 
can be thought of as a photon-electron race, described by the following relation, 


(v.12) 


T(t) = {^k sia(kActh 
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(V.13) 


which says that in the same amount of time it takes the electron to move one Ao, the 
photon will have traveled A + Ao- If we rearrange this equation, we find 

a = a 0 (WA). (V.14) 
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We recall the Lorentz factor and the same assumptions from the previous 


derivations for /3 X and /3 y to find 


P 2 = (% + $} +ft 

K 2 

= — [sin 2 (k 0 z) + cos 2 (k 0 z)] + /3 2 

= — + Pi ■ (V.15) 

7 

We put this result into the definition of 7 and solve for (3 Z , 

2 1 

7 = 

7 2 

7 2 - K 2 - 72/32 
7 2 (l-/ 3 2 ) = l + K 2 

I 1 + A' 2 

P z = ^1-^-. (V.16) 

We now have (3 Z to put into equation (V.14). We see in equation (V.14) that 
/3 Z is in both the numerator and denominator. If we have a highly relativistic beam, 
then /3j ~ 1 and 7 > 1. We expand equation (V.16) in a Taylor series and substitute 
into equation (V.14) to find 


A = A, 


1 -Pz 

1 p z 


^ An 
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(V.17) 


From this relation, we see that to adjust the output light from the FEL, assuming 
the period of the undulator is fixed, we can either change the undulator parameter 
(move magnets closer/further away to increase/decrease B) or change the electron 
beam energy. In practice, individual FELs have demonstrated an ability to tune the 
output light over about an order of magnitude from their designed wavelength. 


4. Current Effects 

As the bunch charge passes through the undulator, it generates an optical wave 
due to the synchrotron radiation released as the electrons are accelerated transverse 
to their motion. The strength of the optical held is directly impacted by the amount 
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of current present in the beam. Without going through the derivation steps, we know 
that an electromagnetic wave equation, using the auxiliary held (and cgs units) can 
be written [50] 


1 d 


47T 




(V.18) 


This equation says that an electromagnetic wave is generated perpendicular to the 
direction of the motion of the charges in J. This bodes well for us as our undulator 
provides the necessary motion to obtain radiation along the axis of propagation. It 
can be shown that the sum over the electron bunch which constitutes J can be written 
as a dimensionless current, defined [50] 


8N(enKL) 2 p 

J =- 3—2 - 

7' 5 mc 


(V.19) 


where N is the number of periods in the undulator, L is the total length of the 
undulator, and p is the electron density. The source term of equation (V.18) can then 
be written in dimensionless terms as 


— J± oc (-je * c > , 
c x ' 


(V. 20) 
) is an ensemble 


where ( is the phase of an electron relative to the optical held and 
average over the electron phases. 

From equation (V.19), we see that if we desire to drive the optical held using 
the bunch charge, we can increase the source term by increasing the bunch charge or 
maintaining the charge, but increasing the charge density within the bunch. We do 
this by focusing the electron beam to a tight spot radially and compressing the beam 
longitudinally just before the bunch enters the undulator. When one also takes into 
consideration the gain of the FEL and can get to the high gain regime (j 3 > vr), then 
the gain of the FEL becomes [50] 


G(t) = 


t\P\ \ [i- 
e v V 2 


9 


(V.21) 


where r is the dimensionless time defined ct/L. In the case of a high power FEL, as is 
desired for military applications, the interaction region will necessarily be required to 
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have high gain as out-coupling of the laser light will be significant. It is this concept 
that drives our desire for a high bunch charge, low emittance gun for use with the 
NPS FEL. 
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VI. 


HERMITE-GAUSSIAN DECOMPOSITION 


Free electron lasers (FELs) generate electromagnetic radiation during the 
electron-wiggler interaction. This radiation then interacts with the electron beam 
to provide the feedback used in the FEL interaction. In oscillator FELs, the mirrors 
will interact with the optical mode which is fed back into the FEL interaction to 
be amplified. In an amplifier configuration, there is no feedback of radiation back 
into the interaction; instead, a seed laser or spontaneous emission provides an optical 
pulse and this radiation held is amplified. 

Knowing the spatial mode composition of the radiation held is important in 
determining the propagation properties of the output optical held. There are numer¬ 
ous measures of beam quality which describe how well a beam propagates compared 
with a pure Gaussian, such as M 2 measurements and Strehl ratios (see [62] and [63] 
for examples). We desire a suitable basis set of normal modes with which to decom¬ 
pose the radiation held of interest from an FEL. Since oscillator FELs typically use a 
spherical mirror cavity, they lend themselves to a Hermite-Gaussian (HG) basis set. 
We desire to use them for the ampliher case as well. This application of the HG basis 
set to the ampliher is not without problems, and our procedure for overcoming these 
shortcomings in order to determine useful information will be discussed. 

Through methods developed by Siegman [62], the HG modes are demonstrated 
as a suitable basis set of normal modes with which to decompose optical helds. Initial 
analysis of the more structured FEL oscillator case was explored and demonstrated by 
Vigil [64], We further demonstrate the analysis is also applicable to the ampliher FEL 
and develop a “least modal composition” optimization scheme for fixing the two free 
basis scaling parameters. This theory is then coded and merged by the author with 
the NPS FEL4D simulation code and applied to simulations of current and proposed 
FEL experiments [1]. 
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A. HG SOLUTIONS TO THE WAVE EQUATION 


Let us begin from first principles and derive a wave equation from Maxwell’s 
equations. We will then work through the reduced wave equation and into the paraxial 
case, which is applicable to FEL systems of interest. We will then show that the HG 
basis set is applicable to problems of this type and show how we can recast the wave 
equation for a preferred solution. 

1. The Electromagnetic Wave Equation 

Let us begin our wave equation search with Maxwell’s equations (see [20] or 
[19]), stated below in MKS units, 


V D 

= P ? 

(VI.1) 
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(VI.2) 

V B 

= o, 

(VI.3) 


d D 


V x H 

= J+ m • 

(VI.4) 


where D is the electric displacement held, p is the free charge density, E is the electric 
held, B is the magnetic held, H is the auxiliary magnetic held, and J is the free 
current density. 

We assume that our helds of interest exist in vacuum, so we can use the 
constitutive relations D = e$E and B = H where eo is the permittivity of free 
space and p 0 is the permeability of free space. We can obtain a wave equation 
from either the electric or magnetic held equations. We choose the electric held for 
illustration and note that the derivation using the magnetic held follows in a similar 
manner. 
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Assuming E and B are point-wise continuous in space and continuous in time, 
we begin by taking the curl of equation (VI.2), 
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(VI.5) 


where in the last step, we used v /pueo = 1/c. If we now conhne our discussions only 
to areas outside the interaction region, then there are no free charges and no currents. 
Thus, the entire right hand side becomes 0, leaving us with the homogeneous partial 
differential equation, 

Now that we have a wave equation, we assume the form 


E(x, y ,z,t) = E(x,y)e^ kz -^ , (VI.7) 

where k is the wavenumber, and put this into the wave equation, we generate a much 
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more complicated differential equation. Calculating term by term, we find 
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Putting these into the wave equation from equation (VI.6) results in 
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dx 2 dy 2 dz 2 

If we now apply the slowing varying phase and amplitude approximation, 
assuming dE/dz and dE/dt are small on the scale of an optical wavelength, then we 
can say 
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Applying this approximation and defining = (<9 2 / dx 2 + d 2 /dy 2 ) notation from 
before, we have a paraxial wave equation, 
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(VI.13) 
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We now change variables to a system that is commonly used in studying FEL 
interactions, altering the equation to “follow the light” by substituting u = z — ct. We 
also make dimensionless time by scaling r to the time required for light to propagate 
through the undulator (r = ct/L, where L is the length of the undulator). We now 
must change variables in the operator terms inside the parenthesis, 
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dz 
d_ 

Ft 

d Id 
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Our wave equation now looks like 

E = 0 , (VI.15) 

which contains all longitudinal effects in the r derivative term. 

Since we have made dimensionless time, let us make the transverse spatial 
terms dimensionless as well. We choose the scale factor k/2L so that 
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and our dimensionless wave equation becomes 
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(VI.16) 


2. HG Solutions 

Following the development of Siegman [62] and Vigil [64], we establish the 
HG polynomials as solutions to equation (VI. 16). Using a general approach that will 
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be used for the higher-order modes later, we first show that the zeroth order HG 
polynomial (which is a “pure” Gaussian) is a solution. Let us define a dimensionless 
held a(r, r) such that 

o(r,r) = a 0 e“( p(T)+1 ^M) , (VI.17) 


where a(x, y , r) oc E(x, y, r), r 2 = x 2 + y 2 , zq = Zq/L, Z 0 is the Rayleigh length, q{r) 
is a function describing the complex radius of curvature, and p(r) is a function that 
describes a distance dependent scaling factor. We stipulate the initial condition that 
the held is Gaussian at some time r = 0, implying p(0) = 0 and g(0) = 1. Plugging 
into equation (VI. 16), we find 
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Since f can vary across the entire image plane at all times r, then the only way this 
solution can be valid for all r and r is if both sides equal 0 simultaneously. This 
provides us with a system of two equations for solving for p and q. Starting with the 
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q equation, 
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We use this result in our solution to the p equation, 
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, since p 0 = p(0) = 0. (VI.20) 


It is apparent from the the fraction in the exponent of the field form (VI.17) 
and the solution (VI. 19) that we can define a dimensionless mode radius g(r) such 
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If we call p(0) = go = V^o, we can recast onr equation for q{r) in terms of g, 
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Making the appropriate substitutions into equation (VI.17), we find 
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which has the functional amplitude and phase form we desire. We can now define the 
phase 
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where the unclerbrace indicates the Gaussian amplitude we expected from equation 
(VI. 17) and the overbrace indicates the evolving phase term of the complex field. 


3. Higher Order HG Modes 

To account for higher order modes, we must adjust our general form for the 
dimensionless field, equation (VI. 17). Let us introduce two new functions, h and g, 
which control the transverse behavior of the field in their respective arguments and 
define 
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We apply the dimensionless wave equation (VI. 16) to this more general held, but take 
a useful detour to go through the appropriate partial derivatives, 
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and the h terms are derived in a similar fashion 
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Putting these terms into the dimensionless wave equation (VI. 16) and per¬ 
forming some algebra, 
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We now rearrange these terms into something a bit easier to analyze, 


v +4 


a = 0 = —— 


1 g" Ax g' 


Q 2 9 

V 


zoqg 9 


Aix 9 -\ 
9 Q 


only x dependence (= — a) 

1 h" _ Ay h! 
g 2 h ZqQQ h 


...h! q 

4*y-r~2 

h q z 


only y dependence (= —/3) 


Ai o 4 
Q A 


zo q- 


z 2 q 2 


only r dependence (=0) 

4 A 

— +4iP) . 
z 0 q 


no r dependence 


(VI.32) 


To have a solution to equation (VI.32) at some time r, the coefficient of r 2 
should be zero since it describes the field in the plane, and the remaining terms (“only 
x”, “only y”, and “no r dependence”) should sum to zero. We will let the x terms 
equal a constant, —a, and the y terms equal a different constant, —f3. We now solve 
these four differential equations starting with the terms proportional to r. 
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after factoring out and dividing both sides by 


0 

0, 


4 

zoq 2 ’ 

This is the same equation we solved for q in the zeroth order solution, so it has the 
same solution based on the same conditions, 
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(VI.33) 
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Next, let us focus on the x terms, 
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after substituting for q and 1?. 
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We recognize equation (VI.34) as Hermite’s equation[65]. 
(VI.34) in the explicit form of Hermite’s equation, we introduce 

xV2 
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and 
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d\ 2 

After substituting into (VI.34), we have the form we desire, 

<9 2 g „ dg ag 2 

W 2 ~ X dx + ^ 9 ~ ’ 

whose solutions are the Hermite polynomials with argument y: 
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where m e {JJ > 0} and J represents the set of integers. 


(VI.34) 


To recast equation 


(VI.35) 
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Similarly, the y terms from equation (VI.32) are solved to find 
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(VI.36) 


where we used x = x\/2/g, h! — \[2 dh/d\, h" = 2 d 2 h/dx 2 , and n G {J > 0}. 
Additionally, for Hermite solutions, a and fi from equations (VI.32), (VI.35) and 
(VI.36) are constrained by the Hermite equation form to 

a Q 2 n , fig 2 
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(VI.37) 


Moving to the last term in equation (VI.32), 
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We finally have the pieces we need to state the HG solution to our dimension¬ 
less wave equation that includes higher-order modes. Putting everything together, 
we have 


a(r,r) 
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(VI.39) 
(VIGO) 


where the Hermite polynomials through order rri — 7 are shown in Table 8 for refer¬ 
ence. 


Table 8. Hermite polynomials, H n , through 7th order, where n indicates the order. 


H 0 (x) = 1 

H\ (x) = 2x 

H 2 {x) = Ax 2 — 2 

H 3 (x) = 8x 3 — 12x 

H 4 {x) = 16a; 4 — 48a: 2 + 12 

H 5 (x) = 32a; 5 — 160a; 3 + 120a; 

H 6 (x) = 64a; 6 - 480a; 4 + 720a; 2 - 120 

H 7 (x) = 128a; 7 - 1344a; 5 + 3360a; 3 - 1680a; 


At this point, we have shown the general and complete solution to the wave 
equation. We are almost ready to begin applying it to useful problems. However, it 
would be nice if we could adjust the form so that it had a more intuitive feel. Define 
A = 7r q 2 and D = t/zq. These definitions will allow us to recast our solution in terms 
of two parameters, the mode area A and the mode phase D, for which we have an 
intuitive feel for how they should affect the field’s behavior. After these substitutions, 
our solutions to the wave equation become, 


a (r, D) 
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- Q ° H m (x\f^\ H n (yM) (VI.41) 
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The HG modes defined by the above equations form a complete, orthogonal 
basis [65], from which we can decompose and generate any optical field. We will 
normalize this equation in the next section. Figure 116 shows intensity patterns for 
a number of HG modes. Of note is that the order of the mode indicates the number 
of nulls in that coordinate direction. 
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Figure 116. Intensity patterns for all HG modes through third-order in each coordi¬ 
nate. The top left image is the 0,0 mode with the mode number increasing to the left 
and down. The bottom right image is the 3,3 mode. 
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B. FIELD DECOMPOSITION 


An arbitrary field can be decomposed with the HG basis set that we have 
shown to be a solution to the appropriate wave equation. This is done using an 
overlap integral to determine the complex coefficients, representing the “amount” of 
the (m, n ) HG mode present in the field. However, for the decomposition to be 
deterministic, the A and H scaling factors must be fixed. Since A and are related to 
p and q in our derivation, we have only recast them in more intuitive terms. We now 
address how to apply this decomposition method to representative fields and describe 
the process of fixing the scaling factors such that the minimum number of modes are 
required to describe an optical field of interest. 

Since the overall amplitude of the field is not of importance for the decompo¬ 
sition, we scale the field such that the integrated intensity (aa*) equals unity. This 
normalized field must then be projected against a set of normalized basis functions. 
In the previous section, we determined our basis functions, but we have not yet nor¬ 
malized them. To normalize equation (VI.41), we continue following the Vigil process 
[64] and start by determining the wavefront energy or power in the fundamental (0,0) 
mode, 

roc 

Pop = a o,o a o,o dx d y . (VI.43) 

J J —OO 

Putting equation (VI.41) into the integral, and substituting H 0 = 1, 



(VI.44) 


where we have made use of H 0 (x) = 1. Recalling that 



(VI.45) 
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we continue, setting b = 2tt/A, 


al / Fi F\ / [ 7r\ 

1 + Q 2 VV b) \\b) 
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6(l + fi 2 ) 

OgA 
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(VI.46) 


This is sufficient for the zeroth order HG polynomial, but does not address 
the higher-order modes. We now turn our attention there and recall that the HG 
polynomials comprise an orthonormal basis set with respect to a certain weighting 
function [65], 


Hi ( x) Hj (x) e x2 dx = Sij2^j\y/n , 


(VI.47) 


where is the Kronecker delta function. 


Proceeding for P m ,n as we did for P 0 ,o, 


P, 


m,n 



(VI.48) 


where we note that the integrals in equation (VI.48) are of the same form as equation 
(VI.47) with the exception of the coefficients. If we define x = x^2tt/A, then dx = 
a/A/2vt dx- 

Looking again at the first integral in equation (VI.48), and making note of the 
orthogonality and normalization relation of the HG polynomials, equation (VI.47), 
we solve 
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I = 



Q ~ 2-7T 

[H m (x)] 2 e A 





(VI.49) 


Using this result in equation (VI.48), we continue 



= P 0 ,o m!n!2 m 2 n . (VI.50) 


Thus, we can normalize the power in each higher order mode to the power in the (0,0) 
mode simply dividing by m\n\2 m 2 n . We can now define the mode normalized field as 
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Q"m,n 

Vm\n\2 m 2 n 


(VI.51) 


The next step is to determine the mode content of an unknown field. In 
our case, given an unknown field, a(r,Q), we know we can describe this field as 
a summation of numerous modes from our basis set where each mode has its own 
coefficient, C m ^ n . We can write 


OO 

a (r, fi) = C m,nA m , n (■ r , Q) . (VI.52) 

m,n =0 


To determine the mode content, we must determine the coefficients of interest. 
This is done through application of “Fourier’s trick” (a simple description can be 
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found in [66]). Start by multiplying both sides of equation (VI.52) by the complex 
conjugate of A m ^ n i and integrate over all space, 


aA* n , y da; d y = / / ^ C nhn A mtn A* m , y dx d y 


m,n 


f [ V'' n V l+o 2 JJ JJ ~2-Kr 2 / K 

. 1.1 x “ Vm\n\ 2 m 2" 


x 


m,n 


._ 1+ ^ 2 iJ m /iJ n /e~ 27rf2/A e~^ m '- n/ 

Vrn'\n'\2 m '2 n ' 


(VI.53) 


We pause to look at the imaginary exponential terms in equation (VI.53). Examining 
these terms, 

Q^m.n ,n' — g* — </> m ' >n ' ) 

we expand using equation (VI.42) and considering only the argument of the exponen¬ 
tial, 
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f VI.541 


Since we know the HG polynomials are orthogonal, equation (VI.53) shows that 
m! = m and n’ = n for a non-zero result then the imaginary exponent of equation 
(VI.54) is equal to unity. Using these observations, we continue solving for C m ^ n , so 
that 
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(VI.55) 


To calculate the C mjfl , one need only calculate the overlap integral on the right- 
hand side of equation (VI.55) and multiply by a coefficient that is dependent on the 
scaling factors chosen. In general, C' m ,n is complex. If the intensity of an optical held 
is E*E (oc a*a), which is a power per unit area and we have already normalized it so 
that the integral under the intensity curve is unity, then the percentage of energy or 
power in the (m, n) mode is described by n C m ^ n . 


C. 


CHARACTERIZING FEL FIELDS 


With our basis set properly developed and the mathematical framework of the 
process established, we can now describe the methodology for analyzing an optical 
held and then apply it to real and proposed FEL fields as calculated using the NPS 
FEL simulations [1]. For our purposes, we ignore the source of the held and assume 
that an optical held has been generated and is provided for our analysis. 

The hrst step is to normalize the optical intensity of the held. Since this is a 
simple scaling of the entire held, no modal decomposition information is lost, but it 
makes the computations and results much easier to understand and interpret. This 
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is accomplished in the usual manner, where a is the optical held provided, 


1 


[ba* (x, y)\ [ba (, x , y)} dx d y , 


(VI.56) 


where solving for b gives the scaling required to ensure a normalized optical intensity. 

The next step is to find the center of intensity. This accomplished by calcu¬ 
lating the first moment of the intensity (dropping the tildes for brevity), 


f x \a ( x , y) | 2 dx 
I |o(a;, 2 /)|- > da; 


(VI.57) 


and similarly for y. 

The second-order moment of the held, which will be useful momentarily, is 
found by [63] 


R 2 = 2 ff\ a ( x ,y)\ 2 [(® - xf + (y- yf] da; d y 


(VI.58) 


II l« (a;, 2/) r’ da; dy 

We have yet to really describe the arbitrary scaling factors 11 and A introduced 
in equation (VI.41), except to say that we wish to choose them such that the number 
of modes required to describe a held is minimized. Although we cannot state initially 
what the values of these scaling factors are without some manipulation, we can bound 
their values based on their definitions. 

We recall that and A are defined 


11 
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7T Q 


The curvature of the wavefront is tracked by 11 as it is related to the Gouy phase. 
Because the mode can be anywhere along its propagation path, there is no restriction 
on H’s value and it falls in the interval, — oo < H < oo. A negative value for H 
indicates the held is propagating to a waist; a positive value indicates it is propagating 
away from a waist. On the other hand, we know that the mode waist radius must be 
some positive non-zero number by definition, so A > 0. As an upper bound on A, we 
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consider the second-order moment of the held of interest. This is the 1/e 2 weighted 
radius of the intensity. We know that if the held of interest is composed only of the 
fundamental mode that the area of the mode using the 1/e point of the held is dehned 
by A and this would be the same as nR 2 . For the same R 2 , if there are higher-order 
modes contained in the held of interest, it implies that the scaling factor must be 
something less than R 2 as the higher-order modes tend to push the intensity away 
from the centroid, as seen in Figure 117. Thus, we can establish nR 2 as the upper 
bound for A, putting A in the interval 0 < A < nR 2 . 


-4 

x 10 

8 




Figure 117. Intensity patterns for two modes generated with the following dimen¬ 
sionless parameters: zo = 0.5, r = 0.625 or A = 4.0252, = 1.25. Intensities are 

normalized so that the integrated intensity = 1. Transverse axis labels are dimen¬ 
sionless. Intensity is in arbitrary units. The left image is the intensity pattern of the 
fundamental mode with the 1/e 2 radius, f = 1.1319. The right image is the intensity 
pattern of the 2-2 mode with the 1/e 2 radius, r = 2.5311. 


We must determine a discrimination method that can help us find the desired 
12 and A in order to minimize the number of modes necessary to describe the held. 
To compare choices of the scaling factors, we establish a “coefficient area” to describe 
the number of modes and the amount of content in those modes. To do this, we 
find the mode with the greatest content and normalize the \C m)V \ 2 for this mode to 
1. All the other \C m>n | 2 are scaled appropriately and then all \C m ^ are summed. 
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To illustrate this concept, suppose we had a held composed of 50% (0,0), 25% (2,0) 
and 25% (0,2). This mode would have a coefficient area of 2.0. We can then perform 
additional decompositions with varying scaling factors and compare their coefficient 
areas and choose the least coefficient area as “better” as fewer contributing modes are 
necessary. Another simple illustration of this process is to consider a mode made up 
of only the fundamental. We can perform the decomposition for the correct scaling 
factors and again for incorrect scaling factors, and see how higher-order modes are 
populated contributing to increased coefficient area as demonstrated in Figure 118. 

IA—| I A m , n | 



Figure 118. Modal composition histogram of a purely (0,0) mode (coefficient area = 
1.0) against the correct basis, = 1.25 and A = 4.0252 (left) and a different scaling 
factor, A = 3.0, resulting in a coefficient area = 1.06 and finding 94.7% in the (0,0) 
mode and ~2.5% in the (0,2) and (2,0) modes (right) with the combination of the 
remaining modes constituting less than 0.3%. Histograms are plotted as \C m ,n\ to 
better show the higher-order modes. 


We can use this coefficient area to build up a topology map showing the coef¬ 
ficient area found for many values of scaling factors. Figure 119 shows the topology 
when many scaling factors are examined for their resulting coefficient areas. Looking 
at this coefficient area map, we see that the behavior we expected: when the scaling 
factors are chosen carefully, there is a minimum value. Since we expect to be working 


220 


with fields that are primarily Gaussian as an assumption, we should expect to see 
coefficient area maps similar to Figure 119. The search area typically has a lower 
limit of about 10% of the nR 2 value estimated as the maximum A value discussed 
previously to minimize numerical effects that occur when the algorithm attempts to 
fit the field with a basis set with nearly zero cross-sectional area. 



- 9.4248 0.40252 


Figure 119. Topology of the coefficient area map for a purely (0,0) mode with A = 
4.0252 and G = 1.25. The scaling factors G and A are varied over their respective 
ranges, demonstrating that the improper choice results in a greater coefficient area 
value. 


Let us now see what happens when we generate coefficient area maps for modes 
that are nearly Gaussian to verify that there are no topological features that would 
give a simple minimization algorithm difficulties. We can see from Figure 120 that, 
in the nearly Gaussian case, there are no local minima that would cause difficulty 
for a simplistic minima-seeking algorithm. If the topology is complicated, as in the 
right hand subfigure of Figure 120, more robust or brute force search methodologies 
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must be used. In this more complicated case, the valley that rings the central peak 
would cause a minimum seeking algorithm to have to walk its way around once on 
the valley floor, greatly increasing computation time to achieve the minimum scaling 
factor solution. 



■ 943*6 040426 - 9.4248 0.77116 


Figure 120. The coefficient area topology map on the left is generated from the modes 
shown in the right hand part of Figure 118 containing 95% (0,0) mode and 2.5% (2,0) 
and (0,2) modes. The map on the right is generated for a field of interest consisting 
of 50% (0,0), 20% (0,2), 20% (2,0), and 10% (2,2) modes. 


Once we have minimized the number of modes required to represent a mode of 
interest, we must recognize that our algorithm, in actual practice, must truncate the 
number of modes used at some point. This truncation presents possible problems as 
we will discuss later when using this algorithm on actual FEL simulations. However, 
it is useful to have an independent method of measuring how well our coefficients 
represent the mode of interest. This can be done with a cross-correlation. 

The cross-correlation is a tool used frequently in Fourier analysis and can be 
used as a figure of merit. The cross-correlation can be thought of as a matched- 
filter: it looks through an unknown signal (or field in this case) for a known signal 
(the output optical field from the FEL simulation). In mathematical terms, the one 
dimensional cross-correlation (*), h , between two signals, / and g can be expressed 
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as 


■oo 


h{s) = f*g= f*(x)g(s + x) dx , (VI.59) 

J —OO 

where the asterisk indicates the complex conjugate and s is an offset variable allowing 
the signals to “shift” across one another. A graphical demonstration of the cross¬ 
correlation between two 3 by 3 arrays is shown in Figure 121 for illustrative purposes. 
With large arrays, this direct process can become computationally expensive. 
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Figure 121. This is a simple demonstration of how a two-dimensional cross-correlation 
of two 3x3 arrays is computed in the spatial domain. Starting from the top left, 
the test image, T. is sequentially overlapped with the image of interest, a, as shown 
in brown. The sum of the element products are placed in the corresponding element 
of the resulting 5x5 cross-correlation array, shown in blue. This process continues 
until the last element of the resulting cross-correlation array is filled as shown in the 
bottom right-most array. 


When large arrays are necessary, the Fourier convolution theorem offers us an 
alternative approach (see [65], for example). The convolution theorem states that the 
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Fourier transform, denoted as T { }, of a convolution in the spatial domain is the 
same as multiplication in the frequency domain. To illustrate this, let us examine a 
general convolution (*) for which we desire to know h, 


h = f*g 


h(s) = 


f(x)g(s — x ) da: 


(VI.60) 


now taking the Fourier transform of both sides, 
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H(k) = V2n(G (k) F (k)) , 


d u 


V2n J - c 


f(x)e~ zkx da: 


(VI.61) 


we thus show that a convolution in the spatial domain is multiplication in the fre¬ 
quency domain. To get back to the spatial domain, we only need to take the inverse 
Fourier transform of H(k ). While this is illustrative, our process involves the cor¬ 
relation, not the convolution. The primary difference from equation (VI.60) is that 
f(x) —> f*(x) and g(s — x) —> g(s + x). If we follow the same development as for the 
convolution, we prove the Wiener-Kinchin theorem and find that the cross-correlation 
in the spatial domain can be written 

h(s) — T~ x {H(k)} — T~ x (G(k) [F(/c)]*)| , (VI.62) 

where J-~ l is the inverse Fourier transform and * denotes the complex conjugate. 
Since Fourier transforms have been highly optimized in computational codes, it is 
this approach that is used to determine the cross-correlation between the test held 
and the original held of interest. 

Returning to our discussion of the cross-correlation in terms of the modal 
analysis algorithm, once the algorithm determines the correct scaling factors and the 
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coefficients, C m . n , are determined, we can use them to generate a normalized test 
field, T. If the normalized test held, T, is an exact copy of a, the cross-correlation 
becomes an auto-correlation which has a maximum value of 1.0. If T is only a “near” 
match to a, then the cross-correlation will have some maximum value less than but 
“near” 1.0. It can be used to indicate whether the decomposition process worked well. 
This two-dimensional cross-correlation generates a surface that describes how well the 
test held and original held match up. The location of greatest cross-correlation value 
corresponds to the location of best match. Figure 122 shows the typical surface that 
results after one performs a two-dimensional cross-correlation. 



Figure 122. A two-dimensional auto-correlation surface of a pure Gaussian mode 
where the center of intensity is at the origin. If the mode were offset, the cross¬ 
correlation peak would correspond to the center of intensity. The expanded window 
(double the mode window of 8 in dimensionless units) is a result of the shifting- 
theorem used in the cross-correlation. 
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In a similar fashion to the coefficient area maps, we can generate maps of cross¬ 
correlation values where the scaling factor space can be mapped and the maximum 
two-dimensional cross-correlation value for each test held generated can be captured. 
Figure 123 shows cross-correlation maps generated for our two example fields used 
previously to generate the coefficient area maps in Figure 120. It might seem that 
it would be more efficient to simply use the cross-correlation as the discrimination 
factor and maximize this value. Computationally the cross-correlation is much more 
time consuming. For a 200 by 200 complex arrray (held of interest), the 200 by 200 
point coefficient area maps in Figure 120 were completed in about 17 hours, or about 
1.5 seconds per scaling factor considered. By contrast, the cross-correlation maps in 
Figure 123 took 22 hours total, or about 33% longer. Additionally, if we compare 
Figures 120 and 123, the topology for finding an extremum appears to have more dis¬ 
cernible features in the coefficient area case. Thus, our decision to use the coefficient 
area as the primary discriminator for choosing appropriate scaling factors is validated 
on two accounts. The coefficient area approach has less of a computational cost than 
the cross-correlation and the topology found with the coefficient area map (of which, 
we have only shown two cases, but was generally noted from a much larger sample 
size) has more topological features allowing for better extremum discrimination. 

The approach is now fully described and we can turn our attention to holds 
generated using real FEL parameters. For these studies, the NPS 4D FEL amplifier 
simulation, described in [67], was used to simulate four FELs based on the best pa¬ 
rameters available, shown in Table 9. These particular FELs were chosen, primarily 
because the published data allowed for generation of the necessary input parameters 
for simulation. The Brookhaven National Laboratory (BNL) designs include a pro¬ 
posed untapered undulator design and moderate decreasing linear taper used in their 
Source Development Lab (SDL) and the Los Alamos National Laboratory (LANL) 
designs provide a gentler decreasing taper and a step taper for analysis. 
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Figure 123. The two cross-correlation maps are those for the same modes as shown in 
the coefficient area maps in Figure 120. The map on the left is for a nearly Gaussian 
field: 95 % (0,0) and 2.5% (2,0) and (0,2) modes. The map on the right is for a field 
containing 50% (0,0), 20% (0,2), 20% (2,0), and 10% (2,2) modes. 

The general algorithm described above was coded to perform the necessary 
calculations in a program separate from the NPS 4D FEL simulation code. The NPS 
FEL code provided the output of the optical field, real and imaginary parts, and 
necessary parameters for scaling for use as input into the analysis code. All optical 
analysis was performed just after the interaction region (at the end of the undulator) 
using up to the 7th order HG mode in each dimension (a total of 64 modes). The 
modal analysis results are tabulated in Table 10. 

For the four cases presented in Table 10, we find the amplifier FELs generate 
optical fields primarily in the fundamental mode. Using the minimum number of 
modes to describe the optical field, we see that the FELs are primarily composed of 
the fundamental mode, with cross-correlation values near 0.99, indicating excellent 
agreement between the simulation’s output data and a comparison field generated 
using the modes found from this analysis approach. The BNL SDL FEL is an inter¬ 
esting case as it was the only FEL that the analysis tool had problems decomposing. 
Regardless of the search starting point, the algorithm would find a minimum that 
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Table 9. FEL parameters used for HG decomposition. 



BNL 

BNL 

LANL 

LANL 


SDL [68] 

Proposed [69] 

Proposed [70] 

Proposed [70] 




step taper 

linear taper 

N 

256 

120 

110 

110 

^0 

3.9 cm 

3.275 cm 

2.18 cm 

2.18 cm 

K 

±v rms 

0.78 

0.7 

1.2 

1.2 

E 

102 MeV 

80 MeV 

81 MeV 

81 MeV 

q 

0.35 nC 

1.4 nC 

1 nC 

1 nC 

h 

A B 

B 

1 ps 
-5% 

2.8 ps 

0% 

1 ps 
-2.5% 

1 ps 
-8.0% 


had a cross-correlation value of around 0.73, indicating a poor match to the original 
field. It became apparent that the field of interest required additional higher-order 
modes to properly decompose the field. By doubling the number of modes in use in 
both dimensions (quadrupling the total number of modes), the algorithm was able to 
successfully find a set of scaling factors that generated a high cross-correlation value, 
which is reported in Table 10. 

If we examine the FEL simulation for this particular case, it is evident that 
there are possible reasons for needing more modes. Figures 124 and 125 show the 
relevant output from the FEL simulation. From Figure 124, we can see that there 
is significant structure surrounding the central lobe. This structure becomes signif¬ 
icantly more prominent as the mode is allowed to propagate further, although the 
mode content remains the same. In Figure 125, we see that the FEL reached satu¬ 
ration before even half the undulator was traversed as the gain stabilized at about 
r = 0.5. Saturation effects can drive power into higher-order modes in FELs at the 
expense of the fundamental mode. Thus, we were justified in extending the coeffi¬ 
cient series to higher-order modes, because if there had been no higher order mode 
content, there would have been no improvement in the the cross-correlation value 
for the minimum coefficient area solution. Had we chosen a basis set that assumed 
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Table 10. Modal analysis results of four amplifier FELs. The five greatest contributing 
modes (> 1%) and the percentage of their composition are shown. Cross-correlation 
values are reported as a figure of merit (FOM) for a field generated from the de¬ 
composed mode coefficients and corresponding scaling factors and A, with the 
simulation output field. (* Used 14 x 14 = 196 modes to perform decomposition.) 



Modes 

% Composition 

n 

A 

FOM 

BNL SDL* 

(0,0) 

32.3% 

1.053 

2.374 

0.983 


(2,0) 

15.0% 





(0,2) 

15.0% 





(0,4) 

6.4% 





(4,0) 

6.3% 




BNL 

(0,0) 

95.3% 

1.129 

1.138 

0.996 

proposed 

(2,2) 

1.0% 





(0,4) 

1.3% 





(4,0) 

1.3% 




LANL Step 

(0,0) 

99.4% 

3.675 

5.763 

0.991 

LANL 

(0,0) 

94.4% 

2.545 

3.748 

0.987 

Linear 

(0,2) 

1.5% 





(2,0) 

1.4% 





(0,4) 

1% 





(4,0) 

1% 





axisymmetric modes, such as the Laguerre-Gaussian polynomials, fewer modes may 
have been required to decompose this particular field of interest, but the ability to 
handle non-axisymmetric beams would have been lost. 

This methodology provides researchers and designers another method of mea¬ 
suring the quality of their FEL designs. By monitoring the mode content of the 
output optical field, the propagation characteristics can be predicted and optimized 
for the desired application. We have developed a robust and adaptable algorithm 
for decomposing arbitrary FEL optical fields as well as a well-known figure of merit 
for determining the success of the decomposition. Additionally, the decomposition 
method has provided confirmation that many amplifier FELs operate nearly in the 
fundamental mode regardless of their undulator design and that the mode content 
observed reflects well known FEL interaction phenomena such as saturation. 
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a (x,y, 1) 


a (x,y, 2) 



-5x5 -15 x 15 


Figure 124. The dimensionless optical field amplitude, \a(x,y)\. The image on the 
left is the field at the end of the undulator. The field on the right is the field after 
propagating another undulator length. Intensities are in arbitrary units and do not 
represent the same scale between pictures. 


ln(l+G(x)) 6.52 



Oxl 


Figure 125. The gain, G, experienced by the optical field as it traverses the undulator 
of the BNL SDL FEL. 
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VII. 


CONCLUSION 


In this work, we examined a free electron laser (FEL) from start to end. We 
started with an examination of some of the cathode types available for use in an 
electron source. Electron sources in the form of guns, boosters, and injectors were 
then discussed with simulation results presented for many different configurations and 
electron beam shapes. From simulation, we moved to discussion about the design and 
construction of the first superconducting radio frequency (SRF) gun in the United 
States. Experimental results showing the performance of the first prototype SRF 
gun were presented and compared with simulations showing that the gun is currently 
capable of generating 10% of the desired bunch charge at 25% design energy in only 
its first beam generation experiment. 

In the DC gun, the thermionic cathode has proven why it is the workhorse of 
electron beam generation in many applications. Unfortunately, the basic thermionic 
cathode’s applicability to short pulse regimes like those found in SRF guns is limited 
without additional beam line space being used for buncher cavities and compressor 
chicanes. If the desire is to decrease the overall FEL footprint, some other mechanism 
for pulsing the thermionic cathode, such as photogating, will be necessary, as has been 
shown to be practicable in our experiments. 

Bare metal cathodes have been proven in many other systems and demon¬ 
strated in the NPS SRF gun. They are, in principle, robust and do not require a 
great deal of “care and feeding” to maintain performance. There are preparation 
schemes for metal cathodes that promise an increase in quantum efficiency (QE) 
that could drive the gun bunch charge to levels approaching 1 nC per bunch. Alkali 
coated semiconductor cathodes currently in use in DC guns could also migrate into 
SRF guns. Unfortunately, not much is known about the effects of high QE materials, 
such as cesium, migrating onto the interior surfaces of an SRF cavity. SRF cavity 
cleanliness requirements are already extremely stringent, but as improved models of 
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the NPS SRF gun are built, previous editions could be used as test platforms for the 
cesiated cathodes without worry of impact to ongoing NPS experimental programs. 
Once the gun fields in the cathode region approach those necessary for turn-on of 
held emission arrays, the NPS SRF gun is an ideal platform for testing them in a 
realistic gun environment, especially with the ability to adjust the cathode position 
and the impinging electric fields almost at will. 

The DC gun received from Stanford was originally envisioned as the electron 
source for the NPS FEL and as a test bed for cathode studies. As a test bed, it 
has proven to be capable, but the full use of the device has been hindered by local 
radiation policy issues forcing a change in experimental direction that has impacted 
data collection from this useful tool. As an exercise in returning a significant piece 
of a decommissioned linear accelerator to service, the lessons learned in returning 
the DC gun to service and processing it back to near peak operating voltage were 
invaluable to the experimental group at NPS. It also serves as a case study for how 
rugged DC guns and thermionic cathodes are, even when stored in less than ideal 
conditions. 

Simulation studies of the NPS SRF cavity in the booster configuration showed 
that the booster design is successful. When coupled with the former Stanford DC 
gun, simulations indicate that transport of the beam through the booster cavity and 
into a merge can be accomplished, but would require beam focusing prior to entry 
into the booster, otherwise charge transport would suffer. Coupling the booster to 
the DC gun as an injector would be less than ideal due to the (i) extended beam line 
necessary to provide focusing to counter space charge effects while the beam is still 
at low energy, as well as (ii) the addition of a buncher cavity to compress the bunch 
sufficiently to correct for longitudinal beam growth. 

With the addition of a cathode and stalk to place the cathode, the NPS SRF 
booster cavity can be transformed into an electron source in its own right. The 
prototype cavity has demonstrated through simulation and experiment that it would 
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be a suitable gun. Transport acceptance is wide, allowing for nearly 180° of relative 
phase for which the beam can be fully transported out of the gun. Normalized 
emittances of 7 mm-mrad from 1 nC bunch charges at the design energy of 1.2 MeV 
have been simulated but await verification by experiment. The coaxial coupler has 
proven less problematic as experience with it has increased, but it still requires a 
significant drift space for insertion. The feared aperture impact of the coupler on 
the electron beam proved to be less of an issue than expected; however, it may be 
advantageous to consider other methods of coupling RF to the cavity that would 
allow for a more compact device. Experimentally, we verified the deleterious impact 
of the cathode stalk on the cavity fields, but we have not been able to quantify it as 
yet. It does indicate that a different cathode insertion system may be warranted. 

As part of the merger studies, a double quarter wave (or highly re-entrant 
pillbox) cavity is proposed as a booster. When coupled with the NPS SRF gun, it 
could constitute a > 3 MeV injector for the NPS FEL linear accelerator. The double 
quarter wave booster, when optimized for small spot size or best emittance, was capa¬ 
ble of delivering a suitable beam for injection into a linear accelerator in simulation. 
This booster design has the advantage of requiring only minimal development as the 
drawings and procedures for production would be based on a proven cavity design. 
To overcome problems with cooling of the superconducting solenoid in the NPS SRF 
cavity, the solenoids should be immersed in the liquid helium cryostat encompassing 
both the gun and booster. Both the long and short versions of the booster showed 
little emittance growth and were successful in transporting 1 nC bunches without 
issues. 

Shifting from experimental to theoretical studies, an analysis tool for FEL 
optical holds has been developed. The methodology is based upon Hermite-Gaussian 
decomposition, which has proven useful in the analysis of oscillator FELs. LInlike 
the oscillator, in which the basis set scaling factors can be determined from the os¬ 
cillator cavity, amplifier FELs have no unique factors that would fix the scaling for 
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decomposition. Under the assumption that the simplest answer is the best answer, a 
methodology for choosing scaling factors such that the number of Hermite-Gaussian 
polynomial terms necessary to decompose an optical held of interest was developed 
and implemented in the NPS FEL4D simulation code. Using the simulation code 
and parameters from the literature, the FEL optical fields of operating and proposed 
amplifiers were analyzed and demonstrated to operate primarily in the fundamen¬ 
tal mode. By using a decomposition methodology similar to those used by other 
laser types, it allows for a more judicious side-by-side comparison between FELs and 
traditional lasers when discussing propagation characteristics. 

In the future, as the NPS FEL4D code is tied to particle accelerator simulation 
codes and the NPS FEL is built, opportunities to examine an entire FEL system in 
simulation and experiment may become possible. It may soon be possible for a 
researcher to propose a system component, model it in a particle accelerator code 
such as GPT, perform FEL simulations, and predict the optical transport of the laser 
light and then perform the experiment for validation. Exciting times are in store for 
future students and researchers as the NPS FEL facilities and capabilities continue 
to grow. 
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